发明名称 一种火电机组给水流量测定方法
摘要 本发明公开了一种火电机组给水流量测定方法,取用机组设计工况给水流量、机组功率、给水管道内径、待测工况给水压力、给水温度、给水动力粘度;测点位置、毕托管测得的测点处给水动压。计算得到各测点处给水流速,计算得到给水流量估算值、给水平均流速、给水管内流动雷诺数,确定给水沿管道截面流速分布。利用给水沿管道截面流速分布,计算测点处给水流速与平均给水流速的关系,确定各测点处流速权数。利用各测点处测得的给水流速,结合各测点处流速权数,计算平均给水流速和给水流量,并与给水流量初始值比较,若偏差不满足精度要求,则进行迭代计算,直至偏差满足精度要求,确定给水流量。
申请公布号 CN105091957A 申请公布日期 2015.11.25
申请号 CN201510279363.8 申请日期 2015.05.27
申请人 东南大学 发明人 王培红;王德鹏;赵刚;苏志刚;郝勇生
分类号 G01F1/36(2006.01)I 主分类号 G01F1/36(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 一种火电机组给水流量测定方法,其特征在于,包括如下步骤:(1),测点布置与测点流速测算:在给水管道任一截面内沿半径方向选取距离截面圆心处长度r=0、r=0.4R、r=0.8R、r=0.95R的四个位置为测点,利用毕托管测量给水的动压Δp<sub>i</sub>,则各测点处给水流速测量值v<sub>mi</sub>为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>v</mi><mi>mi</mi></msub><mo>=</mo><msqrt><mfrac><mrow><mn>2</mn><mi>&Delta;</mi><msub><mi>p</mi><mi>i</mi></msub></mrow><mi>&rho;</mi></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000011.GIF" wi="915" he="141" /></maths>式中:v<sub>mi</sub>——i取1~4,分别对应r=0、r=0.4R、r=0.8R、r=0.95R测点处测得的给水流速测量值;Δp<sub>i</sub>——i取1~4,分别对应r=0、r=0.4R、r=0.8R、r=0.95R测点处测得的给水动压测量值;ρ——给水密度;R为给水管道半径;(2),预估给水流量,得到给水流量初始值:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>w</mi><mn>1</mn></mrow></msub><mo>&ap;</mo><msub><mi>D</mi><mrow><mi>w</mi><mn>0</mn></mrow></msub><mfrac><msub><mi>P</mi><mi>el</mi></msub><msub><mi>P</mi><mrow><mi>el</mi><mn>0</mn></mrow></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000012.GIF" wi="926" he="110" /></maths>式中:D<sub>w1</sub>——给水流量初始值;P<sub>el</sub>——待测工况机组功率;D<sub>w0</sub>——取用的设计工况给水流量;P<sub>el0</sub>——取用的设计工况机组功率;(3),计算各测点流速权数:31),利用所述给水流量初始值D<sub>w1</sub>、给水管道内径和给水密度,预估给水平均流速初始值v<sub>a</sub>’:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><msub><mi>v</mi><mi>a</mi></msub><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><mn>4</mn><msub><mi>D</mi><mrow><mi>w</mi><mn>1</mn></mrow></msub></mrow><mrow><mi>&rho;&pi;</mi><msup><mi>d</mi><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000013.GIF" wi="926" he="108" /></maths>式中:v<sub>a</sub>’——给水平均流速初始值;d——给水管道内径;32),根据给水平均流速初始值,计算给水管道内流动的雷诺数Re:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>Re</mi><mo>=</mo><mfrac><mrow><mi>&rho;</mi><msup><msub><mi>v</mi><mi>a</mi></msub><mo>&prime;</mo></msup><mi>d</mi></mrow><mi>&eta;</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000014.GIF" wi="920" he="112" /></maths>式中:η——给水动力粘度;33),在紊流情况下,即Re&gt;2000时,计算给水管道截面径向的速度分布v为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>v</mi><mo>=</mo><msub><mi>v</mi><mi>max</mi></msub><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>r</mi><mi>R</mi></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>n</mi></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000015.GIF" wi="918" he="132" /></maths>式中v<sub>max</sub>为给水管道轴心处给水流速,紊流速度分布常数n计算如式(7)所示:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>n</mi><mo>=</mo><mn>10</mn><mo>-</mo><msup><mrow><mn>3</mn><mi>e</mi></mrow><mrow><mo>-</mo><mn>1</mn><msup><mn>0</mn><mrow><mo>-</mo><mn>6</mn></mrow></msup><mi>Re</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000021.GIF" wi="923" he="79" /></maths>34),根据所述给水管道截面径向的速度分布v以及管道截面面积,根据式(8)建立沿给水管道截面平均速度v<sub>a</sub>的模型:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>v</mi><mi>a</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><mi>R</mi></msubsup><mn>2</mn><mi>&pi;rvdr</mi></mrow><mrow><mi>&pi;</mi><msup><mi>R</mi><mn>2</mn></msup></mrow></mfrac><mo>=</mo><mn>2</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><msub><mi>v</mi><mi>max</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000022.GIF" wi="931" he="120" /></maths>35),在给水管道截面径向的速度分布v的情况下,估计各测点对应的平均流速:将各测点距离截面圆心处长度r代入式(6),得到各测点处预测流速v<sub>1</sub>、v<sub>2</sub>、v<sub>3</sub>、v<sub>4</sub>:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>v</mi><mn>1</mn></msub><mo>=</mo><msub><mi>v</mi><mi>max</mi></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mn>2</mn></msub><mo>=</mo><msup><mn>0.6</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup><msub><mi>v</mi><mi>max</mi></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mn>3</mn></msub><mo>=</mo><msup><mn>0.2</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup><msub><mi>v</mi><mi>max</mi></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mn>4</mn></msub><mo>=</mo><msup><mn>0.05</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup><msub><mi>v</mi><mi>max</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000023.GIF" wi="917" he="323" /></maths>联立(8)和式(9),消去v<sub>max</sub>,得到各测点对应的平均流速预测值v<sub>a1</sub>、v<sub>a2</sub>、v<sub>a3</sub>、v<sub>a4</sub>:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>v</mi><mrow><mi>a</mi><mn>1</mn></mrow></msub><mo>=</mo><mn>2</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><msub><mi>v</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>a</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.6</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac><msub><mi>v</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>a</mi><mn>3</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.2</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac><msub><mi>v</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mrow><mi>a</mi><mn>4</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.05</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac><msub><mi>v</mi><mn>4</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000024.GIF" wi="934" he="526" /></maths>36),根据各测点对应的平均流速预测值v<sub>a1</sub>、v<sub>a2</sub>、v<sub>a3</sub>、v<sub>a4</sub>,计算沿给水管道截面平均速度值,如式(11)所示:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>v</mi><mi>a</mi></msub><mo>=</mo><mfrac><mrow><mmultiscripts><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></msubsup></mmultiscripts><msub><mi>v</mi><mi>ai</mi></msub></mrow><mn>4</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000025.GIF" wi="935" he="125" /></maths>将式(10)代入(11),得<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>v</mi><mi>a</mi></msub><mo>=</mo><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow><msub><mi>v</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.6</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac><msub><mi>v</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.2</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac><mrow><msub><mi>v</mi><mn>3</mn></msub><mo>+</mo><mfrac><mrow><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.05</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac></mrow><msub><mi>v</mi><mn>4</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000026.GIF" wi="1696" he="144" /></maths>根据式(13),计算各测点流速权数K<sub>1</sub>、K<sub>2</sub>、K<sub>3</sub>、K<sub>4</sub>:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>K</mi><mn>1</mn></msub><mo>=</mo><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.6</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.2</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>4</mn></msub><mo>=</mo><mfrac><mrow><mn>0.5</mn><mrow><mo>(</mo><mfrac><mi>n</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mfrac><mi>n</mi><mrow><mn>2</mn><mi>n</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>)</mo></mrow></mrow><msup><mn>0.05</mn><mfrac><mn>1</mn><mi>n</mi></mfrac></msup></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000031.GIF" wi="934" he="524" /></maths>(4),迭代逻辑:由所述各测点的给水流速测量值v<sub>mi</sub>和各测点的流速权数K<sub>i</sub>,计算基于测点流速测量值的给水管道截面的平均流速:v<sub>ma</sub>=∑K<sub>i</sub>v<sub>mi</sub>   (14)则所测工况的给水流量为:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>D</mi><mi>w</mi></msub><mo>=</mo><mfrac><mrow><mi>&rho;</mi><msub><mi>v</mi><mi>ma</mi></msub><mi>&pi;</mi><msup><mi>d</mi><mn>2</mn></msup></mrow><mn>4</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000725510350000032.GIF" wi="936" he="115" /></maths>将计算得到的给水流量D<sub>w</sub>与初始计算值D<sub>w1</sub>进行误差校验,若相对误差大于预设值,则以计算流量为初始流量,重新计算,直至流量相对误差小于等于预设值,此时得到的给水流量即为待测工况的给水流量。
地址 211189 江苏省南京市江宁区东南大学路2号