一种提取井间测距信号的方法及装置的制造方法

一种提取井间测距信号的方法及装置的制造方法
【专利摘要】本发明的实施方式提供了一种提取井间测距信号的方法及装置，该方法包括：去除原始信号中的直流分量；根据井间测距信号的频带确定上限数字频率和下限数字频率；在以下限数字频率和上限数字频率为边界的频域采样范围中，对去除直流分量后的原始信号进行离散傅里叶变换，得到频域信号；在得到的全部频域信号中搜索幅值最大的频域信号；计算幅值最大的频域信号对应的时域信号的频率和幅度；将该时域信号的频率和幅度确定为井间测距信号的频率和幅值。本发明具有计算速度快、频率分辨率高、计算结果精确度高等优点，为实时测量、精确计算井间距离提供有利条件，为丛式井等复杂结构井的钻井轨迹提供精确导向，能够有效防止邻井相碰事故发生。
【专利说明】
-种提取井间测距信号的方法及装置
技术领域
[0001] 本发明的实施方式设及电磁随钻测量技术领域，更具体地，本发明的实施方式设 及一种提取井间测距信号的方法及装置。
【背景技术】
[0002] 本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的 描述不因为包括在本部分中就承认是现有技术。
[0003] 丛式井W其经济高效、方便管理等诸多优势而被海洋、陆地油田所广泛采用。随着 大型油田的开发，丛式井定向钻井技术得到了长足的发展和提高，井距越来越小，加上井身 复杂结构，钻井过程中极易发生邻井井眼碰撞事故。井眼碰撞会带来严重的后果，轻则造成 套管变形，重则钻穿套管，影响油井的正常生产，造成巨大的经济损失。因此，解决丛式井的 井眼防碰问题，对保证丛式井钻井作业的安全具有极为重要的意义。
[0004] 丛式井正常钻进过程中，电磁随钻测量系统的磁场传感器会采集到一组交流信 号，该交流信号的频率与正钻井周围的已钻井数量等因素相关，可作为计算丛式井邻井距 离和方位的基础。在其他类型油井的随钻测井中，也能采集到类似的交流信号，可用于计算 正钻井与已钻井的距离和方位，本文将运种能用于计算邻井距离和方位的交流信号称为井 间测距信号。
[0005] 正常钻井过程中，电磁随钻测量系统的磁场传感器除了采集到井间测距信号外， 同时还会采集到地磁场的直流信号和噪声干扰信号。由于井间测距信号是计算邻井距离和 方位的关键，且非常微弱，可能为ICTinT数量级，因此在采集到的所有信号中快速、准确提取 出井间测距信号是邻井测距的基础。
[0006] 申请号为CN20121 Ol65021 .X的中国发明专利公开了一种基于传统的DFT (Discrete Fo urier Transform,离散傅里叶变换）的井间测距信号提取方案。该方案对电 磁随钻测量系统采集到的信号进行传统DFT处理，得到DFT幅度谱;通过对DFT幅度谱捜索得 到最大幅值和次大幅值，然后根据最大幅值及其对应的数字频率和次大幅值及其对应的数 字频率重构频域信号的数字频率和幅值初值;最后利用牛顿迭代公式，逼近井间测距信号 的频率和幅度。

【发明内容】

[0007] 在研究本发明的过程中
【申请人】发现，受转盘的旋转频率、正钻井周围的已钻井数 量等限制，用于计算邻井距离和方位的井间测距信号的频率往往固定于一特定范围内。
[000引例如，在丛式井正常钻进过程中，井下有动力钻具驱动钻头钻井，动力钻具本身会 有12化/min~18化/min的转速，而转盘转速一般都控制在2虹/min~10化/minW内。钻柱旋 转时，电磁随钻测量系统的磁场传感器采集到的井间测距信号是一组正弦信号，其信号周 期与转盘的旋转周期相关。转盘的旋转周期大约是0.6s~2.4s，当正钻井周围没有已钻井 时，井间测距信号的周期与转盘的旋转周期相同，也就是说，井间测距信号的频率大约为 0.4化~1.7Hz ;当正钻井周围有S 口已钻井时，井间测距信号的频率会提高S倍。而实际的丛 式井钻井中，同一口正钻井周围的IOm内一般不会超过20 口井，因此井间测距信号的频率最 大一般不会超过34化。根据W上研究发现可知，实际的丛式井钻井中，电磁随钻测量系统所 采集的井间测距信号是一种介于0.4Hz~34化的窄带正弦信号。
[0009] 此外，
【申请人】还发现，传统的DFT方法是将信号从时域变换到频域，是在整个数字 频率域0~1上进行等间隔均匀采样，但运种方法对于微弱窄带信号而言，会采集到很多无 用的频率点;并且，在数据量较小的条件下，采样间隔较大，也很难实现高的频率分辨率，进 而无法精确锁定有用信号频率，导致信号幅度检测的误差较大。例如申请号为 CN201210165021.X的中国发明专利所公开的井间测距信号提取方案，就是基于传统的DFT， 该方案存在计算结果精确度低、计算量较大等缺点。
[0010] 为了克服在电磁随钻测量中利用传统DFT方法提取井间测距信号的频率和幅度时 存在的精确度低、计算量大等问题，本发明提供一种提取井间测距信号的方法及装置。
[0011] 在本发明实施方式的第一方面中，提供了一种提取井间测距信号的方法，包括：
[0012] 获取电磁随钻测量系统采集的原始信号，去除所述原始信号中的直流分量；
[0013] 根据转盘的旋转频率的变化范围，W及正钻井周围一预设范围内的已钻井的最多 数量，确定井间测距信号的频带；
[0014] 根据所述井间测距信号的频带，W及预设的时域采样频率，确定上限数字频率和 下限数字频率；
[0015] 在W所述下限数字频率和所述上限数字频率为边界的频域采样范围中，对去除直 流分量后的原始信号进行离散傅里叶变换，将所述去除直流分量后的原始信号从时域转换 至频域，得到频域信号；
[0016] 在得到的全部频域信号中捜索幅值最大的频域信号；
[0017] 利用所述预设的时域采样频率，计算所述幅值最大的频域信号对应的时域信号的 频率和幅度；
[0018] 将该时域信号的频率和幅度确定为井间测距信号的频率和幅值。
[0019] 在本发明实施方式的第二方面中，提供了一种提取井间测距信号的装置，包括：
[0020] 获取模块，用于获取电磁随钻测量系统采集的原始信号；
[0021] 去直流模块，用于去除所述原始信号中的直流分量；
[0022] 频带确定模块，用于根据转盘的旋转频率的变化范围，W及正钻井周围一预设范 围内的已钻井的最多数量，确定井间测距信号的频带；
[0023] 边界频率确定模块，用于根据所述井间测距信号的频带，W及预设的时域采样频 率，确定上限数字频率和下限数字频率；
[0024] 第一变换模块，用于在W所述下限数字频率和所述上限数字频率为边界的频域采 样范围中，对去除直流分量后的原始信号进行离散傅里叶变换，将所述去除直流分量后的 原始信号从时域转换至频域，得到频域信号；
[0025] 捜索模块，用于在得到的全部频域信号中捜索幅值最大的频域信号；
[0026] 第二变换模块，用于利用所述预设的时域采样频率，计算所述幅值最大的频域信 号对应的时域信号的频率和幅度；
[0027] 结果确定模块，用于将该时域信号的频率和幅度确定为井间测距信号的频率和幅 值。
[0028] 借助于上述技术方案，本发明基于井间测距信号的窄带特点，提出W井间测距信 号的频带作为频域采样范围的离散傅里叶变换，即具有频率聚焦特性的离散傅里叶变换 (FFDFT)方法，从电磁随钻测量系统采集的原始信号中提取出井间测距信号。本发明可快 速、准确地在强干扰大噪声背景下提取出井间测距信号，具有计算速度快、频率分辨率高、 计算结果精确度高等优点，为实时测量、精确计算井间距离提供有利条件，为丛式井等复杂 结构井的钻井轨迹提供精确导向，能够有效防止邻井相碰事故发生。
【附图说明】
[0029] 通过参考附图阅读下文的详细描述，本发明示例性实施方式的上述W及其他目 的、特征和优点将变得易于理解。在附图中，W示例性而非限制性的方式示出了本发明的若 干实施方式，其中：
[0030] 图1为本发明提供的提取井间测距信号的方法的流程示意图；
[0031] 图2为实施例S中分别利用传统的DFT方法和FFDFT方法对计算机产生的纯净正弦 信号X的处理结果；
[0032] 图3为实施例S中分别利用传统的DFT方法和FFDFT方法对计算机产生的纯净正弦 信号Y的处理结果；
[0033] 图4为实施例S中分别利用传统的DFT方法和FFDFT方法对计算机产生的纯净正弦 信号Z的处理结果；
[0034] 图5为实施例S中分别利用传统的DFT方法和F抑FT方法对计算机产生的含白噪声 的正弦信号X'的处理结果；
[0035] 图6为实施例S中分别利用传统的DFT方法和FFDFT方法对计算机产生的含白噪声 的正弦信号Y'的处理结果；
[0036] 图7为实施例S中分别利用传统的DFT方法和F抑FT方法对计算机产生的含白噪声 的正弦信号Z'的处理结果；
[0037] 图8为本发明提供的提取井间测距信号的装置的结构框图；
[0038] 在附图中，相同或对应的标号表示相同或对应的部分。
【具体实施方式】
[0039] 下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出运 些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非W任何 方式限制本发明的范围。相反，提供运些实施方式是为了使本公开更加透彻和完整，并且能 够将本公开的范围完整地传达给本领域的技术人员。
[0040] 本领域技术技术人员知道，本发明的实施方式可W实现为一种系统、装置、设备、 方法或计算机程序产品。因此，本公开可W具体实现为W下形式，即：完全的硬件、完全的软 件(包括固件、驻留软件、微代码等），或者硬件和软件结合的形式。
[0041] 根据本发明的实施方式，提出了一种提取井间测距信号的方法及装置。
[0042] 在本文中，需要理解的是，所设及的术语"原始信号"是指钻井过程中电磁随钻测 量系统采集的、未经处理的信号;所设及的术语"井间测距信号"是指从"原始信号"中提取 出的能够用于计算邻井距离和方位的交流信号;所设及的术语"具有频率聚焦特性的离散 傅里叶变换"，即F抑FT(Frequency Fo州sing Discrete Fourier Transform),是指采样范 围限制于一预设频带的DFT。
[0043] 此外，附图中的任何元素数量均用于示例而非限制，W及任何命名都仅用于区分， 而不具有任何限制含义。
[0044] 示例性方法
[0045] 本发明提供一种提取井间测距信号的方法，如图1所示，该方法包括：
[0046] 步骤SI,获取电磁随钻测量系统采集的原始信号，去除所述原始信号中的直流分 量；
[0047] 步骤S2,根据转盘的旋转频率的变化范围，W及正钻井周围一预设范围内的已钻 井的最多数量，确定井间测距信号的频带；
[0048] 步骤S3,根据所述井间测距信号的频带，W及预设的时域采样频率，确定上限数字 频率和下限数字频率；
[0049] 步骤S4,在W所述下限数字频率和所述上限数字频率为边界的频域采样范围中， 对去除直流分量后的原始信号进行离散傅里叶变换，将所述去除直流分量后的原始信号从 时域转换至频域，得到频域信号；
[0050] 步骤S5,在得到的全部频域信号中捜索幅值最大的频域信号；
[0051] 步骤S6,利用所述预设的时域采样频率，计算所述幅值最大的频域信号对应的时 域信号的频率和幅度；
[0052] 步骤S7,将该时域信号的频率和幅度确定为井间测距信号的频率和幅值。
[0053] W下对步骤Sl~S7分别进行详细介绍。
[0054] 步骤SI,获取电磁随钻测量系统采集的原始信号，去除原始信号中的直流分量。
[0055] 由于电磁随钻测量系统采集的原始信号中不仅包含井间测距信号，还会包含地磁 场直流信号和噪声干扰信号等，为了提取井间测距信号，本步骤首先利用去直流方法去除 原始信号中包含的地磁场直流信号。
[0056] 具体实施时，该步骤可W采用常规的去直流方法。
[0057] 例如，xo(n)表示通过时域采样得到的原始信号，设采集的数据长度（即采集的数 据点数)为N个，其中n = 0，l，2…，N-1。去直流的过程为：
[0化引 m苗先管断右原始信号的平均值
[0化9]
[0060] (2)然后从米集到的原始信号XO (n)中去除之
[0061 ]玉-X
[0062] 其中，x(n)表示去除直流分量后的原始信号。
[0063] 步骤S2,根据转盘的旋转频率的变化范围，W及正钻井周围一预设范围内的已钻 井的最多数量，确定井间测距信号的频带。
[0064] 本步骤中，正钻井周围的预设范围可W是W正钻井为圆屯、、半径为一预设值(例如 IOm或20m)的圆区域。
[0065] 令fha表示井间测距信号的最大频率;fla表示井间测距信号的最小频率;fhz表示转 盘的最大旋转频率，flZ表示转盘的最小旋转频率;S表示正钻井周围预设范围内的已钻井的 最多数量，则有：
[0066] fha = fhz X S
[0067] fla = flz
[0068] 例如，在丛式井的钻井实况中，转盘的旋转周期介于0.6s~2.4s，相应的，转盘的 旋转频率的变化范围为0.4Hz~1.7Hz，而W每口正钻井为圆屯、、半径为1 Om的圆形区域内最 多有20 口已钻井，因此，井间测距信号的频带为0.4化~34化。
[0069] 步骤S3,根据井间测距信号的频带，W及预设的时域采样频率，确定上限数字频率 和下限数字频率。
[0070] 本步骤中的时域采样频率是指对井间测距信号进行时域采样的频率，根据时域采 样定理，该时域采样频率应大于或等于井间测距信号的最大频率的2倍(一般地，可选为5~ 10倍)。
[0071] 令井间测距信号的频带为fla~fha，则井间测距信号的最大频率为fha，最小频率为 fla。
[0072] 令时域采样频率为fs，上限数字频率为fhd，下限数字频率为fid，则有：
[0073]
[0074]
[0075] 通过研巧可知，井间测距信号受转盘的旋转频率、正钻井周围的已钻井数量等因 素影响，属于窄带信号，例如丛式井的井间测距信号是频带为0.4Hz~％化的窄带正弦信 号。考虑到井间测距信号的窄带特点，本发明可W在井间测距信号的频带上进行离散傅里 叶变换（即具有频率聚焦特性的离散傅里叶变换FFDFT), W达到提高频率分辨率、避免采集 无用的频率点，提高计算速度和精确度的效果。
[0076] 步骤S4,通过对去除直流分量后的原始信号进行离散傅里叶变换，将所述去除直 流分量后的原始信号从时域转换至频域，得到频域信号，其中，本步骤进行的离散傅里叶变 换所采用的频域采样范围是W所述下限数字频率和所述上限数字频率为边界，即在所述下 限数字频率和所述上限数字频率之间（包含边界点）进行具有频率聚焦特性的离散傅里叶 变换(FFDFT)。
[0077] 根据前述介绍，上限数字频率为fhd，下限数字频率为fid，则，本步骤进行的FFDFT 处理，是在下限数字频率fid和上限数字频率fhd之间进行等间隔均匀采集N个点，N的取值根 据计算的时效性和计算结果的精度要求来确定。
[0081 ] k = 0,l, . . . ,M-I
[0078] 本步骤中FFDFT处理将去除直流分量后的原始信号从时域转换至频域，得到频域 信号，可采用如下公式表示：
[0079]
[0080]
[0082] 其中，Xc(fck)表示频域信号;f。读示频域信号Xc(fck)对应的数字频率;x(n)表示去 除直流分量后的原始信号;n表示时域采样点的编号;k表示频域采样点的编号;N表示时域 采样点总数;M表示频域采样点总数。
[0083] 步骤S5,在得到的全部频域信号中捜索幅值最大的频域信号。
[0084] 具体的，本步骤可W按照如下过程执行：
[0085] 步骤S51，分别W频域信号的数字频率和幅值为坐标，标注步骤S4得到的全部频域 信号，得到幅度谱。
[0086] 步骤S52,在幅度谱中捜索幅值最大的频域信号。
[0087] 步骤S6,利用所述预设的时域采样频率，计算所述幅值最大的频域信号对应的时 域信号的频率和幅度。
[008引令/。。表示时域信号的频率，身表示时域信号的幅值，IXc(fck) Imax表示步骤S5中捜 索到的幅值最大的频域信号的幅值，fsk-max表示幅值最大的频域信号的数字频率，则有：
[0089]
[0090]
[0091] 步骤S7,将步骤S6计算得到的时域信号的频率和幅度确定为井间测距信号的频率 和幅值，完成提取井间测距信号。
[0092] 应当注意，尽管在附图中W特定顺序描述了本发明方法的操作，但是，运并非要求 或者暗示必须按照该特定顺序来执行运些操作，或是必须执行全部所示的操作才能实现期 望的结果。附加地或备选地，可W省略某些步骤，将多个步骤合并为一个步骤执行，和/或将 一个步骤分解为多个步骤执行。
[0093] 实施例一
[0094] 本实施例基于丛式井的实际钻进情况，对本发明提供的提取井间测距信号的方法 进行介绍。
[0095] 步骤1，获取电磁随钻测量系统采集的原始信号，去除原始信号中的直流分量。
[0096] 步骤2,根据转盘的旋转频率的变化范围，W及正钻井周围预设范围内的已钻井的 最多数量，确定井间测距信号的频带。
[0097] 在丛式井的实际钻进中，转盘的旋转周期的变化范围是0.6s~2.4s，转盘的旋转 频率的变化范围为0.4Hz~1.7Hz，而每口正钻井周围的IOm内一般最多有20口已钻井，因 此，井间测距信号的频带为0.4Hz~％化，井间测距信号的最大频率为％化，最小频率为 0.4Hz。
[009引为了方便计算，本实施例将井间测距信号的频带设置为0~40化，即化3 = 40化，fia = 0Hz。
[0099] 步骤3,根据井间测距信号的频带，W及预设的时域采样频率，确定上限数字频率 和下限数字频率。
[0100] 本实施例中，时域采样频率为fs = 200Hz，则上限数字频率呆
^限数 字频率为
[0101] 巧W4,仕h限甄字频率0和上限数字频率0.2之间，对去除直流分量后的原始信号 进行F抑FT处理，将去除直流分量后的原始信号从时域转换至频域，得到频域信号：
[0102]
[0103]
[0104]
[0105] 其中，姑(fck)表示频域信号;fck表示频域信号Xc(fck)对应的数字频率;x(n)表示去 除直流分量后的原始信号;n表示时域采样点的编号;k表示频域采样点的编号;N表示时域 采样点总数，M表示频域采样点总数。
[0106] 相邻的两个频域采样点的数字频率的间隔为
[0107] 步骤5,在得到的全部频域信号中捜索幅值最大的频域信号。
[0108] 步骤6,计算幅值最大的频域信号对应的时域信号的频率和幅度，将计算出的时域 信号的频率和幅度确定为井间测距信号的频率和幅值。
[0109] 令i；。表示时域信号的频率，1表示时域信号的幅值，|Xc(fd〇 kx表示捜索到的幅 值最大的频域信号的幅值，fsk-max表示幅值最大的频域信号的数字频率，则有：
[0110]
[0111]
[0112] 雨巧桃明的足，伞及明隊了可W应用于丛式井之外，还可W适用于其他类型的复 杂结构井，例如稠油开采中的SAGD双水平井，煤层气开采中的联通井，U型井等等。本实施例 仅仅是本发明的一种具体实施例而已，并不用于限定本发明的保护范围。
[0113] 实施例二
[0114] 本实施例分别利用传统DFT方法、W及FFDFT方法提取丛式井的井间测距信号，并 对二者进行对比。传统DFT方法中，频域采样点总数M和时域采样点总数N必须一致，即M=N; 而本发明的FFDFT方法中并不要求两者必须一致。但在本实施例中，为了对两种方法进行比 较，令FFDFT方法中的M=N。
[0115] (1)利用传统DFT方法提取丛式井的井间测距信号
[0116] 设原始信号x(n)，n = 0，l，...，N-l。传统的DFT处理是对x(n)的傅里叶变换在数字 频率域0~1上进行N等分均匀采样，其计算公式如下：
[0117]
[011 引
[0119]其中，X化读示经过传统的DFT处理得到的频域信号，f读示X(fk)的数字频率。相 邻的两个频域采样点的数字频率的间隔为
[0120] 对于数字频率为时，幅值为A，数据长度为N的正弦信号，可表示为:x(n)=Asin(23i fon)，n = 0,1,. . . ,N-1。
[0121] 其传统DFT的幅度谱为：
[0122]
[0123]
[0124]
[0125] 在数字频率域0~1的N个频域采样点中，如果有某个频域采样点的数字频率fk恰 好等于时，则通过捜索DFT幅度谱，即可得到最大峰值及其对应的数字频率fo，进而确定正弦 信号的幅度和频率。
[0126] 但在实际情况中，在数字频率域0~1的N个频域采样点中，通常是没有某一个频域 采样点的数字频率fk恰好等于时，于是通过捜索最大的幅值来确定正弦信号的幅度和频率 就会出现误差，其最大误差点将出现在数字频^
,该处的幅值为：
[0127]
[01 %]则该幅值与理论最大幅值的相对误差为：
[0129]
[0130] (2)利用FFDFT方法提取丛式井的井间巧鹏信号
[0131] 参考实施例一，FFDFT方法是在数字频率域0~0.2上进行N等分均匀采样，其中数 字频率域0~0.2是依据丛式井的井间测距信号的频带确定的。相比于传统的DFT方法， FFDFT方法将采样的频率范围从数字频率域0~1压缩到了0~0.2,即本发明的频域采样范 围聚焦于井间测距信号的频带上，可有效避免采集无用点的情况，提高计算精度。
[0132] FFDFT力巧的计貸公式为：
[013引 k = 0，l，. . .，N-1
[013
[013
[0136] FFDFT方法中相邻的两个频域采样点的数字频率间隔女 传统DFT方法中相邻的两个频域采样点的数字频率间隔3
f比可知，本发明显著 提高了频率分辨率。
[0137]考虑误差情况，FFDFT方法的最大误差点出现在数字频马
处，该处的幅值为：
[013 引
[0139] 则该幅值与理论最大幅值的相对误差为：
[0140]
[0141] 相比于传统的DFT方法，FFDFT方法中幅值提取的最大相对误差由36%骤降至 1.6%，数字频率提取的最大误差也从也就是说，本发明的频率定位的精 度比传统DFT方法的频率定位精度提高了 4倍。
[0142] 实施例S
[0143] 本实施例分别利用传统的DFT方法、W及利用FFDFT方法（即采样范围限制于一预 设频带的DFT)，对计算机产生的纯净正弦信号(不含噪声）、计算机产生的含有随机噪声的 正弦信号进行处理。传统DFT方法中，频域采样点总数M和时域采样点总数N必须一致，即M = N;而本发明的F抑FT方法中并不要求两者必须一致。但在本实施例中，为了对两种方法进行 比较，令F抑FT方法中的M=N。处理结果如图2-图7所示(横轴表示频率，单位为化，纵轴表示 幅值，单位为V)。
[0144] 图2-图4是本实例分别利用传统的DFT方法和FFDFT方法对计算机产生的纯净正弦 信号的处理结果。
[0145] 图2-图4中的（a)分别为计算机产生的纯净的不含噪声的正弦信号X、Y、Z的示意 图；其中，正弦信号X的幅度为3V，频率为甜Z ;正弦信号Y的幅度为2V，频率为細Z ;正弦信号Z 的幅度为IV，频率为IHz。
[0146] 图2-图4中的（b)分别为利用传统DFT方法对相应的正弦信号X、Y、Z进行处理得到 的幅度谱示意图。
[0147] 图2-图4中的（C)分别为利用FFDFT方法对相应的正弦信号X、Y、Z进行处理得到的 幅度谱示意图。
[0148] 图2-图4中的（d)分别为利用传统DFT方法进行信号消噪后恢复的正弦信号X、Y、Z 的时域波形示意图。
[0149] 图2-图4中的（e)分别为利用F抑FT方法进行信号消噪后恢复的正弦信号X、Y、Z的 时域波形示意图。
[0150] 本实施例应用FFDFT方法和传统的DFT方法对信号的幅度和频率及其相对误差的 检测对比结果列于表1中。由表1数据可见，FFDFT方法的幅度检测相对误差在1 % W内，频率 检测相对误差在0.1% W内，而传统的DFT方法的幅度检测相对误差高达34.553%，频率检 测相对误差达到2.340%，应用FFDFT方法的检测精度明显优于DFT方法，验证了本发明方法 对正弦信号幅度和频率检测的准确性。
[0151] 表1
[0152]
[0153] 图5-图7是本实例分别利用传统的DFT方法和FFDFT方法方法对计算机产生的含有 随机噪声的正弦信号的处理结果。
[0154] 图5-图7中的(a)分别为计算机产生的含白噪声的正弦信号X'、Y'、Z'的示意图；其 中，正弦信号X'的幅度为3V，频率为5Hz;正弦信号Y'的幅度为2V，频率为細Z;正弦信号Z'的 幅度为IV，频率为IHz。
[0155] 图5-图7中的(b)分别为利用传统DFT方法对相应的正弦信号X'、Y'、Z'进行处理得 到的幅度谱示意图。
[0156] 图5-图7中的k)分别为利用F抑FT方法对相应的正弦信号X'、Y'、Z'进行处理得到 的幅度谱示意图。
[0157] 图5-图7中的(d)分别为利用传统DFT方法进行信号消噪后恢复的正弦信号X'、Y'、 Z'的时域波形示意图。
[015引图5-图7中的（e)分别为利用F抑FT方法进行信号消噪后恢复的正弦信号X'、Y'、Z' 的时域波形示意图。
[0159] 本实施例应用FFDFT方法和传统的DFT方法对信号的幅度和频率及其相对误差的 检测对比结果列于表2中。由表2数据可见，FFDFT方法的幅度检测相对误差在6% W内，频率 检测相对误差在1.6% W内，而传统的DFT方法的幅度检测相对误差高达37%，频率检测相 对误差达到2.340% ,FFDFT方法的检测精度明显优于DFT法，验证了本发明方法对加噪正弦 信号的幅度和频率提取的准确性。
[0160] 表2
[0161]
[0162]综上所述，本发明基于井间测距信号的窄带特点，提出W井间测距信号的频带作 为采样范围的离散傅里叶变换，即FFDFT方法，从电磁随钻测量系统采集的原始信号中提取 出井间测距信号。本发明可快速、准确地在强干扰大噪声背景下提取出井间测距信号，具有 计算速度快、频率分辨率高、计算结果精确度高等优点，为实时测量、精确计算井间距离提 供有利条件，为丛式井等复杂结构井的钻井轨迹提供精确导向，能够有效防止邻井相碰事 故发生。
[01创示例性装置
[0164] 在介绍了本发明示例性方法之后，接下来，参考图8对本发明示例性装置进行介 绍。
[0165] 如图8所示，本发明提供一种提取井间测距信号的装置，包括：
[0166] 获取模块81，用于获取电磁随钻测量系统采集的原始信号；
[0167] 去直流模块82，用于去除所述原始信号中的直流分量；
[0168] 频带确定模块83,用于根据转盘的旋转频率的变化范围，W及正钻井周围一预设 范围内的已钻井的最多数量，确定井间测距信号的频带；
[0169] 边界频率确定模块84,用于根据所述井间测距信号的频带，W及预设的时域采样 频率，确定上限数字频率和下限数字频率；
[0170] 第一变换模块85,用于在W下限数字频率和上限数字频率为边界的频域采样范围 中，对去除直流分量后的原始信号进行离散傅里叶变换，将所述去除直流分量后的原始信 号从时域转换至频域，得到频域信号；
[0171] 捜索模块86,用于在得到的全部频域信号中捜索幅值最大的频域信号；
[0172] 第二变换模块87,用于利用所述预设的时域采样频率，计算所述幅值最大的频域 信号对应的时域信号的频率和幅度；
[0173] 结果确定模块88,用于将该时域信号的频率和幅度确定为井间测距信号的频率和 幅值。
[0174] 可选地，所述的频带确定模块83采用如下公式确定井间测距信号的频带：
[0175] fha = fhz X S
[0176] fla = flz
[0177] 其中，fha表示井间测距信号的最大频率;fla表示井间测距信号的最小频率;fhz表 示转盘的最大旋转频率，f Iz表示转盘的最小旋转频率;S表示正钻井周围所述预设范围内的 已钻井的最多数量。
[0178] 可选地，所述的边界频率确定模块84采用如下公式确定上限数字频率和下限数字 频率：
[0179]
[0180]
[0181] 其中，f S表示预设的时域采样频率;fhd表示上限数字频率;fid表示下限数字频率。
[0182] 可选地，所述的第一变换模块85采用如下公式将所述去除直流分量后的原始信号
从时域车电疏巧脯+击化但否Il脯诚隹具.
[0183]
[0184]
[01 化]k = 0，l, . . . ,M-I
[0186] 其中，姑(fck)表示频域信号;f。读示频域信号Xc(fck)对应的数字频率;x(n)表示去 除直流分量后的原始信号;n表示时域采样点的编号;k表示频域采样点的编号;N表示时域 采样点总数;M表示频域采样点总数。
[0187] 可选地，所述的捜索模块86包括：
[0188] 标注子模块，用于分别W频域信号的数字频率和幅值为坐标标注得到的全部频域 信号，得到幅度谱；
[0189] 捜索子模块，用于在幅度谱中捜索幅值最大的频域信号。
[0190] 可选地，所述的第二变换模块87采用如下公式计算所述幅值最大的频域信号对应 的时域倍耳的麻莖巧嶋底.
[0191]
[0192]
[0193] 其中，/。。表示时域信号的频率;J表示时域信号的幅值；IXc(fck) Ux表示所述幅值 最大的频域信号的幅值;fsk-max表示所述幅值最大的频域信号的数字频率。
[0194] 图8所示的提取井间测距信号的装置与图1所示的提取井间测距信号的方法基于 相同的发明思想实现，其【具体实施方式】可参照前述对提取井间测距信号的方法的介绍，此 处不再寶述。
[0195] 应当注意，尽管在上文详细描述中提及了提取井间测距信号的装置的若干模块或 子模块，但是运种划分仅仅并非强制性的。实际上，根据本发明的实施方式，上文描述的两 个或更多模块的特征和功能可W在一个模块中具体化。反之，上文描述的一个模块的特征 和功能可W进一步划分为由多个模块来具体化。
[0196] 虽然已经参考若干【具体实施方式】描述了本发明的精神和原理，但是应该理解，本 发明并不限于所公开的【具体实施方式】，对各方面的划分也不意味着运些方面中的特征不能 组合W进行受益，运种划分仅是为了表述的方便。本发明旨在涵盖所附权利要求的精神和 范围内所包括的各种修改和等同布置。
[0197] W上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，W上所述仅为本发明的具体实施例而已，并不用于限定本发明的保 护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。
[0198] 本领域技术人员还可W 了解到本发明实施例列出的各种说明性逻辑块 (illustrative logical block),单元，和步骤可W通过电子硬件、电脑软件，或两者的结 合进行实现。为清楚展示硬件和软件的可替换性QnterchangeabiIity)，上述的各种说明 性部件（illustrative components)，单元和步骤已经通用地描述了它们的功能。运样的功 能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员 可W对于每种特定的应用，可W使用各种方法实现所述的功能，但运种实现不应被理解为 超出本发明实施例保护的范围。
[0199] 本发明实施例中所描述的各种说明性的逻辑块，或单元，或装置都可W通过通用 处理器，数字信号处理器，专用集成电路(ASIC)，现场可编程口阵列或其它可编程逻辑装 置，离散口或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功 能。通用处理器可W为微处理器，可选地，该通用处理器也可W为任何传统的处理器、控制 器、微控制器或状态机。处理器也可W通过计算装置的组合来实现，例如数字信号处理器和 微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类 似的配置来实现。
[0200] 本发明实施例中所描述的方法或算法的步骤可W直接嵌入硬件、处理器执行的软 件模块、或者运两者的结合。软件模块可W存储于RAM存储器、闪存、ROM存储器、EPROM存储 器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介 中。示例性地，存储媒介可W与处理器连接，W使得处理器可W从存储媒介中读取信息，并 可W向存储媒介存写信息。可选地，存储媒介还可W集成到处理器中。处理器和存储媒介可 W设置于ASIC中，ASIC可W设置于用户终端中。可选地，处理器和存储媒介也可W设置于用 户终端中的不同的部件中。
[0201] 在一个或多个示例性的设计中，本发明实施例所描述的上述功能可W在硬件、软 件、固件或运=者的任意组合来实现。如果在软件中实现，运些功能可W存储与电脑可读的 媒介上，或W-个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可W 是任何通用或特殊电脑可W接入访问的可用媒体。例如，运样的电脑可读媒体可W包括但 不限于341、1?01、66?1?01八0-1?01或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任 何可W用于承载或存储W指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理 器读取形式的程序代码的媒介。此外，任何连接都可W被适当地定义为电脑可读媒介，例 如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双 绞线、数字用户线(D化)或W例如红外、无线和微波等无线方式传输的也被包含在所定义的 电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、错射盘、光盘、DVD、软盘 和蓝光光盘，磁盘通常W磁性复制数据，而碟片通常W激光进行光学复制数据。上述的组合 也可W包含在电脑可读媒介中。
【主权项】
1. 一种提取井间测距信号的方法，其特征在于，包括： 获取电磁随钻测量系统采集的原始信号，去除所述原始信号中的直流分量； 根据转盘的旋转频率的变化范围，W及正钻井周围一预设范围内的已钻井的最多数 量，确定井间测距信号的频带； 根据所述井间测距信号的频带，W及预设的时域采样频率，确定上限数字频率和下限 数字频率； 在W所述下限数字频率和所述上限数字频率为边界的频域采样范围中，对去除直流分 量后的原始信号进行离散傅里叶变换，将所述去除直流分量后的原始信号从时域转换至频 域，得到频域信号； 在得到的全部频域信号中捜索幅值最大的频域信号； 利用所述预设的时域采样频率，计算所述幅值最大的频域信号对应的时域信号的频率 和幅度； 将该时域信号的频率和幅度确定为井间测距信号的频率和幅值。2. 根据权利要求1所述的提取井间测距信号的方法，其特征在于，所述的根据转盘的旋 转频率的变化范围，W及正钻井周围一预设范围内的已钻井的最多数量，确定井间测距信 号的频带，采用如下公式： f ha 二 f hz X S fla=flz 其中，fha表示井间测距信号的最大频率;fla表示井间测距信号的最小频率;fhz表示转 盘的最大旋转频率，flz表示转盘的最小旋转频率;S表示正钻井周围所述预设范围内的已钻 井的最多数量。3. 根据权利要求2所述的提取井间测距信号的方法，其特征在于，所述的根据所述井间 测距信号的频带，W及预设的时域采样频率，确定上限数字频率和下限数字频率，采用如下 公式：其中，f S表示预设的时域采样频率;fhd表示上限数字频率;fid表示下限数字频率。4. 根据权利要求3所述的提取井间测距信号的方法，其特征在于，所述的在W所述下限 数字频率和所述上限数字频率为边界的频域采样范围中，对去除直流分量后的原始信号进 行离散傅里叶变换，将所述去除直流分量后的原始信号从时域转换至频域，得到频域信号， 采用如下公式：其中，Xc(fd〇表示频域信号;f。读示频域信号Xc(fd〇对应的数字频率;x(n)表示去除直 流分量后的原始信号;n表示时域采样点的编号;k表示频域采样点的编号;N表示时域采样 点总数;M表示频域采样点总数。5. 根据权利要求4所述的提取井间测距信号的方法，其特征在于，所述的在得到的全部 频域信号中捜索幅值最大的频域信号，包括： 分别W频域信号的数字频率和幅值为坐标标注得到的全部频域信号，得到幅度谱； 在所述幅度谱中捜索幅值最大的频域信号。6. 根据权利要求5所述的提取井间测距信号的方法，其特征在于，利用所述预设的时域 采样频率，计算所述幅值最大的频域信号对应的时域信号的频率和幅度，采用如下公式：其中，/。。表示时域信号的频率;J表示时域信号的幅值；IXc(fck) Ux表示所述幅值最大 的频域信号的幅值;fsk-max表示所述幅值最大的频域信号的数字频率。7. -种提取井间测距信号的装置，其特征在于，包括： 获取模块，用于获取电磁随钻测量系统采集的原始信号； 去直流模块，用于去除所述原始信号中的直流分量； 频带确定模块，用于根据转盘的旋转频率的变化范围，W及正钻井周围一预设范围内 的已钻井的最多数量，确定井间测距信号的频带； 边界频率确定模块，用于根据所述井间测距信号的频带，W及预设的时域采样频率，确 定上限数字频率和下限数字频率； 第一变换模块，用于在W所述下限数字频率和所述上限数字频率为边界的频域采样范 围中，对去除直流分量后的原始信号进行离散傅里叶变换，将所述去除直流分量后的原始 信号从时域转换至频域，得到频域信号； 捜索模块，用于在得到的全部频域信号中捜索幅值最大的频域信号； 第二变换模块，用于利用所述预设的时域采样频率，计算所述幅值最大的频域信号对 应的时域信号的频率和幅度； 结果确定模块，用于将该时域信号的频率和幅度确定为井间测距信号的频率和幅值。8. 根据权利要求7所述的提取井间测距信号的装置，其特征在于，所述的频带确定模块 采用如下公式确定井间测距信号的频带： f ha 二 f hz X S fla=flz 其中，fha表示井间测距信号的最大频率;fla表示井间测距信号的最小频率;fhz表示转 盘的最大旋转频率，flz表示转盘的最小旋转频率;S表示正钻井周围所述预设范围内的已钻 井的最多数量。9. 根据权利要求8所述的提取井间测距信号的装置，其特征在于，所述的边界频率确定 模块采用如下公式确定上限数字频率和下限数字频率：其中，f S表示预设的时域采样频率;fhd表示上限数字频率;fid表示下限数字频率。10. 根据权利要求9所述的提取井间测距信号的装置，其特征在于，所述的第一变换模 块采用如下公式将所述去除直流分量后的原始信号从时域转换至频域并得到频域信号：其中，扣化〇表示频域信号;f。读示频域信号Xc(fd〇对应的数字频率;x(n)表示去除直 流分量后的原始信号;n表示时域采样点的编号;k表示频域采样点的编号;N表示时域采样 点总数，M表示频域采样点总数。11. 根据权利要求10所述的提取井间测距信号的装置，其特征在于，所述的捜索模块包 括： 标注子模块，用于分别W频域信号的数字频率和幅值为坐标标注得到的全部频域信 号，得到幅度谱； 捜索子模块，用于在所述幅度谱中捜索幅值最大的频域信号。12. 根据权利要求11所述的提取井间测距信号的装置，其特征在于，所述的第二变换模 块采用如下公式计算所述幅值最大的频域信号对应的时域信号的频率和幅度：其中，完。表示时域信号的频率;名表示时域信号的幅值；|Xc(fd〇 kx表示所述幅值最大 的频域信号的幅值;fsk-max表示所述幅值最大的频域信号的数字频率。
【文档编号】E21B47/00GK105909233SQ201610284086
【公开日】2016年8月31日
【申请日】2016年4月29日
【发明人】梁华庆, 窦新宇, 耿敏, 高德利
【申请人】中国石油大学(北京)