大庆油田已进入三次采油阶段，注聚合物驱油技术已经成为油田提高原油采收率的重要手段[1-2]，为使聚合物溶液进入预先设定油层并能得到一个较为均匀的聚合物驱前缘，需要准确确定从注聚井中进入各油层聚合物的注入量，因此，注聚井中流量测量是其中的重要测试内容[3-6]。大庆油田外流式四电极电磁流量计在注聚井中流量测井实践表明，电磁流量计在聚合物溶液及清水中测量响应有较大差别。由于仪器倾斜及流场中存在悬浮颗粒等各种因素，都会给电磁流量计响应带来影响，这些测井环境因素对电磁流量计响应影响须要从数值模拟角度给予理论解释与分析，从而为注聚井中电磁流量计流量测井提供理论分析基础。在分析四电极电磁流量计流场及磁场分布的基础上，笔者重点考察了流速剖面分布对电磁流量计响应输出特性影响，取得了对电磁流量计实验结果较好的分析效果。
1 四电极电磁流量计理论模型
根据Shercliff电磁流量计测量理论[3]，有
式中：U为电磁流量计感应电势输出，w为权重函数，v为局部流速分布，B为磁感应强度分布，S为流道内流体占据的截面积。
权重函数w含义为，磁场分布不同时，流体流过磁场时在管道截面上流体微元切割磁力线时产生的感应电势对总感应电势贡献大小不同。由式（1）可知，当权重函数w及流速剖面v为非轴对称分布时，流量测量会产生较大误差。
首先采用ANSYS有限元分析软件对电磁流量计测量区域内磁场分布特点进行考察。
在二维平面场（X-Y平面）中，矢量磁势A和电流密度J相互平行且只有z方向分量，即：Ax=Ay=0，Az=A；Jx=Jy=0，Jz=J。模型中介质为线性介质，磁导率μ为常数。由麦克斯韦方程导出的分矢量泊松方程为
模型有两种边界条件：①Dirichlet条件（AZ约束），即磁通量平行于模型边界；②Neumann条件（自然边界条件），即磁通量垂直于模型边界。第2种条件为默认的边界条件。对于电磁流量计在管道中的模型，只须满足自然边界条件，故施加电流密度后可进行磁场计算。
设定水的相对磁导率μr=1，由于聚合物溶液磁导率小于水的磁导率[7]，在模拟计算时，假定聚合物溶液相对磁导率μr=0.5，如图1所示。从图1可以看出：聚合物溶液产生的磁通线稀疏，而水中磁通线密度则较大，其磁场强度也越大。由式（1）可知，磁场强度越大，其电磁流量计灵敏度越高。由于仪器结构尺寸非常对称，仪器位于管道中心，通电后4个线圈相当于交替放置的N极与S极，故产生的磁场也是对称分布的。流体从仪器与油管环形空间流过，切割磁力线产生感生电势，通过4个对称分布的电极即可进行流量测量。
2 电磁流量计在聚合物溶液及清水管流中动态实验
外流式四电极电磁流量计在聚合物溶液中的动态实验是在大庆油田测试技术服务分公司多相流流动环路中进行的，实验目的是为了认识电磁流量计在清水及聚合物溶液中响应特性。使用模拟井内径为125mm，井筒高度为13m，内衬有长为6m、内径为62mm的油管，实验时用扶正器使仪器居中。实验用聚合物溶液为相对分子质量为（3～8）×106的聚丙烯酰胺，配制不同质量分数、密度约为1.0g/cm3聚合物。在清水中标定的流量范围为0～200m3/d，在聚合物溶液中标定的流量范围为0～150m3/d。聚合物溶液黏度与聚合物质量分数实验关系可通过化验分析获得。
图2为井下四电极电磁流量计及测量响应曲线图。从实验结果可以看出，在清水中的测量响应线在相同流量时总是高于聚合物响应线，尤其是流量增大时差异更明显。此外，当聚合物溶液黏度变化时，每支仪器对黏度变化的响应值差别很小，表明电磁流量计受流体黏度的影响不显著。
3 四电极电磁流量计响应数值模拟
3.1 计算测量场域内流速分布的标准κ-ε模型
在标准κ-ε模型中[8]，κ和ε为两个基本未知量，与之相对应的输运方程为
式中：ui表示时均速度；Gk为由于平均速度梯度引起的湍动能k的产生项；Gb为由于浮力引起的湍动能k的产生项；Ym代表可压湍流中脉动扩张的贡献；μ为因分子黏性而引入的流体动力黏度；μt为湍动黏度（空间坐标函数，取决于流动状态）；C1ε，C2ε和C3ε为经验常数；σk和σε分别为与湍动能k和耗散率ε对应的Prandtl数；Sk和Sε为用户定义的源项。
在此模型中，根据Launder等推荐值，模型常数取值分别为：C1ε=1.44；C2ε=1.92；C3ε=1；σk=1.0；σε=1.3。
3.2 油管与仪器环形空间流速分布
仿真模型油管管径为62mm，油管长度设定为2000mm，仪器直径为35mm。为了在仪器使用过程中测量电极不受损坏，测量处专门设计成凹槽形状。所以测量处仪器直径为33.8mm，仪器长度为1200mm，仪器测量电极段长度为44.5mm。计算仪器与油管环形空间流场计算时单元网格剖分数为（10～60）×104。由于油管长度较长且与管径比例相差太大，无法完整显示网格，这里忽略变化很小的部分，仅取仪器顶部、仪器测量段及仪器尾部网格，合并后的整体网格剖分情况如图3所示。若设定聚合物质量分数变化范围为0.05%～0.2%，则计算出的测量电极处流速分布呈中心对称，随着黏度的增大，流速分布中心流速也随之增大。清水流速剖面分布比聚合物流速剖面分布平坦些，其中心流速比聚合物小，而环形空间两侧存在比聚合物流速大的区域。尽管聚合物中心流速比清水要大，但是清水在中心区域外的两边流速要比对应的聚合物流速要大。
3.3 四电极电磁流量计响应数值模拟
四电极电磁流量计权重函数可认为两电极权重函数的叠加。两电极权重函数的表达式为[3]
则四电极权重函数可以表示为
式中：a为归一化的管子半径；x与y分别为测量场域内平面直角坐标位置变量。
根据前面采用有限元分析方法计算的磁场强度分布，结合测量场域内流速剖面分布及权重函数分布，最后由式（1）就可以计算出四电极电磁流量计响应输出，如图4所示。
由图4可以看出：假设清水磁导率与聚合物不同，清水中电磁流量计响应大于聚合物中的响应；尽管聚合物黏度变化，但在相同聚合物流量时，电磁流量计响应输出结果差别不很明显，说明聚合物在轴对称流速分布时，黏度变化对流量测量影响不大。以上数值模拟结果也与图2所示的在清水与聚合物溶液中实验测量结果相吻合，较好地解释了在清水与聚合物溶液中电磁流量计响应差异的原因，并对聚合物黏度变化对流量测量影响不大给出了理论计算依据。
4 测井环境对电磁流量计响应影响
4.1 仪器在油管中倾斜时的流速分布
从图5可以看出：当仪器向左侧倾斜1°时，仪器顶部及测量段流场分布有很大变化，说明仪器倾斜时对流场分布影响很大；当仪器在油管中倾斜时，电磁流量计权重函数会发生较大的非对称分布，仪器倾斜带来的流速非对称分布会给仪器测量精度带来较大影响。实际测井时，应尽可能保持仪器的居中扶正效果。
4.2 存在悬浮物颗粒时流速分布
图6（a）为注入聚合物流量为100m3/d时流场有悬浮物颗粒时的流速分布，其中聚合物黏度为14mPa·s，油管直径为62mm，仪器直径为38mm，悬浮物颗粒直径为4mm。注入不同流量聚合物时测量电极处流速分布如图6（b）所示。可以看出：有悬浮物颗粒时流速剖面发生很大畸变，尤其是高流量时畸变更为严重，这种情况下电磁流量计测量聚合物流量会有较大误差。
5 结论
（1）从磁场与流速剖面分布角度，对清水及聚合物溶液中的电磁流量计测量响应存在的差异进行了数值模拟计算，验证了清水中仪器常数高于聚合物中的仪器常数。
（2）根据理论上计算的电磁流量计测量区域内流速剖面在不同黏度聚合物中的分布规律可知，电磁流量计在轴对称流动的聚合物溶液中仪器常数变化不显著。
（3）理论上计算了仪器倾斜及悬浮颗粒存在时电磁流量计测量区域内流速剖面分布，得到流速剖面非轴对称分布，进而影响电磁流量计测量响应特性。
参考文献
[1]HanDK，YangCZ，ZhangZQ，etal。Recentdevelopmentofen2hancedoilrecoveryinChina[J]。JournalofPetroleumScienceandEngineering，1999，22（1/3）：181-188。
[2]王德民，程杰成，吴军政，等。聚合物驱油技术在大庆油田的应用[J]。石油学报，2005，26（1）：74-78。 WangDemin，ChengJiecheng，WuJunzheng，etal。ApplicationofpolymerfloodingtechnologyinDaqingOilfield[J]。ActaPetroleiSiNIca，2005，26（1）：74-78。
[3]ShercliffJA。Thetheoryofelectromagneticflowmeasurement[M]。Cambridge：CambridgeUniversityPress，1962：282-311。
[4]BevirMK。Thetheoryofinducedvoltageelectromagneticflowme2ters[J]。JournalofFluidMechanics，1970，43（3）：577-590。
[5]王国强，吕波，张小章。非对称流动对电磁流量计输出的影响[J]。计量技术，2003，47（10）：3-5。
WangGuoqiang，LüBo，ZhangXiaozhang。Effectoftheasymmet2ricalflowonanelectromagneticflowmeter[J]。MeasurementTechnique，2003，47（10）：3-5。
[6]徐立军，王亚，乔旭彤，等。多对电极电磁流量计传感器电极阵列设计[J]。仪器仪表学报，2003，24（4）：335-339。
XuLijun，WangYa，QiaoXutong，etal。DesignofelectrodearrayofmultiOelectrodeelectromagneticflowmeter[J]。ChineseJournalofScientificInstrument，2003，24（4）：335-339。
[7]孙作达，李晓军，余萍。电磁感应式磁化率定性分析装置原理与结构[J]。现代仪器，2006，12（5）：45-48。SunZuoda，LiXiaojun，YuPing。rincipleandconstructionofaqualitativeanalysisapparatusbasedonelectromagneticinductionsusceptibility[J]。ModernInstruments，2006，12（5）：45O48。
[8]LaunderBE，SpaldingDB。Lecturesinmathematicalmodelsofturbulence[M]。London：AcademicPress，1972：90O110。
