CN115479867B - 一种非牛顿流体快速判别方法 - Google Patents

一种非牛顿流体快速判别方法 Download PDF

Info

Publication number
CN115479867B
CN115479867B CN202210923445.1A CN202210923445A CN115479867B CN 115479867 B CN115479867 B CN 115479867B CN 202210923445 A CN202210923445 A CN 202210923445A CN 115479867 B CN115479867 B CN 115479867B
Authority
CN
China
Prior art keywords
fluid
viscosity
measured
characteristic parameter
frequency
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
CN202210923445.1A
Other languages
English (en)
Other versions
CN115479867A (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of 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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202210923445.1A priority Critical patent/CN115479867B/zh
Publication of CN115479867A publication Critical patent/CN115479867A/zh
Application granted granted Critical
Publication of CN115479867B publication Critical patent/CN115479867B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • G01N2011/0026—Investigating specific flow properties of non-Newtonian fluids

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种非牛顿流体快速判别方法，该方法包括以下步骤：获取待检测流体的谐振频率和反谐振频率；绘制反谐振频率与黏度非线性特性参数组合之间的关系曲线；绘制谐振频率与黏度非线性特性参数组合之间的关系曲线；通过绘制得到的反谐振频率与黏度非线性特性参数组合之间的关系曲线和谐振频率与黏度非线性特性参数组合之间的关系曲线进行比较，完成流体性质的判别。本发明的方法为基于多特征频率参数对其流体性质进行判别，仅需要通过单次对待测流体进行谐振频率和反谐振频率的测量便可实现，测量速度快、测量稳定性高、能够有效快速地实现非牛顿流体的测量，有效解决了现有技术中对采用单一参数对非牛顿流体的判别导致准确率较低的问题。

Description

一种非牛顿流体快速判别方法
技术领域
本发明涉及流体性质判别领域，尤其涉及一种非牛顿流体快速判别方法。
背景技术
牛顿流体是指牛顿1687年提出的一种理想粘性液体。即指具有层流特征的流体，相邻的两层平行流动的液体间产生的剪切应力与垂直于流动方向的速度梯度成正比时，这种液体即为牛顿流体。
非牛顿流体，是指不满足牛顿黏性实验定律的流体，即其剪应力与剪切应变率之间不是线性关系的流体。常见的具有简单结构的液体是牛顿液(如水，甘油，油等)。然而，具有复杂结构的流体(聚合物熔体或溶液、悬浮液等)通常是非牛顿流体。黏度作为流体最重要的物理特性之一，其特性的表征对于工程应用和基础研究都极为重要。传统的流体黏度测量方法有毛细管法、旋转法、落体法、平板法、黏度杯法等，这些传统的黏度测量方法耗时长，不能实时在线地获得流体的黏度特性，且所需要的样品量一般较大，不能满足小样品量的测量要求。尽管黏度及其非线性问题的提出已经有数百年的时间，期间研究者提出了各式各样的模型来描述非线性，例如剪切增稠液、剪切稀化液、宾汉塑性体、结构流体等等，其中又有多种模型来描述，如Carreau model,Ellis model,Sisko model,Power law model,Carreau-Yasuda model,Meter model等。但每一个模型都有其局限性，仅能在比较有限的范围内来描述黏度特性。
目前，非牛顿流体与牛顿流体之间的关系主要用黏度是否等于恒定值用来判断流体的性质。牛顿流体能够在剪切应力参数增大的过程中，其黏度依然保持一个恒定数值，例如：水、油等；非牛顿流体的流体切应力变化性质目前还没有深入的研究，随着剪切应力的增大，黏度与剪切应力呈现非线性的关系，至于这种非线性关系与牛顿流体之间的线性关系之间的联系还没有进行深入探索。
目前，已有的非牛顿流体检测方法是只利用一种参数对流体性质进行判别。例如：使用声波传感器仅对待测流体的单一特征频率进行测量，用于代表牛顿流体的黏度特性从而进行判别；使用磁致伸缩条仅对待测流体的谐振频率进行测量，最终得到单一参数代表牛顿流体的黏度特性进行判别，其具有判别时间长、判别工序繁多、对操作人员要求高、计算量大等缺点。而且单一参数不适用于非牛顿流体的快速判别，造成检测的准确率较低。
发明内容
本发明的目的在于提出一种非牛顿流体快速判别方法，以解决上述背景技术中存在的一个或多个技术问题。
为达此目的，本发明采用以下技术方案：
一种非牛顿流体快速判别方法，该方法包括以下步骤：
获取待检测流体的谐振频率和反谐振频率；
以反谐振频率fr′作为横坐标，以黏度非线性特性参数组合作为纵坐标来绘制反谐振频率与黏度非线性特性参数组合之间的关系曲线；
以谐振频率fr作为横坐标，以黏度非线性特性参数组合作为纵坐标来绘制谐振频率与黏度非线性特性参数组合之间的关系曲线；
通过绘制得到的反谐振频率fr′与黏度非线性特性参数组合之间的关系曲线和谐振频率fr与黏度非线性特性参数组合/>之间的关系曲线进行比较，完成流体性质的判别。
优选的，该方法还包括获取待测流体剪切力与黏度随着时间的变化曲线。
优选的，待测流体剪切力与黏度随着时间的变化曲线中，若待测流体随着剪切应力地增大，其黏度特性无发生任何变化，则该待测流体为牛顿流体；若待测流体随着剪切应力的增大，其黏度展现出正相关或负相关的非线性变化趋势，则该待测流体为非牛顿流体。
优选的，反谐振频率fr′与黏度非线性特性参数组合之间的关系曲线中，若待测流体的黏度非线性特性参数组合/>与反谐振频率fr′之间为正相关的线性关系，则该待测流体为牛顿流体；若待测流体黏度非线性特性参数组合/>与反谐振频率fr′之间不存在线性关联，则该待测流体为非牛顿流体。
优选的，谐振频率fr与黏度非线性特性参数组合之间的关系曲线中，若待测流体的谐振频率fr与黏度非线性特性参数组合/>之间为正相关的线性关系，则该待测流体为牛顿流体；若待测流体的谐振频率fr与黏度非线性特性组合/>的关系并非线性相关关系，且其变化曲线的曲率大于0，其二次导数从小于0向大于0的趋势进行过渡，并具有指数式增加的趋势，则该待测流体为非牛顿流体。
优选的，待测流体的谐振频率和反谐振频率通过阻抗分析仪获取。
本发明的有益效果为：本发明的方法为基于多特征频率参数对其流体性质进行判别，仅需要通过单次对待测流体进行谐振频率和反谐振频率的测量便可实现，测量速度快、测量稳定性高、测量成本低、求测量人员要求不高、能够有效快速地实现非牛顿流体的测量，有效解决了现有技术中对采用单一参数对非牛顿流体的判别导致准确率较低的问题。
附图说明
附图对本发明做进一步说明，但附图中的内容不构成对本发明的任何限制。
图1是非牛顿液体中谐振频率，反谐振频率，实验共振频率之间的关系曲线；
图2是不同性质流体的黏度与剪切速率之间的变化曲线；
图3是不同流体中的反谐振频率fr′与黏度非线性特性参数组合的关系曲线；
图4是不同流体中的谐振频率fr黏度非线性特性参数组合的理论关系曲线
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。
本实施例的一种非牛顿流体快速判别方法该方法包括以下步骤：
获取待检测流体的谐振频率和反谐振频率；
以反谐振频率fr′作为横坐标，以黏度非线性特性参数组合作为纵坐标来绘制反谐振频率与黏度非线性特性参数组合之间的关系曲线；
以谐振频率fr作为横坐标，以黏度非线性特性参数组合作为纵坐标来绘制谐振频率与黏度非线性特性参数组合之间的关系曲线；
通过绘制得到的反谐振频率fr′与黏度非线性特性参数组合之间的关系曲线和谐振频率fr与黏度非线性特性参数组合/>之间的关系曲线进行比较，完成流体性质的判别。
由此，本实施例的方法为基于多特征频率参数对其流体性质进行判别，仅需要通过单次对待测流体进行谐振频率和反谐振频率的测量便可实现，测量速度快、测量稳定性高、测量成本低、求测量人员要求不高、能够有效快速地实现非牛顿流体的测量，有效解决了现有技术中对采用单一参数对非牛顿流体的判别导致准确率较低的问题。
若待测流体为牛顿流体，由于黏度与剪切速率无关，则在测试的过程中只需要用到一种实验参数，例如谐振频率fr足以用于描述液体黏度。
若待测流体为非牛顿流体，由于黏度随剪切速率和速度而变化。则在本发明技术方案中需要引入额外的参数来呈现这个非线性。本发明中引入谐振频率fr，反谐振频率fr′，实验共振频率f0三种参数进行描述。其中，谐振频率fr通常被选择在最大信号点，反谐振频率fr′被选择在在最小信号点处，实验共振频率f0被选择在这两点之间的某处被拾取。以上三种频率分布在一个相当宽的范围内，在该区间内位移、速度和剪切速率变化很大，如图1所示。在这个范围内，非牛顿流体黏度的非线性特性也非常明显，因此可利用这几种频率的特殊函数组合来表征黏度的非线性特性，从而快速识别待测流体的性质。
优选的，该方法还包括获取待测流体剪切力与黏度随着时间的变化曲线。
进一步的，待测流体剪切力与黏度随着时间的变化曲线中，若待测流体随着剪切应力地增大，其黏度特性无发生任何变化，则该待测流体为牛顿流体；若待测流体随着剪切应力的增大，其黏度展现出正相关或负相关的非线性变化趋势，则该待测流体为非牛顿流体。参考附图2，图中的Pseudoplastic，Bingham Plastic，Structural，Dilatant均为经典的非牛顿流体黏度变化曲线。从图中可以看出代表常见的非牛顿流体的Bingham Plastic，Structural曲线在一开始随着剪切应力增大的过程中，其黏度会在短时间内维持不变；随着剪切应力的继续增大，其黏度会持续下降；但是最终当剪切应力增大到一定的程度时，最后黏度的变化趋向于平稳，总体呈现出“阶梯型”的变化趋势。因此，该特征可以作为判别牛顿流体和非牛顿流体的参考。
优选的，反谐振频率fr′与黏度非线性特性参数组合之间的关系曲线中，若待测流体的黏度非线性特性参数组合/>与反谐振频率fr之间为正相关的线性关系，则该待测流体为牛顿流体；若待测流体黏度非线性特性参数组合/>与反谐振频率fr之间不存在线性关联，则该待测流体为非牛顿流体。如图3中，可以从fr′与黏度非线性特性参数组合/>之间的关系进行分析，牛顿流体与BP2380(一种参考牛顿流体)的黏度非线性特性参数与反谐振频率之间呈现正相关的线性关系；而对于非牛顿流体而言，黏度非线性特性参数与反谐振频率之间线性关联较小，在该区间内有明显的非线性特性。
优选的，谐振频率fr与黏度非线性特性参数组合之间的关系曲线中，若待测流体的谐振频率fr与黏度非线性特性参数组合/>之间为正相关的线性关系，则该待测流体为牛顿流体；若待测流体的谐振频率fr与黏度非线性特性组合/>的关系并非线性相关关系，且其变化曲线的曲率大于0，其二次导数从小于0向大于0的趋势进行过渡，并具有指数式增加的趋势，则该待测流体为非牛顿流体。如图4，根据其变化趋势以及对其曲率的分析，可以用于直观的判别待测流体的性质。
优选的，待测流体的谐振频率和反谐振频率通过阻抗分析仪获取。由此，通过通过阻抗分析仪单次对待测流体进行测量便可实现，利用阻抗分析仪的快速测量与快速计算能力，通过绘制谐振频率fr，反谐振频率fr′与代表黏度非线性特性参数组合之间的关系曲线与牛顿流体的经验曲线变化趋势进行比较，快速的实现非牛顿流体的判别。
以上结合具体实施例描述了本发明的技术原理。这些描述只是为了解释本发明的原理，而不能以任何方式解释为对本发明保护范围的限制。基于此处的解释，本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式，这些方式都将落入本发明的保护范围之内。

Claims (4)

1.一种非牛顿流体快速判别方法，其特征在于，该方法包括以下步骤：
获取待检测流体的谐振频率和反谐振频率；
以反谐振频率作为横坐标，以黏度非线性特性参数组合/>作为纵坐标来绘制反谐振频率与黏度非线性特性参数组合之间的关系曲线；
以谐振频率作为横坐标，以黏度非线性特性参数组合/>作为纵坐标来绘制谐振频率与黏度非线性特性参数组合之间的关系曲线；
通过绘制得到的反谐振频率与黏度非线性特性参数组合/>之间的关系曲线和谐振频率/>与黏度非线性特性参数组合/>之间的关系曲线进行比较，完成流体性质的判别；
反谐振频率与黏度非线性特性参数组合/>之间的关系曲线中，若待测流体的黏度非线性特性参数组合/>与反谐振频率/>之间为正相关的线性关系，则该待测流体为牛顿流体；若待测流体黏度非线性特性参数组合/>与反谐振频率/>之间不存在线性关联，则该待测流体为非牛顿流体；
谐振频率与黏度非线性特性参数组合/>之间的关系曲线中，若待测流体的谐振频率/>与黏度非线性特性参数组合/>之间为正相关的线性关系，则该待测流体为牛顿流体；若待测流体的谐振频率/>与黏度非线性特性组合/>的关系并非线性相关关系，且其变化曲线的曲率大于0，其二次导数从小于0向大于0的趋势进行过渡，并具有指数式增加的趋势，则该待测流体为非牛顿流体。
2.根据权利要求1所述的一种非牛顿流体快速判别方法，其特征在于，该方法还包括获取待测流体剪切力与黏度随着时间的变化曲线。
3.根据权利要求2所述的一种非牛顿流体快速判别方法，其特征在于，待测流体剪切力与黏度随着时间的变化曲线中，若待测流体随着剪切应力地增大，其黏度特性无发生任何变化，则该待测流体为牛顿流体；若待测流体随着剪切应力的增大，其黏度展现出正相关或负相关的非线性变化趋势，则该待测流体为非牛顿流体。
4.根据权利要求1所述的一种非牛顿流体快速判别方法，其特征在于，待测流体的谐振频率和反谐振频率通过阻抗分析仪获取。
CN202210923445.1A 2022-08-02 2022-08-02 一种非牛顿流体快速判别方法 Active CN115479867B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210923445.1A CN115479867B (zh) 2022-08-02 2022-08-02 一种非牛顿流体快速判别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210923445.1A CN115479867B (zh) 2022-08-02 2022-08-02 一种非牛顿流体快速判别方法

Publications (2)

Publication Number Publication Date
CN115479867A CN115479867A (zh) 2022-12-16
CN115479867B true CN115479867B (zh) 2024-05-28

Family

ID=84422140

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210923445.1A Active CN115479867B (zh) 2022-08-02 2022-08-02 一种非牛顿流体快速判别方法

Country Status (1)

Country Link
CN (1) CN115479867B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105004633A (zh) * 2015-06-24 2015-10-28 广东工业大学 基于纳米线垂直阵列的流体粘度检测装置及检测方法
CN106290065A (zh) * 2016-07-28 2017-01-04 北京建筑大学 一种确定非牛顿流体在反应器内搅拌死区的方法及装置
CN107271332A (zh) * 2017-07-04 2017-10-20 西安交通大学 一种基于面内谐振的mems流体黏度传感器芯片及其制备方法
CN108871200A (zh) * 2018-06-01 2018-11-23 广东工业大学 一种基于表面纳米轴向光子结构回音壁微腔的探针式微位移传感系统
CN110631959A (zh) * 2019-09-30 2019-12-31 南方科技大学 测量非牛顿流体表观粘度的端面效应误差消除方法及装置
CN111537396A (zh) * 2020-04-24 2020-08-14 西安交通大学 一种多操作模式的压电黏度传感器芯片及其工作方法和制备方法
CN112729152A (zh) * 2020-12-08 2021-04-30 广东工业大学 一种基于双曲线阈值去噪的峰值提取方法
CN214097045U (zh) * 2020-12-17 2021-08-31 广州机械科学研究院有限公司 一种黏度传感器的标定系统
CN113533135A (zh) * 2021-07-15 2021-10-22 湖南工业大学 一种基于非牛顿流体流变特性测试系统的流量控制方法
CN114201900A (zh) * 2021-12-10 2022-03-18 西安石油大学 一种表征低渗储层非达西渗流的方法
WO2022099039A1 (en) * 2020-11-06 2022-05-12 Aramco Services Company Device and method to obtain the viscosity of polymer fluids used for enhanced oil recovery
CN114528780A (zh) * 2022-02-18 2022-05-24 青岛石大恒瑞石油新技术开发有限公司 一种考虑稠油非牛顿特性的数值模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021141957A1 (en) * 2020-01-06 2021-07-15 Saudi Arabian Oil Company Determining the rheological properties of a fluid through a non-linear response

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105004633A (zh) * 2015-06-24 2015-10-28 广东工业大学 基于纳米线垂直阵列的流体粘度检测装置及检测方法
CN106290065A (zh) * 2016-07-28 2017-01-04 北京建筑大学 一种确定非牛顿流体在反应器内搅拌死区的方法及装置
CN107271332A (zh) * 2017-07-04 2017-10-20 西安交通大学 一种基于面内谐振的mems流体黏度传感器芯片及其制备方法
CN108871200A (zh) * 2018-06-01 2018-11-23 广东工业大学 一种基于表面纳米轴向光子结构回音壁微腔的探针式微位移传感系统
CN110631959A (zh) * 2019-09-30 2019-12-31 南方科技大学 测量非牛顿流体表观粘度的端面效应误差消除方法及装置
CN111537396A (zh) * 2020-04-24 2020-08-14 西安交通大学 一种多操作模式的压电黏度传感器芯片及其工作方法和制备方法
WO2022099039A1 (en) * 2020-11-06 2022-05-12 Aramco Services Company Device and method to obtain the viscosity of polymer fluids used for enhanced oil recovery
CN112729152A (zh) * 2020-12-08 2021-04-30 广东工业大学 一种基于双曲线阈值去噪的峰值提取方法
CN214097045U (zh) * 2020-12-17 2021-08-31 广州机械科学研究院有限公司 一种黏度传感器的标定系统
CN113533135A (zh) * 2021-07-15 2021-10-22 湖南工业大学 一种基于非牛顿流体流变特性测试系统的流量控制方法
CN114201900A (zh) * 2021-12-10 2022-03-18 西安石油大学 一种表征低渗储层非达西渗流的方法
CN114528780A (zh) * 2022-02-18 2022-05-24 青岛石大恒瑞石油新技术开发有限公司 一种考虑稠油非牛顿特性的数值模拟方法

Also Published As

Publication number Publication date
CN115479867A (zh) 2022-12-16

Similar Documents

Publication Publication Date Title
Agoston et al. Evaluation of a vibrating micromachined cantilever sensor for measuring the viscosity of complex organic liquids
Cheng Yield stress: a time-dependent property and how to measure it
Giacomin et al. Large-amplitude oscillatory shear
US9494503B2 (en) Inline rheology/viscosity, density, and flow rate measurement
Barnes et al. Rotating vane rheometry—a review
Benmouffok-Benbelkacem et al. Non-linear viscoelasticity and temporal behavior of typical yield stress fluids: Carbopol, xanthan and ketchup
Jones et al. On the extensional viscosity of mobile polymer solutions
Jakoby et al. Physical sensors for water-in-oil emulsions
US11493419B2 (en) Measurement of properties of flowing vibrated yield stress fluids
CN104502235B (zh) 自闭式油液粘度在线检测传感器及其方法
US20160097743A1 (en) Method for measuring the properties of liquid based on a quartz crystal microbalance sensor
CN115479867B (zh) 一种非牛顿流体快速判别方法
Dogan et al. In‐line measurement of rheological parameters and modeling of apparent wall slip in diced tomato suspensions using ultrasonics
Konigsberg et al. Online process rheometry using oscillatory squeeze flow
Laun et al. New laboratory tests to measure rheological properties of paper coatings in transient and steady-state flows
US6711940B2 (en) Method and apparatus for measuring the elasticity of fluids
Macosko et al. Fluid mechanics measurements in non-Newtonian fluids
CN205091283U (zh) 一种超声波石油凝点测量装置
Tao et al. The importance of vane configuration on yield stress measurements of cement slurry
Savins et al. The Differentiation Method in Rheology: I. Poiseuille-Type Flow
CN112943224B (zh) 稠油井动液面的计算方法
CN102435536A (zh) 工程泥浆流变性能测试方法
CN108776083B (zh) 一种磁流变液再分散性定量测试方法
Buiochi et al. Viscosity measuring cell using ultrasonic wave mode conversion
CN207280900U (zh) 一种非接触式电容耦合触发开关

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