CN105822289B - 用于油井动液面检测的频率估算方法 - Google Patents

用于油井动液面检测的频率估算方法 Download PDF

Info

Publication number
CN105822289B
CN105822289B CN201610177653.6A CN201610177653A CN105822289B CN 105822289 B CN105822289 B CN 105822289B CN 201610177653 A CN201610177653 A CN 201610177653A CN 105822289 B CN105822289 B CN 105822289B
Authority
CN
China
Prior art keywords
frequency
window
signal
maximum
value
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
Application number
CN201610177653.6A
Other languages
English (en)
Other versions
CN105822289A (zh
Inventor
罗久飞
周伟
周盼
李晓亮
李太福
易军
张元涛
胡刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201610177653.6A priority Critical patent/CN105822289B/zh
Publication of CN105822289A publication Critical patent/CN105822289A/zh
Application granted granted Critical
Publication of CN105822289B publication Critical patent/CN105822289B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/04—Measuring depth or liquid level
    • E21B47/047—Liquid level
    • 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/14—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 using acoustic waves
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Business, Economics & Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Marketing (AREA)
  • Animal Husbandry (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Agronomy & Crop Science (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种用于油井动液面检测的频率估算方法，包括以下步骤：采集油井动液面的声场信号，得到采样信号x(n)；对采样信号x(n)进行加窗处理，得到加窗函数xw(n)；将加窗函数xw(n)进行离散傅里叶变换，得到频谱Xw(k)；从频谱Xw(k)中寻找幅值最大的频点、最大频谱幅值Xw(l)和第二大频谱幅值Xw(l±1)，其中幅值最大的频点记为第一次迭代时的频率初始值l1；通过插值计算真实频率值λ0；按照公式计算信号的频率fo。有益效果：本算法运用范围广，适用于所有选择的窗函数；计算误差小，受噪声的影响小，具有较好的一致性；计算过程简洁，无需对数据进行预先计算，也不需要存储窗函数的相关参数信息。

Description

用于油井动液面检测的频率估算方法
技术领域
本发明涉及油井动液面测量中信号频率的检测技术领域，具体的说是一种用于油井动液面检测的频率估算方法。
背景技术
目前已提出的基于管柱声场特性的油井动液面检测方法，其原理是在井口发出声波激发井下套管内空气柱的共振，利用检测到的空气柱声波共振频率来估算空气柱长度，进而测量出空气柱下面的动液面深度。其数学模型为
fn为油套环空内空气柱的第n阶共振频率，c为井下声波的传播速度，d为井下管道的直径。由此数学模型可知，在声速一定的情况下，只要准确测量井下空气柱的共振频率fn即可得估算出动液面深度。因此，准确高效的估算出井下空气柱的共振频率是测量动液面深度的关键。
但是由于事先并不知道声波信号的频率，很难做到对信号的整周期采样，因此对采集到信号的傅里叶变换一般都会出现频谱泄漏和栅栏效应，导致对信号频率的估算存在很大的偏差，人们常采用加窗、插值迭代法来克服缺陷，但是对于不同的窗函数的插值公式不一样，没有适用于所用窗函数的通用算法，不能满足人们的需求，而且目前这些算法对谱线错误定位敏感，在大噪声条件下误差较大。
发明内容
针对上述问题，本发明提供了一种用于油井动液面检测的频率估算方法，来计算油井动液面的声场信号频率，来解决声场信号频率计算过程复杂、计算误差大的缺陷。
为达到上述目的，本发明采用的具体技术方案如下：
一种用于油井动液面检测的频率估算方法，包括以下步骤：
S1：采集油井动液面的声场信号，得到M点长度的采样信号x(n)：
其中Ao表示幅值，fo表示信号的频率，fs表示采样频率，表示相位，n表示采样点序数，M为对信号的总的采样点数；
S2：对步骤S1得到的采样信号x(n)进行加窗处理；
构造K点长度的窗函数w(n)，对采样信号x(n)进行加窗，得到函数：
xw(n)＝x(n)w(n)；
S3：将步骤S2得到的加窗函数xw(n)进行离散傅里叶变换，得到频谱Xw(k)：
S4：从步骤S3所得的频谱Xw(k)中寻找幅值最大的频点以及当前最大频谱幅值Xw(l)和第二大频谱幅值Xw(l±1)；其中幅值最大的频点记为第一次迭代时的频率初始值l1；
S5：计算真实频率值λ0；
根据公式lm+1＝lm+δm计算第m+1次迭代时频率估算值，且第m+1次迭代时频率估算值和第m个频率估算值偏差其中：
q为矩形窗与选择的窗函数的扇形损失比，|Xw(lm)|为第m次迭代时频率估算值对应的频谱幅值，|X'w(lm-1)|为|Xw(lm-1)|、|Xw(lm-0.5)|以及|Xw(lm+0.5)|之间的最大值，Xw(lm-1)为第m-1次迭代时频率估算值对应的频谱幅值，Xw(lm-0.5)为第m次迭代时频率估算值减0.5个步长所在频点的幅值，Xw(lm+0.5)为第m次迭代时频率估算值加0.5 个步长所在频点的幅值，当|lm+1-lm|＜τ时，τ为预设阈值，迭代结束，并设定真实频率值λ0＝lm+1；
S6：将步骤S5得到的真实频率值值λ0带入公式计算得到信号的频率fo。
进一步描述，当所加窗函数为矩形窗时，对其进行离散傅里叶变换得到：
再进一步描述，当所加窗函数为三角窗时，对其进行离散傅里叶变换得到：
当所加窗函数为汉宁窗时，对其进行离散傅里叶变换得到：
再进一步描述，步骤S5中矩形窗与选择的窗函数的扇形损失比q为： 其中SLr矩形窗的扇形损失，SLc为选择的窗函数的扇形损失，且窗函数的扇形损失为：其中W(0.5)为信号最差非整周期采样的情形时候的最大谱线幅值，W(0)为信号整周期采样时最大谱线幅值。
本发明的有益效果：本算法运用范围广，适用于所有的窗函数；计算误差小，可以克服噪声对谱线造成的影响；计算过程简洁，无需对数据进行预先计算，也不需要存储选择窗函数的相关参数信息。
附图说明
图1是本发明的声场信号频率计算流程图；
图2是常用窗函数计算信号频率系统误差图；
图3是随机噪声对不同信号频率算法的影响效果图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
从图1可以看出，一种用于油井动液面检测的频率估算方法，包括以下步骤：
S1：采集油井动液面的声场信号，得到M点长度的采样信号x(n)：
其中Ao表示幅值，fo表示信号的频率，fs表示采样频率，表示相位，n表示采样点序数，M为对信号的总的采样点数；
S2：对步骤S1得到的采样信号x(n)进行加窗处理；
构造K点长度的窗函数w(n)，对采样信号x(n)进行加窗，得到函数：xw(n)＝x(n)w(n)；
在本实施例中，所加窗函数为矩形窗，则其表达式为：w(n)＝1(n＝0,1,2…K)，则；xw(n)＝x(n)(n＝1,2,…K)
S3：将步骤S2得到的加窗函数xw(n)进行离散傅里叶变换，得到频谱Xw(k)：
S4：从步骤S3所得的频谱Xw(k)中寻找幅值最大的频点以及当前最大频谱幅值Xw(l)和第二大频谱幅值Xw(l±1)；其中幅值最大的频点记为第一次迭代时的频率初始值l1；
S5：计算真实频率值λ0：
在实际应用中，由于不知道信号的实际频率，因此对信号的截断很难是整周期，那么进行DFT变换就会出现频谱泄露，即归一化的真实频率值λ0几乎都是出现在DFT变换后的两个谱线间，因此真实频率值λ0假设表示为：
λ0＝lm+δm (2)
其中lm和δm(-0.5≤δ≤0.5)分别为整数和小数；δm为选择的窗函数第m+1次迭代时频率估算值和第m个频率估算值偏差，δm具体计算方法为：
将公式(2)带入公式(1)可得到：此处假设5＜λ0＜N/2-5，则：|W(2lm+δm)|＜＜|W(-δm(-δm)|，最大频率幅值的表达式为：
第二大的频谱幅值表达式为：
令最大与第二大频谱幅值之比为：
对于任意窗函数的归一化频谱区间[-0.5,0.5]内有：
[Wc(k)]q≈Wr(k)
其中Wc(k)、Wr(k)表示任意窗函数和矩形窗的频谱，即[Wc(k)]q与Wr(k)在归一化频谱区间[-0.5,0.5]内形状几乎一样,将任意的窗函数逼近矩形窗，则矩形窗对应的最大与第二大频谱幅值之比为：
其中q为矩形窗与选择的窗函数的扇形损失比，且：其中SLr矩形窗的扇形损失，SLc为选择的窗函数的扇形损失，且窗函数的扇形损失为： 其中W(0.5)为信号最差非整周期采样的情形时候的最大谱线幅值，W(0)为信号整周期采样时最大谱线幅值；
对于矩形窗，即wr(n)＝1(n＝0,1,2…N)，对其进行离散傅里叶变换(DFT)可得到：当|k|/N＜＜1，则：
结合公式(3)(4)(5)(6)可得到矩形窗函数存在以下关系：
则频率偏差
在本实施例中，选择所加的窗函数为矩形窗函数，则矩形窗与选择的窗函数的扇形损失比则频率偏差
根据步骤S4寻找的第一次迭代时的频率初始值l1以及当前最大频谱幅值Xw(l1)和第二大频谱幅值Xw(l1±1)，可计算出第2次迭代时频率估算值和第1次迭代时频率估算值偏差其中则第2次迭代时频率估算值为l2＝l1+δ1，若满足δ1＝|l2-l1|＜τ时，τ为预设阈值，此时第2次迭代时频率估算值l2为真实频率值λ0；
若δ1＝|l2-l1|＞τ，则迭代继续，按照公式计算第m+1次迭代时频率估算值和第m次迭代时频率估算值偏差；其中，为了计算第m次迭代时最大与第二大频谱幅值之比αrm，就要要寻找最大频谱幅值和第二大频谱幅值；
首先比较第m-1次迭代时频率估算值对应的频谱幅值Xw(lm-1)、第m次迭代时频率估算值减0.5个步长所在频点的幅值Xw(lm-0.5)以及第m次迭代时频率估算值加0.5个步长所在频点的幅值Xw(lm+0.5)大小，三者中的最大值设为|X'w(lm-1)|；然后比较|X'w(lm-1)|与当第m次迭代时频率估算值对应的频谱幅值|Xw(lm)|的大小，若|X'w(lm-1)|＜|Xw(lm)|，则为|Xw(lm)|为最大频谱幅值，|X'w(lm-1)|为第二大频谱幅值；若|X'w(lm-1)|≥|Xw(lm)|，|X'w(lm-1)|最大频谱幅值，|X'w(lm-1)|为第二大频谱幅值；
即当|X'w(lm-1)|＜|Xw(lm)|时，
当|X'w(lm-1)|≥|Xw(lm)|时，按照公式lm+1＝lm+δm计算下一次迭代时频率估算值，当满足|lm+1-lm|＜τ，此时得到的第m+1个频率估算值为真实频率值λ0；
S6：将步骤S5得到的真实频率值λ0带入公式计算得到信号的频率fo。
计算系统误差：
设置采样信号x(n)的幅值A0设为1，采样频率fs为1024Hz，数据总长度M为1024点，故频率分辨率为1Hz，对每个窗函数的离散傅里叶变换谱线k，在区间[k–0.5，k+0.5)内以步距0.05进行频率扫描。对于每一个步距频率，相位在区间(-π，π]范围内变化，步距为π/72。即对于每个窗函数的离散傅里叶变换谱线产生2880个测试信号，其中与理论信号的最大绝对误差作为该谱线k的系统误差。
从图2可以看出，不同窗函数在该方法下计算出的信号的频率fo的系统误差，其中窗函数包括矩形窗(Rectangle Window)、海明窗(Hamming Window)、三角窗(TriangleWindow)、汉宁窗(Hanning Window)、布莱克曼窗(Blackman Window),从图2可以看出，λ0满足λ0＞5和λ0＜N/2-5系统误差是非常小的， 可被忽略，对于所有窗函数，最大频率误差约为10-3，这足以满足大多数工程应用，同时表明该算法可适用于不同的窗函数。
检测随机噪声对不同信号频率算法的影响：
设置采样信号x(n)的幅值A0设为1，采样频率fs为1024Hz，数据总长度M为1024点，故频率分辨率为1Hz，为了使得谱线定位错误的几率较高，将信噪比设为–5dB，为了充分验证谱线定位错误对本算法的影响仿真信号的频率在[255.5，256.5]区间内逐渐增加，步距为0.025，仿真信号的相位在[-π,π)范围内随机选取。每个频率均独立产生50000个带有加性高斯白噪声干扰的实例，观察不同算法的方差；从图3可以看出，无论是加矩形窗还是汉宁窗，本算法(ACIIpDFT)都具有最小的方差，且方差不随频率偏差δ改变而改变，具有较好的一致性。因此，本算法不但具有较好的噪声性且不会受到谱线定位错误的影响。

Claims (2)

1.一种用于油井动液面检测的频率估算方法，其特征在于包括以下步骤：
S1：采集油井动液面的声场信号，得到M点长度的采样信号x(n)：
其中Ao表示幅值，fo表示信号的频率，fs表示采样频率，表示相位，n表示采样点序数，M为对信号的总的采样点数；
S2：对步骤S1得到的采样信号x(n)进行加窗处理；
构造K点长度的窗函数w(n)，对采样信号x(n)进行加窗，得到函数：xw(n)＝x(n)w(n)；
S3：将步骤S2得到的加窗函数xw(n)进行离散傅里叶变换，得到频谱Xw(k)；
S4：从步骤S3所得的频谱Xw(k)中寻找幅值最大的频点以及当前最大频谱幅值Xw(l)和第二大频谱幅值Xw(l±1)；其中幅值最大的频点记为第一次迭代时的频率初始值l1；
S5：计算真实频率值λ0；
根据公式lm+1＝lm+δm计算第m+1次迭代时频率估算值，且第m+1次迭代时频率估算值和第m个频率估算值偏差其中：q为矩形窗与选择的窗函数的扇形损失比，|Xw(lm)|为第m次迭代时频率估算值对应的频谱幅值，|X'w(lm-1)|为|Xw(lm-1)|、|Xw(lm-0.5)|以及|Xw(lm+0.5)|之间的最大值，Xw(lm-1)为第m-1次迭代时频率估算值对应的频谱幅值，Xw(lm-0.5)为第m次迭代时频率估算值减0.5个步长所在频点的幅值，Xw(lm+0.5)为第m次迭代时频率估算值加0.5个步长所在频点的幅值，当|lm+1-lm|＜τ时，τ为预设阈值，迭代结束，并设定真实频率值λ0＝lm+1；
由于不知道信号的实际频率，因此对信号的截断很难是整周期，那么进行DFT变换就会出现频谱泄露，即归一化的真实频率值λ0几乎都是出现在DFT变换后的两个谱线间，因此真实频率值λ0假设表示为：
λ0＝lm+δm (2)
其中lm和δm(-0.5≤δ≤0.5)分别为整数和小数；δm为选择的窗函数第m+1次迭代时频率估算值和第m个频率估算值偏差，δm具体计算方法为：
将公式(2)带入公式(1)可得到：此处假设5＜λ0＜N/2-5，则：|W(2lm+δm)|＜＜|W(-δm)|，最大频率幅值的表达式为：
第二大的频谱幅值表达式为：
令最大与第二大频谱幅值之比为：
对于任意窗函数的归一化频谱区间[-0.5,0.5]内有：
[Wc(k)]q≈Wr(k)
其中Wc(k)、Wr(k)表示任意窗函数和矩形窗的频谱，即[Wc(k)]q与Wr(k)在归一化频谱区间[-0.5,0.5]内形状几乎一样,将任意的窗函数逼近矩形窗，则矩形窗对应的最大与第二大频谱幅值之比为：
其中q为矩形窗与选择的窗函数的扇形损失比，且：其中SLr矩形窗的扇形损失，SLc为选择的窗函数的扇形损失，且窗函数的扇形损失为：其中W(0.5)为信号最差非整周期采样的情形时候的最大谱线幅值，W(0)为信号整周期采样时最大谱线幅值；
S6：将步骤S5得到的真实频率值λ0带入公式计算得到信号的频率fo；
当所加窗函数为矩形窗时，对其进行离散傅里叶变换为得到：
当所加窗函数为三角窗时，对其进行离散傅里叶变换得到：
当所加窗函数为汉宁窗时，对其进行离散傅里叶变换得到：
2.根据权利要求1所述用于油井动液面检测的频率估算方法，其特征在于：步骤S5中矩形窗与选择的窗函数的扇形损失比q为：其中SLr矩形窗的扇形损失，SLc为选择的窗函数的扇形损失，且窗函数的扇形损失为：其中W(0.5)为信号最差非整周期采样的情形时候的最大谱线幅值，W(0)为信号整周期采样时最大谱线幅值。
CN201610177653.6A 2016-03-25 2016-03-25 用于油井动液面检测的频率估算方法 Active CN105822289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610177653.6A CN105822289B (zh) 2016-03-25 2016-03-25 用于油井动液面检测的频率估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610177653.6A CN105822289B (zh) 2016-03-25 2016-03-25 用于油井动液面检测的频率估算方法

Publications (2)

Publication Number Publication Date
CN105822289A CN105822289A (zh) 2016-08-03
CN105822289B true CN105822289B (zh) 2019-03-19

Family

ID=56525279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610177653.6A Active CN105822289B (zh) 2016-03-25 2016-03-25 用于油井动液面检测的频率估算方法

Country Status (1)

Country Link
CN (1) CN105822289B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107271002A (zh) * 2017-06-19 2017-10-20 重庆邮电大学 一种快速高精度的频谱校正插值算法
CN112664177B (zh) * 2020-12-29 2022-06-21 重庆邮电大学 一种基于归一化卷积与自适应滤波的油井动液面测量方法
CN114923690A (zh) * 2022-04-18 2022-08-19 南京航空航天大学 一种高精度轴承故障特征频率估计诊断方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1895433A1 (fr) * 2006-08-29 2008-03-05 France Télécom Procédé d'estimation de phase pour la modélisation sinusoidale d'un signal numérique
CN103321633A (zh) * 2013-06-25 2013-09-25 重庆科技学院 油井动液面深度检测方法及装置
CN103821499A (zh) * 2014-02-20 2014-05-28 重庆科技学院 用于油井动液面深度检测的声音信号处理方法
CN104597321A (zh) * 2015-01-28 2015-05-06 常洪山 基于四条离散傅里叶复数谱线的信号频率测量方法及装置
CN105137185A (zh) * 2015-07-23 2015-12-09 河海大学 一种基于离散傅里叶变换的频域插值电力谐波分析方法
CN105353215A (zh) * 2015-12-23 2016-02-24 合肥工业大学 基于Nuttall窗四谱线插值FFT的谐波检测方法
CN105372492A (zh) * 2014-08-31 2016-03-02 盛吉高科（北京）科技有限公司 基于三条dft复数谱线的信号频率测量方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1895433A1 (fr) * 2006-08-29 2008-03-05 France Télécom Procédé d'estimation de phase pour la modélisation sinusoidale d'un signal numérique
CN103321633A (zh) * 2013-06-25 2013-09-25 重庆科技学院 油井动液面深度检测方法及装置
CN103821499A (zh) * 2014-02-20 2014-05-28 重庆科技学院 用于油井动液面深度检测的声音信号处理方法
CN105372492A (zh) * 2014-08-31 2016-03-02 盛吉高科（北京）科技有限公司 基于三条dft复数谱线的信号频率测量方法
CN104597321A (zh) * 2015-01-28 2015-05-06 常洪山 基于四条离散傅里叶复数谱线的信号频率测量方法及装置
CN105137185A (zh) * 2015-07-23 2015-12-09 河海大学 一种基于离散傅里叶变换的频域插值电力谐波分析方法
CN105353215A (zh) * 2015-12-23 2016-02-24 合肥工业大学 基于Nuttall窗四谱线插值FFT的谐波检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Frequency estimation of the weighted real tones or resolved multiple tones by iterative interpolation DFT algorithm;Jiufei Luo ,et al.;《Digital Signal Processing》;20150310;第41卷;第118-129页

Also Published As

Publication number Publication date
CN105822289A (zh) 2016-08-03

Similar Documents

Publication Publication Date Title
US10472957B2 (en) Entropy based multiphase flow detection
CN105822289B (zh) 用于油井动液面检测的频率估算方法
US7209752B2 (en) Error estimate concerning a target device's location operable to move in a wireless environment
US10451475B2 (en) Gauge length optimization in distributed vibration sensing
CN108444653A (zh) 一种面阵型超声波气体泄漏点定位检测装置及方法
CN101852638B (zh) 基于固定频段声波共振频率的液位测量方法
CN113494990A (zh) 一种风洞扰动对超声速层流边界层厚度影响的分析方法
JP6132990B2 (ja) 状態推定装置
Xue et al. Stall recognition for centrifugal compressors during speed transients
EP2957716A2 (en) Determining a quantitative bond using signal attenuation
CN102928713B (zh) 一种磁场天线的本底噪声测量方法
EP3775670A1 (en) Method and system for non-intrusively determining cross-sectional variation for a fluidic channel
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
AU2014208222A1 (en) Methods, systems and devices for generating slowness-frequency projection logs
Yu et al. Adaptive PIV algorithm based on seeding density and velocity information
US8306762B2 (en) Systems and methods for analysis of downhole data
KR101354522B1 (ko) 변분법과 초기추정법을 이용한 고분해능 레이더 바람장 산출방법
CN107271002A (zh) 一种快速高精度的频谱校正插值算法
NL2022689B1 (en) Method and system for non-intrusively determining cross-sectional variation for a fluidic channel
CN101825722B (zh) 一种鲁棒的地震信号瞬时频率的估计方法
CN108106717B (zh) 一种基于声音信号识别机组状态的方法
CN109855536B (zh) 一种基于应变测量的油气管道堵塞检测方法
CN114859404A (zh) 超采样地震波形匹配方法及装置
CN105626039B (zh) 一种生产测井相关流量数据预处理的方法
US10990846B2 (en) Pattern matching to detect defects in rod lift downhole cards

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20160803

Assignee: WUZHOU JINZHENGYUAN ELECTRONIC TECH. Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980053985

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20231227

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20160803

Assignee: Guangzhou Yuming Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000306

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240110

Application publication date: 20160803

Assignee: Guangzhou chuangyixin Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000305

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240110

Application publication date: 20160803

Assignee: Guangzhou Shangqi Network Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000303

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240110

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20160803

Assignee: GUANGZHOU XUNSU PHOTOELECTRIC TECHNOLOGY CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000630

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240119

Application publication date: 20160803

Assignee: GUANGZHOU CHENGKE ELECTRONIC TECHNOLOGY CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000629

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240119

Application publication date: 20160803

Assignee: Zhongke Liben (Chongqing) Technology Development Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000644

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240119

Application publication date: 20160803

Assignee: GUANGZHOU DIYUE NETWORK TECHNOLOGY CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000631

Denomination of invention: Frequency estimation method for dynamic liquid level detection in oil wells

Granted publication date: 20190319

License type: Common License

Record date: 20240119