Modeling and statistical analysis of distribution parameters of random cable bundles based on image recognition technology
-
摘要: 提出了一种对实际弯曲随机捆扎线束的建模方法，该方法首先基于图像识别技术，使用实际线束在侧视和俯视方向的两幅照片来实现弯曲线束轴心三维坐标的重建，然后再基于随机转移路径方法实现弯曲线束的捆扎随机性。基于该建模方法，通过蒙特卡洛模拟来分析弯曲随机线束分布参数的统计特征，发现自电感、互电感和互电容均值沿线变化趋势与线束高度变化趋势一致，自电容均值则趋势相反；自电容、自电感和互电感的变异系数与线束高度存在负相关特征；捆扎随机性不会改变自电感和自电容均值，但是会降低互电容与互电感均值。Abstract: In this paper, a modeling method of actual bending random bundled wire harness is proposed. Firstly, based on image recognition technology, the three-dimensional coordinates of bending wire harness axis are reconstructed by using two photos of actual wire harness in side view and top view; then the random bundled wire harness is realized based on random transfer path method. Based on this modeling method, this paper analyzes the statistical characteristics of distribution parameters of bending random wire harness by Monte Carlo simulation, and finds that the variation trend of self inductance, mutual inductance and mutual capacitance along the line is consistent with the variation trend of wire harness height, while the trend of self capacitance is opposite; the coefficient of variation of self capacitance, self inductance and mutual inductance has negative correlation with wire harness height; the bundling randomness is not obvious It will change the mean value of self inductance and self capacitance, but reduce the mean value of mutual capacitance and mutual inductance.
-
Key words:
- wire-bundles /
- random /
- image recognition /
- distributed parameter
-
自福岛核电站事故发生以来，反应堆的安全性越来越受到公众的关注，成为影响核能发展的关键因素。采用非能动系统来提高反应堆的固有安全性已经在新的反应堆设计中得到了广泛的应用。其中，利用依靠回路内冷热段内冷却剂的密度差来产生驱动压头带走堆芯余热的自然循环成为主要的研究方向。目前，国际上对自然循环的研究以试验和数值模拟为主。试验研究方面，杨祖毛等[1]对单相自然循环进行了大量的试验，研究了稳态自然循环流量与热源功率、进出口温差的关系；黄彦平[2]通过试验研究了自然循环能力的影响因素等。但试验研究大多针对具体的回路及系统，不具有良好的移植性，且需要大量的资金投入。数值模拟方面，国内外通常使用一些系统分析程序，如RELAP5，RETRAN-02等大型商用软件，对于某些反应堆，也会采用独立开发的专用程序，如田文喜等[3]开发了CARR堆的自然循环能力计算程序，黄洪文等[4]开发了池式研究堆的自然循环能力分析程序，Kaminaga等[5]开发了JRR-3M自然循环能力程序。针对JRR-3M类型的研究堆，稳定自然循环形成前存在冷却剂流动方向反转的过程，有必要开展瞬态安全计算分析。本文利用RELAP5/Mod 3.4对JRR-3M进行了细致的堆芯建模，研究了不同事故条件下自然循环建立前的流量反转过程瞬态热工特性，并得到了瞬态时自然循环的载热极限，同时验证了计算模型和计算方法的可行性。
1. 计算模型
1.1 系统概述
RELAP程序是美国爱达荷国家工程实验室开发、美国核管理委员会批准的用于工程审评的大型瞬态热工水力计算程序，是一维瞬态、两流体六方程流体力学及点堆动力学模型[6]，现已广泛用于压水堆的事故模拟分析。JRR-3M是一座重水慢化、轻水冷却的池式研究堆，在正常运行工况下，冷却剂在两台主泵的驱动下从堆芯上部水池进入堆芯，向下流经堆芯带走裂变热后流入下部的下联箱，与换热器内二次侧的冷却水进行热交换冷却后重新流入主泵。在自然循环阶段，下联箱与水池连接的自然循环阀开启，池水在自然循环压头的驱动下进入下联箱，向上流过堆芯带走衰变热后返回水池，形成自然循环回路。应急冷却系统由两台由不间断电源供电的辅助泵组成，它们的额定流量为75.00 kg·s-1，最多能持续工作3 h。在发生场外电源丧失事故时，一台辅助泵投入运行，带走堆芯热量。
在一回路失去强迫流动后，存在一个由于冷却剂流动方向反转所导致的零流量时刻，此时由于没有冷却剂的流动带走衰变热，堆芯温度迅速上升，该时刻堆芯最有可能被烧毁。基于RELAP5/Mod 3.4程序，对JRR-3M进行建模，反应堆节点图如图 1所示。
其中，考虑到堆水池中存在明显的热分层现象以及内部循环，用两个管道(100，102)来模拟堆芯上部水池，两根管间用横向接管连接，最上方的控制体为空气层，时间控制体130表示空气的温度压力边界，取35 ℃和1.013×105 Pa。104和106是两种不同燃料组件的平均通道，108表示热通道，110和112表示旁路流道，堆芯流道均划分为15个控制体。114表示下联箱，116是环绕堆芯的水池，它们之间用自然循环阀连接。正常运行时从堆芯流出的冷却剂在热交换器118中与202中的冷却水进行冷却热交换。冷却剂由主泵流出后在分支部件126汇合后流回堆芯上部水池。
对于堆芯内功率的分布取该池式反应堆的物理设计结果，其功率衰变曲线采用中国行业标准EJ/T 745-2001[7]推荐公式计算给出, 即
P(τ0)=G(t)4∑i=1∫τ00Pi(t)Qif(t)dt (1) Q(t)=1.0+τ0.40φ(3.24×10−6+5.23×10−10t) (2) 式中：P(τ0)是τ0时刻的裂变总功率；G(t)是考虑裂变产物中子俘获的校正因子；Pi(t)是核素i在t时刻的裂变功率；Qi为核素i一次裂变可回收能量(MeV)；f(t)为可裂变核素的衰变热功率函数，它的值由标准附表给出(略)；φ为一个初始裂变原子的裂变数；i取1，2，3，4分别代表235U热中子裂变、239Pu热中子裂变、238U快中子裂变和241Pu热中子裂变对衰变热功率的贡献；t为停堆后时间(s)。
1.2 稳态参数
将系统调节至稳态工况是进行瞬态计算的前提条件。在自然循环的瞬态计算中，系统各个参数应在某个预计工况下稳态运行，随后关闭主泵打开自然循环阀进行瞬态计算，而这个稳态的工况就决定了事故的初始状态。通过调节回路内的物理参数，使得稳态工况下的热工参数与设计的热工参数一致且达到稳定。表 1为该池式研究堆的热工设计参数。
表 1 池式研究堆主要设计参数Table 1. Main design parameters of pool-type research reactortotal power/MW total mass flow rate/(kg·s-1) inlet temperature/℃ outlet temperature/℃ core inlet pressure/MPa 20.00 661.91 35.00 42.19 1.69 1.3 安全准则
按照该堆设计确定的安全标准为：
(1) 使用Sudo临界热流密度经验公式[8]，计算的最小偏离泡核沸腾比(MDNBR)应大于限值1.5。
(2) 燃料元件芯体最高温度低于250 ℃。
2. 计算结果与分析
2.1 无应急冷却的场外电源丧失事故分析
反应堆在满功率正常运行了100 s后发生外电源丧失事故，主泵由于失去电力供应立即失电惰转，二次侧热阱立即丧失，自然循环阀打开，1 s后控制棒落下反应堆进行停堆，此时堆芯所产生的热量主要依靠主泵惰转和自然循环带走。整个过程的事件序列列于表 2。
表 2 无应急冷却系统事件序列Table 2. Accident process without emergency cooling systemevent start time event start time main pump loss of off-site power supply 100 s control rods drop 101 s secondary side loss of hot trap 100 s finish main pump coasting 220 s open the natural circulation valve 100 s reverse coolant flow 221 s 图 2为事故分析的结果曲线。
从图 2(a)和图 2(b)中可以看出，在0~100 s内，反应堆满功率稳态运行，堆芯流量保持不变。在100 s时刻，由于主泵失电惰转流量迅速减少，在t=220 s时主泵惰转停止，堆芯流量减少至0。此时自然循环阀已经打开，在自然循环驱动压头的作用下，由于热通道中热流密度最大，故最先产生流量反转，随后堆芯流量反转，冷却剂由下而上地通过堆芯，自然循环迅速建立。
从图 2(c)看出，在稳态运行时，燃料芯块最热温度为108.27 ℃，包壳最高温度为99.37 ℃, 在101 s时控制棒开始下落后，由于堆芯功率的呈指数式降低，燃料板和冷却剂的温度也随之急剧下降。随后，由于冷却剂流量不足以完全带走燃料所产生的热量，热构件和冷却剂的温度逐渐上升。在220 s时冷却剂流量为0，燃料外表面与冷却剂的对流换热性能最差，燃料板的温度达到峰值，燃料芯块的峰值温度为114.30 ℃, 包壳峰值温度为114.07 ℃。此后，由于自然循环流量的增大，对流换热作用增强，燃料板和冷却剂的温度开始下降。在稳定的自然循环建立后，燃料板和冷却剂的温度由于堆芯功率的缓慢衰减而缓慢降低。
图 2(d)是整个过程中的MDNBR变化曲线。从图中可以看出，热管处的MDNBR在稳态运行时为2.66，在反应堆停堆后由于燃料板热流密度的急剧减少而迅速增加，随后随着流量的降低而降低，在流量反转时刻为最小值，最小值为0.31，远低于安全限值1.50。
2.2 有应急冷却的场外电源丧失事故分析
当该池式研究堆在满功率正常运行了100 s后发生外电源丧失事故，主泵失电惰转，当冷却剂流量低于75.00 kg·s-1时，一台辅助泵投入运行。当t=3 188 s时，衰变功率已降至200 kW，辅助泵停止运行，在辅助泵停止运行前5 s打开自然循环阀。事故过程中的事件时间序列如表 3所示。
表 3 有应急冷却事件序列Table 3. Accident process with emergency cooling systemevent start time event start time main pump loss of off-site power supply 100 s close auxiliary pump 3 188 s control rods drop 101 s auxiliary pump finish coast 3 198 s input auxiliary pump 110 s coolant flow reverse 3 199 s open the natural circulation valve 3 183 s 图 3为事故分析的结果曲线。
从图 3(a)和图 3(b)中可以看出，在100 s以前，反应堆满功率稳态运行，堆芯流量不变。100 s后主泵失电惰转，流量降低。当流量降低至75.00 kg·s-1时，辅助泵投入运行，堆芯流量保持不变。当t=3 183 s时，自然循环阀打开，从图中可以发现堆芯流量平台有一个下降。这时由于此时堆水池与反应堆一回路连接，具有了一定的流量，导致流过堆芯的流量减少。在t=3 188 s时关闭辅助泵，堆芯流量逐渐降低至0，随后由于自然循环的建立，流量发生反转，最后趋向稳定。
通过比较图 2(c)和图 3(c)可以发现，在投入应急冷却系统的事故工况下，由于堆芯拥有一定程度的强迫流量，这一部分的流量带走了部分衰变热，有效地降低了燃料板及冷却剂的温度，相较于未投入应急冷却系统的条件下，燃料芯块峰值温度降低了46.10 ℃，堆芯材料不会发生熔毁。同时可以发现冷却剂的温度和燃料板温度的变化趋势相似，但有一定的滞后。
从图 3(d)可以看出，有应急冷却时事故过程中的MDNBR始终大于安全限值，不会发生偏离核态沸腾，池式研究堆始终处于安全状态。
2.3 最大载热能力计算
为进一步探究池式研究堆自然循环瞬态过程中的最大载热能力，在堆芯功率衰减至不同功率水平时关闭辅助泵并打开自然循环阀，不同功率下自然循环过程中的最高燃料温度和MDNBR如图 4所示。
从图 4可以明显看出，此时反应堆的安全性主要受到来自MDNBR的限制，因此，池式研究堆最大瞬态自然循环能力约为590 kW。
3. 结论
本文利用系统性分析程序RELAP5/Mod 3.4对JRR-3M池式研究堆进行了建模，计算了在无应急冷却和有应急冷却情况下的丧失场外电源事故。RELAP5很好地模拟了堆芯流量的反转过程以及自然循环的建立过程。针对这两种事故工况，计算得到了堆芯各个流道的温度、压力、冷却剂流量等热工水力参数。计算结果表明：若未投入应急冷却系统，反应堆在失去场外电源后热管内MDNBR为0.31，小于1.5的安全限值，会发生偏离核态沸腾，存在烧毁堆芯的风险。应急冷却系统能有效地降低事故后燃料板及冷却剂的温度，提高了系统的安全性。在瞬态条件下，池式研究堆的最大自然循环能力为590 kW，大于设计值200 kW，表明该堆具有良好的固有安全性。
-
图 4 一组弯曲随机线束样本建模效果图
Figure 4. Modeling effect picture of a group of bending random harness samples
图 5 本文分布参数计算与商业软件结果对比图
Figure 5. The comparison between the results of distribution parameter calculation and commercial software in this paper
图 6 一组平直随机线束样本建模效果图
Figure 6. Modeling effect picture of a group of straight random harness samples
图 7 弯曲未随机线束分布参数的沿线变化
Figure 7. Variation of distribution parameters of bending non random wire harness along the line
图 8 平直随机线束分布参数沿线变化
Figure 8. Variation of distribution parameters of straight random wire harness along the line
图 9 弯曲随机线束分布参数沿线变化
Figure 9. Variation of distribution parameters of bending random wire harness along the line
-
[1] Ciccolella A, Canavero F G. Stochastic prediction of wire coupling interference[C]//Proceedings of International Symposium on Electromagnetic Compatibility. 1995: 51-56. [2] Goodchild M F. Fractals and the accuracy of geographical measures[J]. Journal of the International Association for Mathematical Geology, 1980, 12(2): 85-98. doi: 10.1007/BF01035241 [3] Salio S, Canavero F, Lefèbvre J, et al. Statistical description of signal propagation on random bundles of wires[C]//Presented at the 13th International Zurich Symposium and Technical Exhibition on Electromagnetic Compatibility. Zurich, Switzerland: EMC, 1999. [4] Salio S, Canavero F, Lecointe D, et al. Crosstalk prediction on wire bundles by Kriging approach[C]//Proceedings of the IEEE International Symposium on Electromagnetic Compatibility. Washington: IEEE, 2000. [5] Sun Shishuang, Liu Geping, Drewniak J L, et al. Hand-assembled cable bundle modeling for crosstalk and common-mode radiation prediction[J]. IEEE Transactions on Electromagnetic Compatibility, 2007, 49(3): 708-718. doi: 10.1109/TEMC.2007.897142 [6] Bellan D, Pignari S A. Efficient estimation of crosstalk statistics in random wire bundles with lacing cords[J]. IEEE Transactions on Electromagnetic Compatibility, 2011, 53(1): 209-218. doi: 10.1109/TEMC.2010.2093147 [7] Spadacini G, Grassi F, Pignari S A. Statistical properties of low frequency voltages induced by a plane-wave field across the terminal loads of a random wire-bundle[C]//Proceedings of 2015 IEEE International Symposium on Electromagnetic Compatibility. Dresden, Germany: IEEE, 2015: 824-829. [8] 马喜来. 汽车电磁兼容性预估计的研究[D]. 长春: 吉林大学, 2008: 5-12.Ma Xilai. Research on the pre-assessments of electromagnetic compatibility in automobile[D]. Changchun: Jilin University, 2008: 5-12 [9] 王瑞宝, 周诚祥. 船舶线束寄生参数和串扰动态变化估计[J]. 中国修船, 2013, 26(2):34-36, 40. (Wang Ruibao, Zhou Chengxiang. Estimation of parasitic parameters and crosstalk dynamics of ship harness[J]. China Ship Repair, 2013, 26(2): 34-36, 40 doi: 10.3969/j.issn.1001-8328.2013.02.014 [10] 高印寒, 王天皓, 杨开宇, 等. 汽车线束的动态串扰特性预测[J]. 吉林大学学报(工学版), 2014, 44(5):1258-1263. (Gao Yinhan, Wang Tianhao, Yang Kaiyu, et al. Prediction of the dynamic crosstalk of automotive wiring hardness[J]. Journal of Jilin University (Engineering and Technology Edition), 2014, 44(5): 1258-1263 [11] 张昭, 王世山, 赵亮, 等. 多导体线束内串扰概率分布的预测[J]. 电工技术学报, 2017, 32(7):204-214. (Zhang Zhao, Wang Shishan, Zhao Liang, et al. Prediction of crosstalk probability distribution in cable bundles[J]. Transactions of China Electrotechnical Society, 2017, 32(7): 204-214 [12] Huang Chao, Zhao Yang, Yan Wei, et al. A new method for predicting crosstalk of random cable bundle based on BAS-BP neural network algorithm[J]. IEEE Access, 2020, 8: 20224-20232. doi: 10.1109/ACCESS.2020.2969221 [13] 袁赛, 赵海龙, 曹海丽, 等. 基于MATLAB的带噪图像的中值滤波[J]. 电子世界, 2016(18):185. (Yuan Sai, Zhao Hailong, Cao Haili, et al. Median filtering of noisy image based on MATLAB[J]. Electronics World, 2016(18): 185 doi: 10.3969/j.issn.1003-0522.2016.18.132 [14] 曾建华, 黄时杰. 典型图像边缘检测算子的比较与分析[J]. 河北师范大学学报(自然科学版), 2020, 44(4):295-301. (Zeng Jianhua, Huang Shijie. Comparison and analysis of typical image edge detection operators[J]. Journal of Hebei Normal University (Natural Science Edition), 2020, 44(4): 295-301 [15] 纪执楷, 赵翔, 闫丽萍, 等. 随机线缆束的建模及其场线耦合分析[J]. 无线电工程, 2018, 48(9):788-792. (Ji Zhikai, Zhao Xiang, Yan Liping, et al. Modeling of random wire-bundle and analysis of field-to-line coupling[J]. Radio Engineering, 2018, 48(9): 788-792 doi: 10.3969/j.issn.1003-3106.2018.09.14 [16] Danielson G. On finding the simple paths and circuits in a graph[J]. IEEE Transactions on Circuit Theory, 1968, 15(3): 294-295. doi: 10.1109/TCT.1968.1082837 [17] Paul C R. Analysis of multiconductor transmission lines[M]. New York: John Wiley & Sons, Inc., 1994. [18] Haase H. Full-wave field interactions of nonuniform transmission lines[D]. Magdeburg: der Otto-von-Guericke-Universität Magdeburg, 2005. [19] 弗雷德里卡·M·特奇, 米歇尔·V·艾诺茨, 托比杰恩·卡尔松. EMC分析方法与计算模型[M]. 吕英华, 王旭莹, 译. 北京: 北京邮电大学出版社, 2009: 348-350.Tesche F M, Ianoz M V, Karlsson T. EMC analysis methods and computational models[M]. Lü Yinghua, Wang Xuying, trans. Beijing: Beijing University of Posts and Telecommunications Press, 2009: 348-350 [20] 谢处方, 饶克谨. 电磁场与电磁波[M]. 北京: 高等教育出版社, 2003: 67-68.Xie Chufang, Rao Kejin. Electromagnetic field and electromagnetic wave[M]. Beijing: Higher Education Press, 2003: 67-68 期刊类型引用(3)
1. 何进，赵士龙，郭凌伟，华有杰，叶仁广，徐时清. 反射式体布拉格光栅设计及其在激光器中的应用研究. 光电子·激光. 2024(02): 122-127 . 
百度学术2. 李奎奎，张伟桥，董风鑫，张建心，周旭彦. 高功率半导体激光叠阵光纤耦合系统与光束整形设计. 中国激光. 2024(17): 98-103 . 百度学术3. 任慧娟. 激光惯性约束下多帧影视图像视觉传达模型构建. 激光杂志. 2023(09): 103-108 . 百度学术其他类型引用(1)
-
下载:
下载:
点击查看大图
计量
- 文章访问数: 1096
- HTML全文浏览量: 318
- PDF下载量: 64
- 被引次数: 4
