CN108959683B - 一种基于cfd的数字化锅炉构建方法 - Google Patents
一种基于cfd的数字化锅炉构建方法 Download PDFInfo
- Publication number
- CN108959683B CN108959683B CN201810312543.5A CN201810312543A CN108959683B CN 108959683 B CN108959683 B CN 108959683B CN 201810312543 A CN201810312543 A CN 201810312543A CN 108959683 B CN108959683 B CN 108959683B
- Authority
- CN
- China
- Prior art keywords
- value
- initial positioning
- tdoa
- leakage point
- temperature field
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
本发明涉及一种基于CFD的数字化锅炉构建方法,包括:1)基于传感器阵列接收压力管道泄漏点信号,获得各TDOA值;2)将炉内的温度场假定为均匀温度场,并给出初始定位值;3)采用投影方法调整初始定位值坐标,将初始定位值投影至离泄漏点最近的墙壁上;4)以初始定位值为泄漏点,基于炉膛温度场分布数据库,计算由初始定位值所获得的各TDOA值;5)对比测量与计算所得TDOA值,判断初始定位值与实际泄漏点的大致方向;6)设定移动步长,通过初始定位值的位移获得新的泄漏位置,进行TDOA计算;7)重复步骤4)及步骤5),直至TDOA值小于预设误差。本发明克服了未知非均匀温度场中声速不固定而导致被动定位算法时延值无法反推距离差的问题。
Description
技术领域
本发明属于锅炉泄漏定位技术领域,尤其涉及一种基于CFD的数字化锅炉构 建方法。
背景技术
由于燃煤锅炉内环境恶劣,炉内压力管道泄漏事故时有发生,若能够及时 定位出泄漏位置,则能够缩短检修时间,提高经济效益。
虽然已有锅炉燃烧模型及以此为基础构建数据库的研究,但主要用于燃烧 及污染物排放方面;同时,泄漏定位的研究主要集中于冷态方面,声线弯曲的 研究主要集中于已知温度场的重建方面。以锅炉燃烧模型构建温度场数据库用 于炉膛热态非均匀温度场的定位方法尚未有学者进行研究。
发明内容
本发明的目的是提供一种基于CFD的数字化锅炉构建方法,将锅炉温度场 数据库用于锅炉泄漏定位,通过初值迭代寻优模型克服未知非均匀温度场中声 速不固定而导致被动定位算法时延值无法反推距离差的问题。
本发明提供了一种基于CFD的数字化锅炉构建方法,
步骤1,基于传感器阵列接收压力管道泄漏点信号,获得各TDOA值;
步骤2,将炉内的温度场假定为均匀温度场,并给出初始定位值;
步骤3,采用投影方法调整初始定位值坐标,将初始定位值投影至离泄漏点 最近的墙壁上;
步骤4,以初始定位值为泄漏点,基于炉膛温度场分布数据库,计算由初始 定位值所获得的各TDOA值;
步骤5,对比测量与计算所得TDOA值,判断初始定位值与实际泄漏点的大 致方向;
步骤6,设定移动步长,通过初始定位值的位移获得新的泄漏位置,进行TDOA计算;
步骤7,重复步骤4及步骤5,直至TDOA值小于预设误差。
进一步地,步骤5包括:
以传感器的连线与泄漏点所在墙壁平行为原则选取用于判断移动方向的 TDOA值。
借由上述方案,通过基于CFD的数字化锅炉构建方法,将锅炉温度场数据 库用于锅炉泄漏定位,通过初值迭代寻优模型克服了未知非均匀温度场中声速 不固定而导致被动定位算法时延值无法反推距离差的问题。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术 手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附 图详细说明如后。
附图说明
图1是泄漏点定位模型;
图2是本发明基于CFD的数字化锅炉构建方法泄漏点位置更新示意图;
图3是本发明基于CFD的数字化锅炉构建方法传感器阵列坐标示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以 下实施例用于说明本发明,但不用来限制本发明的范围。
本实施例提供了一种基于CFD的数字化锅炉构建方法,包括如下步骤:
步骤1,传感器阵列接收压力管道泄漏点信号,获得各TDOA(时延)值;
步骤2,将炉内的温度场假定为均匀温度场给出初始定位值;
步骤3,由于炉膛内泄漏点均位于墙壁上,因此需要调整初始定位值坐标, 采用投影的方法将初始定位值投影至离其最近的墙壁上;
步骤4,以初始定位值为泄漏点,调用炉膛温度场分布数据库,计算由初始 值所获得的各TDOA值;
步骤5,对比测量与计算所得TDOA值,判断初始定位值与实际泄漏点的大 致方向;
步骤6,设定移动步长,通过初始定位值的位移获得新的泄漏位置,进行 TDOA计算;
步骤7,重复步骤4及步骤5,直至TDOA值小于预设误差。
通过该基于CFD的数字化锅炉构建方法,将锅炉温度场数据库用于锅炉泄 漏定位,通过初值迭代寻优模型克服了未知非均匀温度场中声速不固定而导致 被动定位算法时延值无法反推距离差的问题。
下面对本发明作进一步详细说明。
1、冷态情况下定位方法原理(球型插值法,简称SI)。
参图1所示,1为传感器,2为泄漏点,3为压力管道。
对传感器位置进行标定,另传感器数量为n,对应坐标为(xi,yi,zi),其中 i=1,2,...,n,假设压力管道泄漏点坐标为rleak=(x,y,z)。
建立方程组首先需要选定参考点传感器,将其坐标记为r1=(x1,y1,z1),那么rleak到r1的距离可表示为Rleak=||rleak-r1||,任意传感器到参考点传感器的距离可表示为 Rj=||rj-r1||,其中j=2,...,n。令泄漏点到rj和r1的距离差为dj,其表达式如下:
其中,dj可以通过测量值计算获得,例如泄漏点到rj和r1的时延值(TDOA)t1j由传感器测量后进行互相关计算而得到,将TDOA值乘以声速既得到dj。
式(1)经推导后可得出:
式(2)化简后可得出:
ε=δ-2Rleakd-2rTrleak (4)
其中,
如令Rleak为已知数,则式(4)对于rleak是线性的,为满足其均方误差最小,需 符合下式的情况:
其中,W=(rTr)-1rT。
a=4-4dTWTWd
其中,b=4dTWTWδ。
c=-δWTWδ
式(6)有且仅有一个正根,表示为:
将式(7)代入式(5)可得:
2、温度与声速关系(用于计算热态温度场下时延值)。
声速随温度的变化而变化,在声学基本理论中有详细的推导,在此不做过 多描述,其变化规律如下:
其中,γ表示比热比,无量纲;k=γ/M,M表示气体摩尔质量,kg/mol;P 表示压强,Pa;ρ表示气体密度,kg/m3;R表示气体适普常数,8.31J/mol·K; T表示温度函数,K。
3、热态温度场压力管道泄漏定位模型。
在均匀温度场情况下,声速恒定,因此声波飞渡时间一旦确定,则能够获 知声源距传感器的距离,因此在被动定位算法中能够通过时延值反推出距离差 dj;而非均匀温度中,由于声速不固定,无法根据时延值反推出距离差。
为解决此问题,本发明提出了初值迭代寻优的模型。为显示方便,以二维 定位模型为例进行阐述,如图2所示。由于t24’<t24,即实际泄漏点比初始值更 靠近2号传感器(如图3所示),因此通过判断初始值与测量值的TDOA差异即 可判断出泄漏点位置更新的移动方向。需要注意的是,在选择判断移动方向的 TDOA值时,以传感器的连线与泄漏点所在墙壁平行为原则进行选取。
参图3所示,数字1-8表示1-8号传感器,以泄漏点坐标(5,9,1)为例对 初值寻优模型进行分析,首先根据数据库调取炉内温度场分布情况,计算温度 场平均值为1136.7K为例,声速为658.9m/s,采用SI算法得到初始泄漏位置为 (3.40,4.84,-1.20),由于炉膛内传感器阵列分段定位,因此,泄漏点在z方 向的值为0到5范围内,同时泄漏点更靠近y=9截面的墙壁,因此投影后,泄 漏点坐标变更为(3.4,9,0)。在单峰温度场情况下,投影后的泄漏点坐标与实 际坐标获得的时延值列于表1。
表1 TDOA
通过对比t31可知,实际泄漏点到3号与1号传感器的距离差要大于投影泄 漏点到两传感器的距离差,即实际泄漏点比投影点在x轴上的值更大,因此投 影点在更新坐标时应向x轴正方向移动。由于以传感器1为参考点判断投影点 移动方向效果不明显,因此选取t62进行对比,实际泄漏点与投影点的值分别为 0.002313ms、0.004517ms,实际泄漏点到6号与2号传感器的距离差要小于投 影泄漏点到两传感器的距离差,即实际泄漏点比投影点在z轴上的值更大,因 此投影点在更新坐标时应向z轴正方向移动。不断更新泄漏点位置,直至TDOA 值小于预设误差。
以上所述仅是本发明的优选实施方式,并不用于限制本发明,应当指出, 对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还 可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。
Claims (2)
1.一种基于CFD的数字化锅炉构建方法,其特征在于,包括如下步骤:
步骤1,基于传感器阵列接收压力管道泄漏点信号,获得各TDOA值;
步骤2,将炉内的温度场假定为均匀温度场,将步骤1获得的TDOA值乘以声速得到传感器到泄露点的距离差,根据得到的距离差,基于冷态情况下定位方法,给出初始定位值;
步骤3,采用投影方法调整初始定位值坐标,将初始定位值投影至离泄漏点最近的墙壁上;
步骤4,以初始定位值为泄漏点,根据初始定位值对应的泄露点到传感器的距离差,以及基于炉膛温度场分布数据库、温度与声速关系得出的声速,计算由初始定位值所获得的各TDOA值;
步骤5,对比测量与计算所得TDOA值,判断初始定位值与实际泄漏点的大致方向;
步骤6,设定移动步长,通过初始定位值的位移获得新的泄漏位置,进行TDOA计算;
步骤7,重复步骤4及步骤5,直至TDOA值小于预设误差。
2.根据权利要求1所述的一种基于CFD的数字化锅炉构建方法,其特征在于,所述步骤5包括:
以传感器的连线与泄漏点所在墙壁平行为原则选取用于判断移动方向的TDOA值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810312543.5A CN108959683B (zh) | 2018-04-09 | 2018-04-09 | 一种基于cfd的数字化锅炉构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810312543.5A CN108959683B (zh) | 2018-04-09 | 2018-04-09 | 一种基于cfd的数字化锅炉构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108959683A CN108959683A (zh) | 2018-12-07 |
CN108959683B true CN108959683B (zh) | 2021-03-23 |
Family
ID=64499157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810312543.5A Active CN108959683B (zh) | 2018-04-09 | 2018-04-09 | 一种基于cfd的数字化锅炉构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108959683B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112763399B (zh) * | 2020-12-28 | 2021-09-21 | 西安特种设备检验检测院 | 一种锅炉低温受热面烟气腐蚀风险区域的检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813545A (zh) * | 2010-04-07 | 2010-08-25 | 华北电力大学 | 立体四元阵列电站锅炉承压管泄漏的精确定位方法 |
CN102410525A (zh) * | 2011-09-01 | 2012-04-11 | 华北电力大学 | 一种炉管泄漏声学监测定位方法及装置 |
CN102829925A (zh) * | 2012-08-09 | 2012-12-19 | 华北电力大学 | 平面八元阵列电站锅炉承压管泄漏定位方法 |
CN107063584A (zh) * | 2017-04-22 | 2017-08-18 | 中南大学 | 一种锅炉管泄漏判别与定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110203535A1 (en) * | 2010-02-19 | 2011-08-25 | Nrg Energy, Inc. | Method and System for Sootblower Flow Analyzer |
-
2018
- 2018-04-09 CN CN201810312543.5A patent/CN108959683B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813545A (zh) * | 2010-04-07 | 2010-08-25 | 华北电力大学 | 立体四元阵列电站锅炉承压管泄漏的精确定位方法 |
CN102410525A (zh) * | 2011-09-01 | 2012-04-11 | 华北电力大学 | 一种炉管泄漏声学监测定位方法及装置 |
CN102829925A (zh) * | 2012-08-09 | 2012-12-19 | 华北电力大学 | 平面八元阵列电站锅炉承压管泄漏定位方法 |
CN107063584A (zh) * | 2017-04-22 | 2017-08-18 | 中南大学 | 一种锅炉管泄漏判别与定位方法 |
Non-Patent Citations (1)
Title |
---|
基于多时延识别的电站锅炉多源泄漏被动定位方法;安连锁 等;《动力工程学报》;20151031;第35卷(第10期);第798-804页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108959683A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pipitone et al. | Determination of TDC in internal combustion engines by a newly developed thermodynamic approach | |
Tunestål | TDC offset estimation from motored cylinder pressure data based on heat release shaping | |
CN108763683B (zh) | 一种三角函数框架下新weno格式构造方法 | |
CN104390680B (zh) | 时差式超声波热量表流量误差修正方法 | |
CN101804581A (zh) | 一种机床热变形自动补偿的实现方法 | |
CN112100734B (zh) | 一种真空羽流对任意构型设备的影响确定方法 | |
CN107464591B (zh) | 一种安全壳泄漏在线监测系统传感器组合优化方法 | |
CN108959683B (zh) | 一种基于cfd的数字化锅炉构建方法 | |
CN115292986B (zh) | 一种液化天然气空温式气化器的性能评价方法 | |
JP4321777B2 (ja) | 熱解析方法およびそのプログラム | |
CN110068399A (zh) | 基于径向基函数和正则化的温度场重建算法 | |
CN111832108A (zh) | 基于倾角监测数据的结构位移响应估计方法和装置 | |
CN106295063A (zh) | 建筑物节点等效面传热系数计算方法 | |
Wu et al. | A method for reconstruction of boiler combustion temperature field based on acoustic tomography | |
CN105627916A (zh) | 一种建立跟踪仪地理坐标系及六自由度测量的方法 | |
CN112668224B (zh) | 一种空调房间温度传感器最优安装位置的确定方法及系统 | |
CN104881535A (zh) | 改进的火电厂锅炉温度场重建测温算法 | |
JP5114251B2 (ja) | 真空処理装置 | |
CN116467957A (zh) | 一种基于piv技术的超声速复杂流场的压力场重构方法 | |
Kong et al. | Research on Temperature Field Reconstruction Based on RBF Approximation with Polynomial Reproduction Considering the Refraction Effect of Sound Wave Paths. | |
CN106840516B (zh) | 一种基于多项式拟合的压力计温漂标定方法 | |
CN106052880B (zh) | 基于等效黑体截面的辐射测温距离确定方法 | |
Kyle et al. | Compensating for the effects of refraction in photogrammetric metrology | |
Wang et al. | A thermodynamics model for the compression and expansion process during the engine’s motoring and a new method for the determination of TDC with simulation technique | |
CN112033576A (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 |