CN114814275A - 多相流体流速的动态计算方法及装置 - Google Patents
多相流体流速的动态计算方法及装置 Download PDFInfo
- Publication number
- CN114814275A CN114814275A CN202110108862.6A CN202110108862A CN114814275A CN 114814275 A CN114814275 A CN 114814275A CN 202110108862 A CN202110108862 A CN 202110108862A CN 114814275 A CN114814275 A CN 114814275A
- Authority
- CN
- China
- Prior art keywords
- echo
- amplitude
- data
- multiphase fluid
- relative amplitude
- 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.)
- Pending
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 134
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 52
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000010606 normalization Methods 0.000 claims abstract description 23
- 238000002592 echocardiography Methods 0.000 claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 17
- 230000008859 change Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 239000003921 oil Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 239000003129 oil well Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000004164 analytical calibration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种多相流体流速的动态计算方法及装置，其中该方法包括：接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；对回波串数据进行归一化处理，得到每一回波的相对幅值；根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速。本发明可自动采集回波串数据，降低了多相流体流速的计算误差；同时，也可实时采集回波串数据，实现了多相流体流速的自动化动态计算。
Description
技术领域
本发明涉及多相流体测量技术领域，尤其涉及多相流体流速的动态计算方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
包含多种组分的复杂多相流体在工业领域十分常见，其准确计量一直是困扰工业界的一大难题。多相流体在线检测系统在这一背景下发展起来。核磁共振多相流量计是多相流体在线检测系统中的新成员，其基于核磁共振技术实现对石油工业生产油、气、水、固多相流体在线无损检测。
核磁共振多相流量计的探测方式是通过分别开展相含率(含油率、含水率、含气率、含固率)测量和管内流速(平均流速或分层流速)测量，最终得到各相流量。
与其他多相流量计相同的是，为了实现核磁共振多相流量计的精确测量，使用前的仪器标定是必不可少的，即针对目标使用油井的流体进行采样测量，根据平均流量、流体性质和平均含水率等测量结果，调整仪器参数，使其测量精度最高。但实际使用时，往往会遇到很多问题，诸如不同地质区块的石油井采出多相流体的流体性质差异较大，每测量一口油井就需要进行一次标定工作，增加了现场作业量和时间；此外，即使是同一口油井，在其整个生命周期内，往往流体性质、含量、组分均可能发生较大变化，这意味着隔一段时间就需要对设备重新标定。
多相流体相含率的测量通常基于T2谱或T1谱峰的位置区分不同组分，通常只需要对不同组分的标准样品进行标定后，实际使用过程中不需要进行再次标定。但多相流体流速标定则不同，其测量原理是基于回波串的衰减速率来计算流速，由于回波串本身受多重因素的影响，即使流速为零也存在非线性衰减。
目前，国内外针对核磁共振多相流量计的流速测量，均采用隔段时间重新标定的方式。需要依靠人工手动选择线性度较好的一段回波串(卡一个时间段内的回波串)进行计算；在实际应用时，时间段一旦选好后就不再改变。然而，不幸的是，这一时间段的选取针对性极强，流体性质、含量、组分的变化均会影响该段回波串的线性度。因此传统的核磁共振多相流量计流速测量方法的参数，无法根据流体变化动态调整，最终造成多相流体流速的测量误差变大。
发明内容
本发明实施例提供一种多相流体流速的动态计算方法，用以实现多相流体流速的动态计算，减小多相流体流速计算误差，该方法包括：
接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；
对回波串数据进行归一化处理，得到每一回波的相对幅值；
根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；
将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速。
本发明实施例还提供一种多相流体流速的动态计算装置，用以实现多相流体流速的动态计算，减小多相流体流速计算误差，该装置包括：
数据采集模块，用于通过核磁共振波谱仪，实时采集预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；
归一化处理模块，用于对回波串数据进行归一化处理，得到每一回波的相对幅值；所述归一化处理用于将回波串数据中每一回波的幅值与首个回波的幅值进行比较；
线性拟合处理模块，用于根据每一回波的相对幅值，对回波串数据进行，得到回波串的线性拟合数据；
流速计算模块，用于以核磁共振多相流量计，根据回波串的线性拟合数据，实时计算多相流体流速。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述多相流体流速的动态计算方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述多相流体流速的动态计算方法的计算机程序。
本发明实施例中，接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；对回波串数据进行归一化处理，得到每一回波的相对幅值；根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速，与现有技术对比，可自动采集回波串数据，不需人工手动选择回波串数据，也避免了手动选取时间段对应的回波串数据的线性度易发生改变，导致多相流体流速计算误差较大的问题，降低了多相流体流速的计算误差；同时，也可实时采集回波串数据，实现了多相流体流速的自动化动态计算。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中一种多相流体流速的动态计算方法的流程示意图；
图2为本发明实施例中多相流体在不同流速下所采集得到的回波串数据的示意图；
图3为本发明实施例中一种多相流体流速的动态计算方法的具体示例图；
图4为本发明实施例中应用一种多相流体流速的动态计算方法的实例示意图；
图5为本发明实施例中应用一种多相流体流速的动态计算方法的实例示意图；
图6为本发明实施例中应用一种多相流体流速的动态计算方法的实例示意图；
图7为本发明实施例中一种多相流体流速的动态计算装置的结构示意图；
图8为本发明实施例中一种多相流体流速的动态计算装置的具体示例图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1为本发明实施例中一种多相流体流速的动态计算方法的流程示意图，本发明实施例提供一种多相流体流速的动态计算方法，可以包括：
步骤101：接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；
步骤102：对回波串数据进行归一化处理，得到每一回波的相对幅值；
步骤103：根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；
步骤104：将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速。
本发明实施例中，接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；对回波串数据进行归一化处理，得到每一回波的相对幅值；根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速，与现有技术对比，可自动采集回波串数据，不需人工手动选择回波串数据，也避免了手动选取时间段对应的回波串数据的线性度易发生改变，导致多相流体流速计算误差较大的问题，降低了多相流体流速的计算误差；同时，也可实时采集回波串数据，实现了多相流体流速的自动化动态计算。
具体实施时，首先接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值。
实施例中，对于连续流动的油气水多相流体，不同流速、不同含水率的流体，其曲线衰减速率不同，体现在回波串直线段的斜率不同。图2所示是相同油、气、水、固比例多相流体在不同流速下所采集得到的回波串数据，图2中曲线均是油水两相混合流体的回波串数据示意图，其中衰减最慢的(是指图2最上面的曲线)曲线是100％纯水样品，往下依次是指含水率为90％、80％、70％.....10％的油水两相混合流体。从图2可以看出，流速越快越高，回波串衰减速率越大，在预设的回波串采集时间内所采集的回波串的线性度就越差。
实际测试时发现，采用不同含水率的油、水两相流体，在同一流速下采集回波串数据，随着含水率的下降，回波串的衰减速率变大，线性度越差，呈现和图2类似的现象。因此，参与线性拟合的回波串数据的选取对于流速测量结果的准确性有直接影响。原则是选取线性度好的回波串数据，并在此前提下尽可能多的选取回波数据。
在上述实施例中，可实时采集的预设时长内多相流体在不同流速下的回波串数据，保证了采集的回波串数据的线性度较好，从而降低了多相流体流速的计算误差。
具体实施时，在接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据后，对回波串数据进行归一化处理，得到每一回波的相对幅值。
实施例中，对回波串数据进行归一化处理，得到每一回波的相对幅值，可以包括：将回波串数据中每一回波的幅值除以首个回波的幅值，得到每一回波的幅值与首个回波幅值的比值；将每一回波的幅值与首个回波的幅值的比值设置为每一回波的相对幅值。
在上述实施例中，通过进行归一化处理，可有助于在后续步骤中根据每一回波的相对幅值，对回波串数据进行线性拟合处理。
具体实施时，在对回波串数据进行归一化处理，得到每一回波的相对幅值后，根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据。
实施例中，根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据，可以包括：在回波串数据中最后一个回波的相对幅值大于第一预设相对幅值时，对回波串数据中相对幅值大于第二预设相对幅值的全部回波，进行线性拟合处理，得到第一线性拟合数据；所述第二预设相对幅值大于第一预设相对幅值；在回波串数据中最后一个回波的相对幅值小于等于第一预设幅值时，对回波串数据中，相对幅值大于等于第一预设幅值且小于等于第二预设相对幅值的全部回波，进行线性拟合处理，得到第二线性拟合数据。
在上述实施例中，第一预设相对幅值和第二预设相对幅值可根据多相流体所处工区的实际情况设置，且可根据工作人员的实际使用需求随时修改。
举一例，图3为本发明实施例中一种多相流体流速的动态计算方法的具体示例图，如图3所示，第一预设相对幅值可设置为0.6；第二预设相对幅值可设置为0.95。在该实例中，首先可采集50ms的回波串数据，并对回波串做归一化处理(所有回波幅值除以首个回波的幅值)；然后判断最后一个回波的幅值是否大于0.6，若是，说明回波串衰减的速度慢，则选取幅值<0.95的所有回波串进行线性拟合(>＝0.95的回波数据不稳定，幅值变化较大，线性度不好)；若判断最后一个回波的幅值是否小于等于0.6，则确定回波串在小于0.6的部分线性度较差，不能参与拟合，故选取幅值小于等于0.95，且大于等于0.6的所有回波串进行线性拟合。
具体实施时，在根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据后，将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速。
在上述实施例中，核磁共振多相流量计的流速测量原理是回波串衰减速率正比于流速，因此可选取线性度较好，衰减主要受流速影响(回波串本身的弛豫衰减不明显)的回波串数据段，进行线性拟合，就可以量化计算出多相流体的流速。本发明实施例弥补了现有技术下核磁共振多相流量计流速测量阶段无法根据流体变化动态自动调节参数的问题，实现了核磁共振多相流量测量全过程参数自调节，可助力于多相流体流量的准确检测。
具体实施时，本发明实施例提供的多相流体流速的动态计算方法还可以包括：根据每一回波的相对幅值，确定多相流体的工况属性；所述多相流体的工况属性包括多相流体的回波串衰减速度、流速和含水率。
实施例中，可根据每一回波的相对幅值，对多相流体的工况属性进行初步判断，实现对多相流体的工况的初步概括。
举一例，如图3所示，可根据判断回波串数据中最后一个回波的幅值是否大于0.6，若是，说明回波串衰减的速度慢，此时通常对应慢流速且含水率较高的工况；若判断最后一个回波的幅值是否小于等于0.6，此时通常对应快流速或含水率低的工况。
下面给出一个具体实施例，来说明本发明的方法的具体应用。图3为本发明实施例中应用一种多相流体流速的动态计算方法的实例示意图，如图3所示，该实施例中，可以包括如下步骤：
首先切换阀门至核磁共振多相流量计的流动测量模式，待测流体会连续流动并通过NMR(Nuclear Magnetic Resonance，核磁共振)探头(如图6所示)，此时可以通过回波串衰减速度测流速，根据回波串衰减速度有两种测量方法供选择，若衰减较慢(即判断最后一个回波的幅值是否大于0.6)，则选取卡回波串衰减幅度的方式截取一段直线段回波串进行拟合(如图4所示)；若衰减较快(即判断最后一个回波的幅值是小于等于0.6)，则选取卡回波串测量时间的方式截取一段直线段回波串进行拟合(如图5所示)，其方法选择原则是，尽量选择回波串里面线性度较好的一个直线段进行斜率拟合。
本发明实施例中，接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；对回波串数据进行归一化处理，得到每一回波的相对幅值；根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速，与现有技术对比，可自动采集回波串数据，不需人工手动选择回波串数据，也避免了手动选取时间段对应的回波串数据的线性度易发生改变，导致多相流体流速计算误差较大的问题，降低了多相流体流速的计算误差；同时，也可实时采集回波串数据，实现了多相流体流速的自动化动态计算。
本发明实施例属于低场核磁共振应用方法领域，可在油田现场实现利用核磁共振方法，在线计量油、气、水、固多相流流量的准确测量，解决了现有技术下在流速测量阶段的回波串选取参数，无法根据流体变化动态调整，从而造成误差较大的问题，因此本发明实施例中回波串动态选取方法，可同时通过时间和回波串衰减率选取回波串，可由核磁共振波谱仪自动完成无需人员干预。本发明实施例针对国际公认的多相流在线准确计量难题和核磁共振多相流量计流速测量参数无法根据流体变化动态调整的问题，提出了一种可用于流体流速测量的回波串动态处理方法，进一步提升其智能性和准确度，并可有效减少现场作业次数。
本发明实施例提供的多相流体流速的动态计算方法，可形成一套基于多维低场核磁共振流量装置的流速计算方法及刻度流程，从而可将核磁共振测量技术应用在具有流动流体计量需求的工业领域，进而可形成安装在工业应用现场的复杂流体流量计量设备，有助于多相流体流速的自动化动态计算。
本发明实施例中还提供了一种多相流体流速的动态计算装置，如下面的实施例所述。由于该装置解决问题的原理与多相流体流速的动态计算方法相似，因此该装置的实施可以参见多相流体流速的动态计算方法的实施，重复之处不再赘述。
图7为本发明实施例中一种多相流体流速的动态计算装置的结构示意图，本发明实施例提供的多相流体流速的动态计算装置，可以包括：
数据采集模块01，用于通过核磁共振波谱仪，实时采集预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；
归一化处理模块02，用于对回波串数据进行归一化处理，得到每一回波的相对幅值；所述归一化处理用于将回波串数据中每一回波的幅值与首个回波的幅值进行比较；
线性拟合处理模块03，用于根据每一回波的相对幅值，对回波串数据进行，得到回波串的线性拟合数据；
流速计算模块04，用于以核磁共振多相流量计，根据回波串的线性拟合数据，实时计算多相流体流速。
在一个实施例中，归一化处理模块，具体用于：将回波串数据中每一回波的幅值除以首个回波的幅值，得到每一回波的幅值与首个回波幅值的比值；
将每一回波的幅值与首个回波的幅值的比值设置为每一回波的相对幅值。
在一个实施例中，线性拟合处理模块，具体用于：
在回波串数据中最后一个回波的相对幅值大于第一预设相对幅值时，对回波串数据中相对幅值大于第二预设相对幅值的全部回波，进行线性拟合处理，得到第一线性拟合数据；所述第二预设相对幅值大于第一预设相对幅值；
在回波串数据中最后一个回波的相对幅值小于等于第一预设幅值时，对回波串数据中，相对幅值大于等于第一预设幅值且小于等于第二预设相对幅值的全部回波，进行线性拟合处理，得到第二线性拟合数据。
在一个实施例中，如图8所示，还可以包括：工况属性确定模块05，用于：根据每一回波的相对幅值，确定多相流体的工况属性；所述多相流体的工况属性包括多相流体的回波串衰减速度、流速和含水率。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述多相流体流速的动态计算方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述多相流体流速的动态计算方法的计算机程序。
本发明实施例中，接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；对回波串数据进行归一化处理，得到每一回波的相对幅值；根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速，与现有技术对比，可自动采集回波串数据，不需人工手动选择回波串数据，也避免了手动选取时间段对应的回波串数据的线性度易发生改变，导致多相流体流速计算误差较大的问题，降低了多相流体流速的计算误差；同时，也可实时采集回波串数据，实现了多相流体流速的自动化动态计算。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种多相流体流速的动态计算方法，其特征在于，包括：
接收核磁共振波谱仪实时采集的预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；
对回波串数据进行归一化处理，得到每一回波的相对幅值；
根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据；
将回波串的线性拟合数据提供给核磁共振多相流量计，以使核磁共振多相流量计根据回波串的线性拟合数据实时计算多相流体流速。
2.如权利要求1所述的方法，其特征在于，对回波串数据进行归一化处理，得到每一回波的相对幅值，包括：
将回波串数据中每一回波的幅值除以首个回波的幅值，得到每一回波的幅值与首个回波幅值的比值；
将每一回波的幅值与首个回波的幅值的比值设置为每一回波的相对幅值。
3.如权利要求1所述的方法，其特征在于，根据每一回波的相对幅值，对回波串数据进行线性拟合处理，得到回波串的线性拟合数据，包括：
在回波串数据中最后一个回波的相对幅值大于第一预设相对幅值时，对回波串数据中相对幅值大于第二预设相对幅值的全部回波，进行线性拟合处理，得到第一线性拟合数据；所述第二预设相对幅值大于第一预设相对幅值；
在回波串数据中最后一个回波的相对幅值小于等于第一预设幅值时，对回波串数据中，相对幅值大于等于第一预设幅值且小于等于第二预设相对幅值的全部回波，进行线性拟合处理，得到第二线性拟合数据。
4.如权利要求1所述的方法，其特征在于，还包括：根据每一回波的相对幅值，确定多相流体的工况属性；所述多相流体的工况属性包括多相流体的回波串衰减速度、流速和含水率。
5.一种多相流体流速的动态计算装置，其特征在于，包括：
数据采集模块，用于通过核磁共振波谱仪，实时采集预设时长内多相流体在不同流速下的回波串数据；所述回波串数据包括以时间顺序排列的多相流体在不同流速下的多个回波的幅值；
归一化处理模块，用于对回波串数据进行归一化处理，得到每一回波的相对幅值；所述归一化处理用于将回波串数据中每一回波的幅值与首个回波的幅值进行比较；
线性拟合处理模块，用于根据每一回波的相对幅值，对回波串数据进行，得到回波串的线性拟合数据；
流速计算模块，用于以核磁共振多相流量计，根据回波串的线性拟合数据，实时计算多相流体流速。
6.如权利要求5所述的装置，其特征在于，归一化处理模块，具体用于：将回波串数据中每一回波的幅值除以首个回波的幅值，得到每一回波的幅值与首个回波幅值的比值；
将每一回波的幅值与首个回波的幅值的比值设置为每一回波的相对幅值。
7.如权利要求5所述的装置，其特征在于，线性拟合处理模块，具体用于：
在回波串数据中最后一个回波的相对幅值大于第一预设相对幅值时，对回波串数据中相对幅值大于第二预设相对幅值的全部回波，进行线性拟合处理，得到第一线性拟合数据；所述第二预设相对幅值大于第一预设相对幅值；
在回波串数据中最后一个回波的相对幅值小于等于第一预设幅值时，对回波串数据中，相对幅值大于等于第一预设幅值且小于等于第二预设相对幅值的全部回波，进行线性拟合处理，得到第二线性拟合数据。
8.如权利要求5所述的装置，其特征在于，还包括：工况属性确定模块，用于：根据每一回波的相对幅值，确定多相流体的工况属性；所述多相流体的工况属性包括多相流体的回波串衰减速度、流速和含水率。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110108862.6A CN114814275A (zh) | 2021-01-27 | 2021-01-27 | 多相流体流速的动态计算方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110108862.6A CN114814275A (zh) | 2021-01-27 | 2021-01-27 | 多相流体流速的动态计算方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114814275A true CN114814275A (zh) | 2022-07-29 |
Family
ID=82524782
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110108862.6A Pending CN114814275A (zh) | 2021-01-27 | 2021-01-27 | 多相流体流速的动态计算方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114814275A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115902170B (zh) * | 2022-10-31 | 2023-08-18 | 江苏麦赫物联网科技有限公司 | 多相流流态工况分析方法和设备 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS59182316A (ja) * | 1983-01-27 | 1984-10-17 | シ−メンス・アクチエンゲゼルシャフト | 流速測定装置 |
| AU8167698A (en) * | 1997-06-24 | 1999-01-04 | Southwest Research Institute | Measurement of flow fractions, flow velocities, and flow rates of a multiphase fluid using nmr sensing |
| US20120092006A1 (en) * | 2010-10-19 | 2012-04-19 | Baker Hughes Incorporated | Multiphase Flow Measurement Using Nuclear Magnetic Resonance |
| CN106706050A (zh) * | 2015-08-10 | 2017-05-24 | 杭州思筑智能设备有限公司 | 一种利用超声波流量计测量气体流量的方法 |
| CN107942398A (zh) * | 2017-11-07 | 2018-04-20 | 中国石油天然气股份有限公司 | 一种确定多相流体流速的方法及装置 |
| US20190086249A1 (en) * | 2017-09-19 | 2019-03-21 | Petrochina Company Limited | Method and apparatus for determining flow rates of components of multiphase fluid |
-
2021
- 2021-01-27 CN CN202110108862.6A patent/CN114814275A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS59182316A (ja) * | 1983-01-27 | 1984-10-17 | シ−メンス・アクチエンゲゼルシャフト | 流速測定装置 |
| AU8167698A (en) * | 1997-06-24 | 1999-01-04 | Southwest Research Institute | Measurement of flow fractions, flow velocities, and flow rates of a multiphase fluid using nmr sensing |
| US20120092006A1 (en) * | 2010-10-19 | 2012-04-19 | Baker Hughes Incorporated | Multiphase Flow Measurement Using Nuclear Magnetic Resonance |
| CN106706050A (zh) * | 2015-08-10 | 2017-05-24 | 杭州思筑智能设备有限公司 | 一种利用超声波流量计测量气体流量的方法 |
| US20190086249A1 (en) * | 2017-09-19 | 2019-03-21 | Petrochina Company Limited | Method and apparatus for determining flow rates of components of multiphase fluid |
| CN107942398A (zh) * | 2017-11-07 | 2018-04-20 | 中国石油天然气股份有限公司 | 一种确定多相流体流速的方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 李利品等: "Halbach阵列核磁共振多相流传感器参数优化与仿真", 仪表技术与传感器, no. 02, 15 February 2020 (2020-02-15), pages 5 - 9 * |
| 李玉星等: "多相流流速和相分率测试技术研究进展", 管道技术与设备, no. 03, 20 June 1999 (1999-06-20), pages 6 - 9 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115902170B (zh) * | 2022-10-31 | 2023-08-18 | 江苏麦赫物联网科技有限公司 | 多相流流态工况分析方法和设备 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105804734A (zh) | 一种利用核磁共振测井识别稠油储层的方法 | |
| CN103529237B (zh) | 一种泥沙群体沉速的测量方法及测量装置 | |
| CN107525553B (zh) | 一种确定多相流体组分流量的方法及装置 | |
| CN110346005B (zh) | 基于深度学习的科里奥利质量流量计数字信号处理方法 | |
| AU2013206724A1 (en) | Nuclear magnetic flowmeter | |
| Carlander et al. | Installation effects on an ultrasonic flow meter with implications for self diagnostics | |
| CN114814275A (zh) | 多相流体流速的动态计算方法及装置 | |
| CN115824324A (zh) | 一种基于旋进旋涡mems检测技术的湿气流量测量方法 | |
| EP2504696A1 (en) | Composition of multiphase flow | |
| CN114137015A (zh) | 孔隙度校正方法及装置 | |
| RU2386930C2 (ru) | Способ определения параметров потока многофазной смеси жидкости и газа | |
| CN114165216A (zh) | 一种单井连续多相流量计量系统、存储介质和计算机设备 | |
| CN105571661A (zh) | 基于示值误差拟合的电磁流量计仪表特征系数的计算方法 | |
| Thiers et al. | Use of a Computer Program to Correct for Sample Interaction: A Significant Adjunct to Continuous-Flow Analysis | |
| CN211291499U (zh) | 一种流体输出过程的多组分体积配比的实时在线检测系统 | |
| CN103498658A (zh) | 一种伽马测井数据处理方法 | |
| CN112683347B (zh) | 流体各组分流量测量装置及方法 | |
| Ingram et al. | An automated instrumentation system for flow and loss measurements in a cascade | |
| CN108760773B (zh) | 在线碳氢化合物含水率及纯油密度分析撬座系统 | |
| NO851809L (no) | Fremgangsmaate for bestemmelse av stroemningshastigheten til en flerfasestroem. | |
| CN114739486A (zh) | 多相流核磁共振流量计的刻度方法 | |
| CN114624274A (zh) | 多相流相含率的检测方法 | |
| CN115478835A (zh) | 一种井下流体测量模型获取方法及装置 | |
| CN115902170B (zh) | 多相流流态工况分析方法和设备 | |
| RU2390757C1 (ru) | Способ определения вязкости жидкости |
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 |