CN105738257A - 测量方法及系统 - Google Patents

Abstract

本发明揭示一种测量方法及系统。该测量方法包括：发射若干超声信号至被测流体；接收从被测流体中的颗粒反射或散射回来的回波信号；根据所述回波信号计算与颗粒属性有关的R参数的计算值；确定所述与颗粒属性有关的R参数的理论曲线；及根据所述R参数的计算值和所述R参数的理论曲线确定流体中颗粒的平均直径；根据所述颗粒的平均直径确定流体中颗粒的体积浓度。本发明也提供了一种测量系统，用来测量颗粒的平均直径和体积浓度。

Description

测量方法及系统

技术领域

本发明有关一种测量方法和系统，尤其涉及一种用来测量流体中悬浮颗粒的平均直径和体积浓度的测量方法及系统。

背景技术

在利用多普勒原理测量流体流量的方法中，悬浮于流体中的颗粒反射或散射超声探头发射出的超声信号。反射或散射回来的多普勒回波信号的能量与颗粒的直径和体积浓度等有关。当颗粒的直径较小或体积浓度较低时，回波信号较弱；颗粒的直径增大或体积浓度升高时，回波信号变大。然而，一方面，当颗粒直径太大时，由于自重过大，悬浮颗粒相对于流体本身将会存在较大的相对速度，这时如果仍然将悬浮颗粒的运动速度作为流体本身的运动速度将会降低流量测量的精度；另一方面，当颗粒体积浓度太大时，超声信号在流体中的衰减较大，探头接收到的回波信号减弱，较低的信噪比也会影响流量的测量。如此，在利用多普勒原理测量流体流量的过程中，测量颗粒的直径和体积浓度显得尤其重要。

发明内容

本发明的一个方面在于提供一种测量方法。该测量方法包括：发射若干超声信号至被测流体；接收从被测流体中的颗粒反射或散射回来的回波信号；根据所述回波信号计算与颗粒属性有关的R参数的计算值；确定所述与颗粒属性有关的R参数的理论曲线；及根据所述R参数的计算值和所述R参数的理论曲线确定流体中颗粒的平均直径；根据所述颗粒的平均直径确定流体中颗粒的体积浓度。

本发明的另一个方面在于提供一种测量系统。该测量系统包括：超声探头，用来发射若干超声信号至被测流体且接收从被测流体中的颗粒反射或散射回来的回波信号；及处理器，连接于所述超声探头，且用来根据所述回波信号计算与颗粒属性有关的R参数的计算值、确定所述与颗粒属性有关的R参数的理论曲线、根据所述R参数的计算值和所述R参数的理论曲线确定流体中颗粒的平均直径，且根据所述颗粒的平均直径确定流体中颗粒的体积浓度。

附图说明

通过结合附图对于本发明的实施方式进行描述，可以更好地理解本发明，在附图中：

图1所示为本发明测量系统的一个实施例的示意图；

图2所示为本发明测量方法的一个实施例的流程图；

图3所示为本发明测量方法的R参数的理论曲线的一个实施例的曲线图；

图4所示为安装于管道的多个超声探头的示意图。

具体实施方式

除非另作定义，此处使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明专利申请说明书以及权利要求书中使用的“第一”“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。同样，“一个”或者“一”等类似词语也不表示数量限制，而是表示存在至少一个。除非另行指出，“前部”“后部”“下部”和/或“上部”等类似词语只是为了便于说明，而并非限于一个位置或者一种空间定向。另外，“连接”或者“相连”等类似的词语并非用来区分两个元件之间的直接或间接连接。当然，除非另行说明，此元件间可以直接或间接连接。

图1所示为一个实施例的测量系统10的示意图。测量系统10可用来测量流体20的颗粒大小和颗粒体积浓度。流体20可以沿图中箭头方向在管道24等内流动。流体20中悬浮有颗粒22，例如碎屑、气泡、泥等杂质或污染物。颗粒22的形状、大小和分布一般不确定，往往根据实际运用环境变化。在一个典型的运用中，测量系统10可用于石油或天然气开采钻井中来测量钻井过程中回流的液体的颗粒大小和浓度，该回流的液体中往往悬浮有钻井产生的石头碎屑、泥等杂质。但测量系统10并不限于该运用领域，还可应用于其他领域。

测量系统10包括超声探头11、超声收发装置13、数字转换器17和处理器19。超声探头11可安装于管道24的壁上。超声收发装置13产生并发出电信号给超声探头11。超声探头11将该电信号转换成超声信号(超声波)32，发射至被测流体20。超声探头11发射的超声信号可覆盖流体20的部分区域。一个超声探头11发射的超声信号可透射一个近似于锥形的区域21。超声探头11一般包括压电元件，能够将电信号转换成物理脉冲(本实施例中为超声信号32)，超声信号32的频率由电信号的频率控制。本实施例中，超声探头11发射若干具有不同频率(f₁,f₂,…f_n)的超声信号32至被测流体20。在另一个实施例中，超声探头11发射单一频率的超声信号32至被测流体20。图1中仅示出了一个超声探头11。在一些实施例中，多个超声探头11可用来发射超声信号。

超声信号32可被颗粒22反射或散射回来。至少部分从被测流体20中的颗粒22反射或散射回来的回波信号34被超声探头11接收。超声探头11将回波信号34转换成电信号。超声收发装置13接收超声探头11输出的电信号且将该电信号的幅值放大。数字转换器17将超声收发装置13放大的电信号转换成数字信号给处理器19。

处理器19包括编码在固定存储器36或可拆卸的存储器38内的程序指令，通过处理来自数字转换器17的数字信号产生颗粒的平均直径和/或颗粒的体积浓度。在一些实施例中，处理器19还用来产生流体20的流速、流量等。在一个实施例中，处理器19为具有中央处理器(CPU)、只读存储器(ROM)、随机存储器(RAM)等的微计算机。处理器19连接于显示器40来输出颗粒大小、颗粒的体积浓度等流体20的信息。

图2所示为一个实施例的测量方法50的流程图。测量方法50可用来测量流体中颗粒的平均直径和颗粒的体积浓度。测量方法50可以通过图1所示的测量系统10实现，但不限于此。步骤51中，发射超声信号至被测流体。超声信号可通过超声探头发射至被测流体中的至少一个区域。一个超声探头发出的超声信号可透射流体的一个近似锥形的区域。若干个超声探头可以发射超声信号至流体的若干近似锥形的区域。在一个实施例中，发射单一频率的超声信号。在另一个实施例中，发射若干具有不同频率(f₁,f₂,…f_n)的超声信号。超声信号的频率可以根据被测颗粒的大小进行调整，用低频率的超声信号来测量较大的颗粒，用高频率的超声信号来测量较小的颗粒。

步骤52中，接收从被测流体中的颗粒反射或散射回来的回波信号。步骤51中发射出的超声信号传递至流体中的颗粒，被颗粒反射或散射，至少部分反射或散射回来的回波信号到达超声探头，被超声探头接收。回波信号被转换成电信号，并进行放大且转换成数字信号进行处理。

步骤53中，根据回波信号计算与颗粒属性有关的R参数的计算值。颗粒的属性包括超声探头所覆盖的区域内的颗粒的大小、分布、密度、反射系数，不包括颗粒的浓度。R参数与颗粒的浓度无关。一个超声探头发出的超声信号透射的区域内可以计算出至少一个R参数。在一个实施例中，一个超声探头发出的超声信号透射的区域对应有一个R参数。在另一个实施例中，一个超声信号透射的区域可以进一步划分为多个小区域，每一小区域计算出一个R参数。该计算R参数的计算值的步骤53包括确定回波信号的电压幅度均方值关于颗粒属性的函数。在本实施例中，该函数还与流体的属性、探头的属性和测量系统的属性等有关。然后根据该函数计算R参数的计算值。本实施例中，假设颗粒均匀地悬浮在流体中，且来自各个颗粒的回波信号不相干。采用回波信号的均方根函数，该函数表达式如下面的函数式(1)所示：

$V_{\mathrm{rma}}\left(r\right)=\frac{K_s{K}_t}{\mathrm{r\ψ}}\sqrt{M}{e}^{-2r({\mathrm{\α}}_m+{\mathrm{\α}}_s)}---\left(1\right)$

其中，V_rms是回波信号的电压幅度均方值。r是探头到被测颗粒的距离，可通过采样时间换算获得。K_s、M和α_s为与颗粒大小和/或体积浓度有关的参数。K_s是反映颗粒大小、颗粒分布、颗粒密度和颗粒散射特性的参数。M为颗粒的体积浓度。α_m和α_s均为衰减系数，其中，α_m反映了超声信号在流体中的衰减，可测量获得。α_s反映了流体中的悬浮颗粒对超声信号造成的衰减。在钻井的运用中，α_m主要反映了钻井泥浆对超声信号造成的衰减，α_s反映了钻井泥浆中的悬浮颗粒对超声信号造成的衰减。K_t是与测量系统和超声探头的固定属性等有关的系统常量。Ψ是超声探头的近场校正因子，可根据经验公式计算获得。

进一步地，K_s可通过下面的表达式(2)计算获得：

$K_s=<f>/\sqrt{<a_s>{\mathrm{\&rho;}}_s}---\left(2\right)$

其中，<f>为颗粒的反向散射形态函数平均值(averagebackscatterformfunction)，其反映单位体积内的颗粒反向散射超声的能量的能力，其与颗粒的大小分布相关。<a_s<为颗粒平均直径。ρ_s为颗粒密度。

α_s可通过下面的表达式(3)计算获得：

α_s＝3M<χ>/4<a_s>ρ_s(3)

其中，<χ>为归一化总散射截面平均值(averagenormalizedtotalscatteringcross-section)。

对于一个频率(例如频率f₁)的超声信号，定义参数及参数B＝α_m+3/4A²<χ>。根据自定义的参数A和B以及K_s和α_s的表达式(2)和(3)，将函数式1转换成下面的函数式(4)：

$V_{\mathrm{rma}}\left(r\right)=\frac{A<f>K_t}{\mathrm{r\&psi;}}{e}^{-2\mathrm{rB}}---\left(4\right)$

从函数式(4)可以计算出指数-2rB，其中颗粒到探头的距离r可以通过采样时间换算得到，如此可以计算出参数B。

根据参数B的表达式，进一步定义参数且根据已经计算出的参数B的值可以计算出参数C₁的值。定义参数C₂＝A<f>，且根据参数B的值和函数式(4)计算出C₂的值。定义R参数R＝C₂/C₁，根据C₁和C₂的值可以计算出R参数的计算值。且根据C₁和C₂的表达式可以得到R参数的另一计算公式由该R参数的计算公式可以看出，R参数与<f>和<χ>有关。而<f>和<χ>则与颗粒大小分布、颗粒和流体的密度、颗粒和流体的弹性系数有关。在本实施例中，假设颗粒的分布为高斯分布，且进一步假设高斯分布中的颗粒直径的方差已知，假设颗粒和流体的密度及弹性系数已知，如此R参数随着颗粒的平均直径的变化而变化，其并不反映颗粒浓度。在其他实施例中，可以假设颗粒分布为其他形式的分布。颗粒直径的方差可以根据经验值设定，颗粒和流体的密度及弹性系数可以通过实验等方式获得。如此将只与颗粒平均直径有关的参数独立出来进行后续计算。在一个实施例中，对于不同的超声频率f₁-f_n，可根据上述方法分别计算出相应频率下的R参数的计算值R₁-R_n。

步骤54中，确定与颗粒属性有关的R参数的理论曲线。R参数相对于颗粒的平均直径变化的理论曲线可以根据计算公式计算获得。在一个实施例中，可以计算获得不同频率f₁-f_n下的R参数的理论曲线。在一个实施例中，R参数的理论值R_{theoretical,i}可以预先计算后存储在图1中的处理器19的存储器36、38中。

请参考图3，图3图示了一种实施例的R参数的理论曲线。图3实施例中的颗粒的密度为4300kg/m³，颗粒的平均直径方差为0.1mm，且超声在流体中的速度为1240m/s。图3显示了在此条件下R参数的理论曲线，R参数的理论曲线体现了在不同频率下R参数和颗粒的平均直径之间的关系。设置不同的颗粒的密度、平均直径方差和超声在流体中的速度时，R参数的理论曲线会改变。根据实际应用，设置不同的条件，获得与实际条件相符的理论曲线。

继续参考图2，步骤55中，根据R参数的计算值和R参数的理论曲线确定颗粒的平均直径。

颗粒的平均直径<a_s>的表达式如下面的表达式(5)所示：

$<a_s>=\underset{<a_s>}{\arg \min }{\mathrm{\&Sigma;}}_{i=1}^n{(R_i-R_{\mathrm{theoretical},i})}^2---\left(5\right)$

R_i为计算得到的第i个超声频率f_i下R参数的计算值，R_{theoretical,i}为第i个超声频率f_i下R参数的理论值。拟合R参数的计算值R_i和R参数的理论曲线获得颗粒的平均直径。在一个实施例中，寻求颗粒的平均直径<a_s>的值使得不同频率下的若干R参数的计算值R_i和若干R参数的理论值R_{theoretical,i}之间的差值的平方和最小。在一个实施例中，超声频率为单一频率，此时n＝1。如此，可以获得某一区域内的颗粒的平均直径。

步骤56中，根据颗粒的平均直径确定颗粒的体积浓度M。根据颗粒的平均直径<a_s>和参数A计算颗粒的体积浓度M。如此可以获得相应区域内的颗粒的体积浓度M。该测量方法50的计算步骤简单，不需要进行迭代，可以将颗粒大小相关的参数和颗粒体积浓度相关的参数分开进行计算。且可测量尺寸较大的颗粒，适用的颗粒大小范围较广。如图3所示，在该实施例中，本方法可以测量平均直径达到0.02m的颗粒。此外，在一些实施例中，本方法也可以仅使用单个频率的超声探头完成颗粒平均直径和体积浓度的测量，利于实际应用和部署。

在一个实施例中，测量方法50进一步包括步骤57，确定管道区域内的颗粒大小分布。结合参考图4，若干超声探头11间隔排布于管道24，超声探头11发射超声信号至被测流体20的若干区域21。在一个实施例中，多个超声探头11沿管道24的纵长方向间隔排布。在另一个实施例中，多个超声探头11沿管道24的圆周方向间隔排布。在再一个实施例中，多个超声探头11沿管道24的纵长方向和圆周方向间隔排列。超声探头11的排布并不限于此，可以根据实际应用设置。每一个超声探头11可发射超声信号覆盖近似圆锥体的区域21。多个超声探头11透射的区域不重叠或部分重叠。在一个实施例中，同一个超声探头11可以发射不同的频率。在另一个实施例中，不同的超声探头11发射不同的频率。根据步骤51-56，分别测得每一区域内的颗粒的平均直径和体积浓度，然后根据该些区域内的颗粒的平均直径和体积浓度确定管道区域内的颗粒大小分布。可以确定多个超声探头11覆盖到的管道内的区域的颗粒大小分布。在一些实施例中，可以测得整个管道区域内的颗粒大小分布。

在一个实施例中，颗粒大小、体积浓度和/或颗粒大小分布的数据信息可被显示。在另一个实施例中，颗粒大小、体积浓度和/或颗粒大小分布可用于其他计算处理和/或控制步骤，例如，颗粒大小、体积浓度和/或颗粒大小分布的数据信息可以用来辅助流体流速的计算，或者根据颗粒大小、体积浓度和/或颗粒大小分布的情况调整对钻井系统的控制。

方法50的动作以功能模块的形式图示，图2所示的模块的先后顺序和模块中动作的划分并非限于图示的实施例。例如，某些模块可以按照不同的顺序进行；一个模块中的动作可以与另一个或多个模块中的动作组合，或拆分为多个模块。例如，计算R参数的计算值的步骤53和确定R参数的理论曲线的步骤54的先后顺序可以调换或同时进行。图1所示的超声探头11可以用来发射超声信号并接收回波信号。回波信号可被图1中的超声收发装置13、数字转换器17和处理器19进行处理产生颗粒大小、体积浓度和颗粒大小分布的信息，由显示器40显示。处理器19可用来实施步骤53-57。

虽然结合特定的实施方式对本发明进行了说明，但本领域的技术人员可以理解，对本发明可以作出许多修改和变型。因此，要认识到，权利要求书的意图在于涵盖在本发明真正构思和范围内的所有这些修改和变型。

Claims (10)

1.一种测量方法，其特征在于，所述测量方法包括：

发射若干超声信号至被测流体；

接收从被测流体中的颗粒反射或散射回来的回波信号；

根据所述回波信号计算与颗粒属性有关的R参数的计算值；

确定所述与颗粒属性有关的R参数的理论曲线；及

根据所述R参数的计算值和所述R参数的理论曲线确定流体中颗粒的平均直径；

根据所述颗粒的平均直径确定流体中颗粒的体积浓度。

2.如权利要求1所述的测量方法，其特征在于：所述测量方法进一步包括利用间隔排布在管道上的若干超声探头发射超声信号至被测流体的若干区域，确定每一所述区域内的颗粒的平均直径和体积浓度，且根据若干所述区域内的颗粒的平均直径和体积浓度确定管道区域内的颗粒大小分布。

3.如权利要求1所述的测量方法，其特征在于：所述确定颗粒的平均直径的步骤包括拟合所述R参数的计算值和所述R参数的理论曲线获得颗粒的平均直径。

4.如权利要求1所述的测量方法，其特征在于：所述计算R参数的计算值的步骤包括：

确定所述回波信号的电压幅度均方值关于颗粒属性的函数；及

根据所述函数计算所述R参数的计算值。

5.如权利要求1所述的测量方法，其特征在于：所述确定R参数的理论曲线的步骤包括根据颗粒的反向散射形态函数平均值和归一化总散射截面平均值确定所述R参数的理论曲线。

6.一种测量系统，其特征在于，所述测量系统包括：

超声探头，用来发射若干超声信号至被测流体且接收从被测流体中的颗粒反射或散射回来的回波信号；及

处理器，连接于所述超声探头，且用来根据所述回波信号计算与颗粒属性有关的R参数的计算值、确定所述与颗粒属性有关的R参数的理论曲线、根据所述R参数的计算值和所述R参数的理论曲线确定流体中颗粒的平均直径，且根据所述颗粒的平均直径确定流体中颗粒的体积浓度。

7.如权利要求6所述的测量系统，其特征在于：所述测量系统包括管道，若干所述超声探头间隔排布于所述管道上且用来发射超声信号至被测流体的若干区域，且所述处理器用来确定每一所述区域内的颗粒的平均直径和体积浓度，并根据若干所述区域内的颗粒的平均直径和体积浓度确定管道区域内的颗粒大小分布。

8.如权利要求6所述的测量系统，其特征在于：所述处理器用来拟合所述R参数的计算值和所述R参数的理论曲线获得颗粒的平均直径。

9.如权利要求6所述的测量系统，其特征在于：所述处理器用来确定所述回波信号的电压幅度均方值关于颗粒属性的函数，且根据所述函数计算所述R参数的计算值。

10.如权利要求6所述的测量系统，其特征在于：所述处理器用来根据颗粒的反向散射形态函数平均值和归一化总散射截面平均值确定所述R参数的理论曲线。