《万方数据-数字化期刊群》全文上网期刊
CNKI《中国学术期刊(网络版)》全文收录期刊
《中文科技期刊数据库》(维普网)全文收录期刊
超星期刊域出版平台、博看网全文收录期刊
日本科学技术振兴机构数据库收录

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析

冯桃君 焦子龙 姜利祥 郑慧奇 姜海富 刘宇明 王鹏 刘学超 李涛

冯桃君, 焦子龙, 姜利祥, 等. 宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析[J]. 航天器环境工程, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
引用本文: 冯桃君, 焦子龙, 姜利祥, 等. 宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析[J]. 航天器环境工程, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
FENG T J, JIAO Z L, JIANG L X, et al. Method for information extraction of a cosmic dust trajectory sensor and analysis of influencing factors on inversion precision[J]. Spacecraft Environment Engineering, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
Citation: FENG T J, JIAO Z L, JIANG L X, et al. Method for information extraction of a cosmic dust trajectory sensor and analysis of influencing factors on inversion precision[J]. Spacecraft Environment Engineering, 2023, 40(6): 692-701 doi:  10.12126/see.2023024

宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析

doi: 10.12126/see.2023024
基金项目: 中国航天科技集团有限公司自主研发项目
详细信息
    作者简介:

    冯桃君,工程师,主要从事空间环境及反演算法研究

  • 中图分类号: V242.2; V423.6

Method for information extraction of a cosmic dust trajectory sensor and analysis of influencing factors on inversion process

  • 摘要: 利用基于多层金属丝栅网的尘埃轨迹探测器(DTS)可获得宇宙尘埃粒子的速度及电荷信息,对确定粒子起源、揭示带电尘埃粒子与行星际空间环境的相互作用具有重要意义。文章针对利用DTS探测信号提取尘埃粒子速度和电荷信息的方法进行研究,并仿真分析反演精度的影响因素。结果表明,尘埃粒子电荷量和速度大小的探测误差随尘埃粒子速度的增大而增大,随探测器电荷噪声比(QNR)的增大而减小,不受尘埃粒子入射角θy和尘埃带电量的影响;当QNR大于10时,对尘埃粒子的电荷量、速度和运动方向的反演精度分别优于1.5%、0.8%和0.6°。以上结果可为DTS设计及其后续空间应用提供参考。
  • 图  1  DTS结构示意[5]

    Figure  1.  Structure schematic of DTS[5]

    图  2  带电尘埃粒子在DTS中的穿越轨迹示意[5]

    Figure  2.  Trajectory schematic of charged particles travelling in DTS[5]

    图  3  DTS数值模型

    Figure  3.  Numerical model of DTS

    图  4  Auer在尘埃实验室记录到的DTS信号[1]

    Figure  4.  DTS signals recorded in the dust laboratory by Auer[1]

    图  5  参数提取反演流程

    Figure  5.  Inversion process for parameter extraction

    图  6  尘埃粒子从(13, -22, 100)到(-7, 38, -100)穿越DTS的模拟信号

    Figure  6.  Simulated signals generated by dust particles travelling from (13, -22, 100) to (-7, 38, -100) through DTS

    图  7  χ2随迭代次数的变化

    Figure  7.  Variation of χ2 with iteration steps

    图  8  不同参数估计对应的金属丝信号

    Figure  8.  Metal wire signals corresponding to different parameter estimates

    图  9  反演误差均值及其绝对值均值随vz的变化

    Figure  9.  Variations of mean inversion error and mean absolute inversion error with vz

    图  10  反演误差均值及其绝对值均值随θy的变化

    Figure  10.  Variations of mean inversion error and mean absolute inversion error with θy

    图  11  反演误差均值及其绝对值均值随Q的变化

    Figure  11.  Variations of mean inversion error and mean absolute inversion error with Q

    图  12  反演误差均值及其绝对值均值随QNR的变化

    Figure  12.  Variations of mean inversion error and mean absolute inversion error with QNR

    表  1  尘埃电荷与轨迹的反演结果

    Table  1.   Inversion results of the dust charge and trajectory

    参数真值反演值绝对误差相对误差/%
    Q/C1.000.99-0.01-1.00
    v/(km·s-1)5.244.84-0.40-7.63
    θx/(°)5.715.800.091.58
    θy/(°)-16.70-16.71-0.01-0.06
    下载: 导出CSV

    表  2  各控制变量的取值

    Table  2.   Values for control variables

    控制变量 取值
    vz/(km·s-1) 2.5、5.0、10.0、20.0、40.0、80.0
    θy/(°) 0、10、20、30
    Q/C 1、2、3、4、5
    QNR 10、20、30、40、50、60、70、80
    下载: 导出CSV

    表  3  不同θy对应的尘埃轨迹

    Table  3.   Dust trajectories for different values of θy

    θy/(°) xin, yin, zin/mm xout, yout, zout/mm
    0 10.0, 16.0, 100.0 10.0, 16.0, -100.0
    10 10.0, 30.1, 100.0 10.0, -5.2, -100.0
    20 10.0, 45.1, 100.0 10.0, -27.7, -100.0
    30 10.0, 62.2, 100.0 10.0, -53.3, -100.0
    下载: 导出CSV
  • [1] AUER S, GRÜN E, KEMPF S, et al. Characteristics of a dust trajectory sensor[J]. Review of Scientific Instruments, 2008, 79: 084501 doi:  10.1063/1.2960566
    [2] AUER S, LAWRENCE G, GRÜN E, et al. A self-triggered dust trajectory sensor[J]. Nuclear Instruments and Methods in Physics Research, 2010, A622: 74-82
    [3] MANN I, KIMURA H, BIESECHER D A, et al. Dust near the sun[J]. Space Science Reviews, 2004, 110: 269-305 doi:  10.1023/B:SPAC.0000023440.82735.ba
    [4] KRÜGER H, STRUB P, SRAMA R, et al. Modelling DESTINY+ interplanetary and interstellar dust measurements en route to the active asteroid (3200) Phaethon[J]. Planetary and Space Science, 2019, 172: 22-42 doi:  10.1016/j.pss.2019.04.005
    [5] XIE J, STERNOVSKY Z, GRÜN E, et al. Dust trajectory sensor: accuracy and data analysis[J]. Review of Scientific Instruments, 2011, 82: 105104 doi:  10.1063/1.3646528
    [6] 姜利祥, 焦子龙, 郑慧奇, 等. 地月平动点尘埃动力学特性原位探测方法研究[J]. 航天器环境工程, 2022, 39(4): 413-418 doi:  10.12126/see.2022.04.013

    JIANG L X, JIAO Z L, ZHENG H Q, et al. Method for in-situ detection of dynamical properties of dust cloud particles at Earth-Moon liberation point[J]. Spacecraft Environment Engineering, 2022, 39(4): 413-418 doi:  10.12126/see.2022.04.013
    [7] GRÜN E, GUSTAFSON B A S, DERMOTT S, et al. Interplanetary dust[M]. Berlin: Springer Verlag, 2001: 324-330
    [8] 焦子龙, 姜利祥, 李涛, 等. 宇宙尘星载原位探测技术研究[J]. 装备环境工程, 2018, 15(5): 72-78

    JIAO Z L, JIANG L X, LI T, et al. Spacebone in-situ detection technologies of cosmic dust[J]. Equipment Environmental Engineering, 2018, 15(5): 72-78
    [9] KHALISI E, SRAMA R, GRÜN E. Counter data of the cosmic dust analyzer aboard the Cassini spacecraft and possible “dust cloud” at Saturn[J]. Advances in Space Research, 2015, 55: 303-310 doi:  10.1016/j.asr.2014.09.002
    [10] MALASPINA D M, O’BRIEN L E, THAYER F, et al. Revisiting STEREO interplanetary and interstellar dust flux and mass estimates[J]. Journal of Geophysical Research: Space Physics, 2015, 120: 6085-6100 doi:  10.1002/2015JA021352
    [11] WOOD S R, MALASPINA D M, ANDERSSON L, et al. Hypervelocity dust impacts on the Wind spacecraft: correlations between Ulysses and Wind interstellar dust detections[J]. Journal of Geophysical Research: Space Physics, 2015, 120: 7121-7129 doi:  10.1002/2015JA021463
    [12] GRÜN E, SRAMA R, HELFERT S, et al. Prospects of dust astronomy missions[C]//Workshop on Dust in Planetary Systems (ESA SP-643). Kauai, HI, USA, 2005: 245-249
    [13] STERNOVSKY Z, GRÜN E, DRAKE K, et al. Novel instrument for dust astronomy: dust telescope[C]//IEEE Aerospace Conference. Big Sky, MT, USA, 2011
    [14] SRAMA R, SROWIG A, RACHEV M, et al. Development of an advanced dust telescope[J]. Earth, Moon, Planets, 2005, 95: 211-220
    [15] STERNOVSKY Z, AMYX K, BANO G, et al. Large area mass analyzer instrument for the chemical analysis of interstellar dust particles[J]. Review of Scientific Instruments, 2007, 78: 014501 doi:  10.1063/1.2431089
    [16] LI Y W, KEMPF S, SIMOLKA J, et al. Instrument concept of a single channel dust trajectory detector[J]. Advances in Space Research, 2017, 59(6): 1636-1641 doi:  10.1016/j.asr.2016.12.037
    [17] VORONOV K E, TELEGIN A M, RYAZANOV D M. The concept of a grid-based instrument for measuring micrometeoroid parameters[J]. Applied Physics, 2021, 4: 73-80
  • [1] 王领华, 张皓, 王骞, 吕建伟, 吴勇, 赵允宁, 刘辰.  航天器动力管路热控设计与试验研究 . 航天器环境工程, 2023, 40(1): 1-5. doi: 10.12126/see.2022093
    [2] 董洪建, 李小军, 封国宝, 李韵.  基于光子晶体的太赫兹环行器设计及其内带电效应研究 . 航天器环境工程, 2022, 39(2): 161-168. doi: 10.12126/see.2022.02.007
    [3] 陈旭, 杨晓宁, 王晶, 秦家勇, 毕研强, 蒋山平.  基于太阳模拟器的网状天线热变形仿真与试验研究 . 航天器环境工程, 2022, 39(2): 140-147. doi: 10.12126/see.2022.02.004
    [4] 刘岩, 陈瑞勋, 孙犇.  一种基于单机EMC试验结果的航天器系统级电磁兼容性分析方法 . 航天器环境工程, 2021, 38(4): 395-400. doi: 10.12126/see.2021.04.003
    [5] 郑雪松, 丁丽娜, 王文炎, 张洪伟, 罗磊, 刘婧宇, 吉俐, 王荣飞, 刘贲.  空间用长寿命自浮动电连接器原子氧侵蚀效应仿真研究 . 航天器环境工程, 2021, 38(4): 407-412. doi: 10.12126/see.2021.04.005
    [6] 焦子龙, 赵明, 姜利祥, 李涛, 刘宇明, 姜海富, 孙继鹏, 郑慧奇, 朱云飞.  一种航天器空间环境效应仿真分析软件架构 . 航天器环境工程, 2021, 38(5): 555-562. doi: 10.12126/see.2021.05.010
    [7] 陈锴, 耿兴宁, 李吉宁, 钟凯, 徐德刚, 蒋山平, 张景川, 姚建铨.  太赫兹波在高速飞行器等离子体鞘套中的传输特性 . 航天器环境工程, 2020, 37(5): 421-427. doi: 10.12126/see.2020.05.001
    [8] 高奔, 韩佩彤, 石光达, 盆洪民, 冉江南, 李祺.  “嫦娥五号”月球探测器着陆姿态模拟装置设计 . 航天器环境工程, 2020, 37(1): 65-72. doi: 10.12126/see.2020.01.011
    [9] 丁亮, 彭毓川, 郑慧奇, 唐振宇, 任琼英, 赵华, 秦玮.  基于环境工质的超低轨吸气式螺旋波电推进器仿真分析 . 航天器环境工程, 2020, 37(4): 330-335. doi: 10.12126/see.2020.04.003
    [10] 谷运双, 霍岩.  载人航天器密封舱内强制通风下的火灾温度场及流场特性仿真分析 . 航天器环境工程, 2020, 37(1): 31-36. doi: 10.12126/see.2020.01.005
    [11] 邢卓异, 王彤, 舒燕, 张旭辉.  地外天体上航天器起飞瞬时羽流对主发动机干扰效应影响 . 航天器环境工程, 2019, 36(6): 565-570. doi: 10.12126/see.2019.06.006
    [12] 韩继广, 陶晶亮, 盖照亮, 周国锋, 陈丽, 彭光东.  航天器热平衡试验用大面阵外热流动态模拟系统设计及应用验证 . 航天器环境工程, 2019, 36(5): 495-501. doi: 10.12126/see.2019.05.014
    [13] 高海洋, 王栋, 李栋, 刘明辉, 韩晓健, 樊世超.  航天器微振动工装动力学特性仿真分析与试验验证 . 航天器环境工程, 2019, 36(3): 252-256. doi: 10.12126/see.2019.03.009
    [14] 张彧, 赵吉喆, 张翔, 刘冈云.  环月轨道一体式星敏感器热设计及仿真验证 . 航天器环境工程, 2016, 33(5): 516-520. doi: 10.12126/see.2016.05.010
    [15] 邓佳欣, 易忠, 孟立飞, 陈金刚, 唐小金, 王松.  卫星磁力矩器磁场效应仿真分析 . 航天器环境工程, 2016, 33(1): 52-57. doi: 10.12126/see.2016.01.009
    [16] 刘金龙1, 郄殿福1,2, 裴一飞1, 陶涛1,2使用太阳模拟器进行热平衡试验的附加外热流影响分析 . 航天器环境工程, 2015, 32(2): 182-186. doi: 10.12126/see.2015.02.007
    [17] 李唯丹1, 闫荣鑫1,2, 孙伟1, 綦磊1载人航天器内超声在轨检漏的应用基础研究 . 航天器环境工程, 2014, 31(2): 178-181. doi: 10.12126/see.2014.02.013
    [18] 王伟臣, 石泳, 张健, 郭娟.  基于羽流效应分析的航天器发动机布局设计研究 . 航天器环境工程, 2014, 31(4): 374-377. doi: 10.12126/see.2014.04.006
    [19] 魏传锋, 张伟, 曹剑峰, 郝平.  载人航天器密封舱噪声控制与试验 . 航天器环境工程, 2013, 30(1): 91-93.
    [20] 叶红, 王志瑾.  高超声速飞行器热防护结构参数优化及对比分析 . 航天器环境工程, 2013, 30(5): 516-521.
  • 加载中
图(12) / 表ll (3)
计量
  • PDF下载量:  6
  • 文章访问数:  40
  • HTML全文浏览量:  3
文章相关
  • 中图分类号:  V242.2; V423.6
  • 收稿日期:  2023-03-02
  • 修回日期:  2023-11-20
  • 网络出版日期:  2023-12-28
  • 刊出日期:  2023-12-25

宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析

doi: 10.12126/see.2023024
    基金项目:  中国航天科技集团有限公司自主研发项目
    作者简介:

    冯桃君,工程师,主要从事空间环境及反演算法研究

  • 中图分类号: V242.2; V423.6

摘要: 利用基于多层金属丝栅网的尘埃轨迹探测器(DTS)可获得宇宙尘埃粒子的速度及电荷信息,对确定粒子起源、揭示带电尘埃粒子与行星际空间环境的相互作用具有重要意义。文章针对利用DTS探测信号提取尘埃粒子速度和电荷信息的方法进行研究,并仿真分析反演精度的影响因素。结果表明,尘埃粒子电荷量和速度大小的探测误差随尘埃粒子速度的增大而增大,随探测器电荷噪声比(QNR)的增大而减小,不受尘埃粒子入射角θy和尘埃带电量的影响;当QNR大于10时,对尘埃粒子的电荷量、速度和运动方向的反演精度分别优于1.5%、0.8%和0.6°。以上结果可为DTS设计及其后续空间应用提供参考。

English Abstract

冯桃君, 焦子龙, 姜利祥, 等. 宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析[J]. 航天器环境工程, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
引用本文: 冯桃君, 焦子龙, 姜利祥, 等. 宇宙尘埃轨迹探测器信号提取方法及反演精度影响因素分析[J]. 航天器环境工程, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
FENG T J, JIAO Z L, JIANG L X, et al. Method for information extraction of a cosmic dust trajectory sensor and analysis of influencing factors on inversion precision[J]. Spacecraft Environment Engineering, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
Citation: FENG T J, JIAO Z L, JIANG L X, et al. Method for information extraction of a cosmic dust trajectory sensor and analysis of influencing factors on inversion precision[J]. Spacecraft Environment Engineering, 2023, 40(6): 692-701 doi:  10.12126/see.2023024
    • 微米和亚微米级的宇宙尘埃粒子通常来自彗星、小行星、月球、星体间的碰撞碎片和行星际物质[1-2]。宇宙尘埃暴露在等离子体、太阳紫外辐射、高能粒子、电磁场等组成的复杂空间环境中,运动过程中会携带电荷,因此尘埃轨迹是行星和星际磁场、太阳和行星引力、等离子体阻力及辐射压力等因素共同作用的结果[3-4]。Xie等[5]综述并研究宇宙尘埃的运动方式发现:小行星尘埃由于Poynting-Robertson效应会以螺旋的方式缓慢接近太阳;β-流星体主要受到太阳辐射压的作用,以双曲线轨迹逃逸出太阳系;Ulysses号航天器上的尘埃探测器发现木星尘埃粒子受到磁层电场的加速并以大于100 km/s的速度逃逸出木星。并且指出,宇宙尘埃的轨迹和电荷信息有助于确定粒子起源,辨别不同起源的粒子,揭示带电尘埃粒子与行星际空间环境的相互作用。

      对宇宙尘埃粒子进行探测的方法很多是利用尘埃探测器的原位测量[6-8],高速粒子在探测器内部发生碰撞电离,可通过电离产生的总电荷量和电脉冲信号的上升时间分别推算出粒子质量和碰撞速度;并通过对粒子飞行时间的分析得出粒子的化学成分。该项技术已在Vega、Wind、Cassini、STEREO等航天器上用于探测尘埃的质量、成分和速度大小[9-11],但速度方向只能通过统计数据分析得到,要求探测器进行全方位扫描,限制了碰撞探测器的应用。Grün 等[12]提出了尘埃天文学概念,涉及宇宙尘埃粒子的起源以及粒子化学成分、电荷、速度、轨迹的测量。尘埃轨迹探测器(dust trajectory sensor, DTS)是尘埃天文学的关键探测组件,当带电粒子穿过DTS内部的金属丝阵列时,会产生诱导电荷,可通过分析金属丝输出的电信号提取粒子的电荷量和轨迹信息[1, 13]。DTS与不同的仪器结合可构成尘埃望远镜[14]、离子质量分析器[15]和静电质量探测器[1]等,适用于不同速度、不同尺寸的尘埃颗粒探测。

      DTS在几十年前就已被提出[5],其几何构型已经过不断优化。Auer等[1-2]利用COULOMB软件对不同几何构型的DTS进行数值仿真,分析金属丝的几何参数、间距及数量等对电荷分布、信号灵敏度的影响,获得了适用于速度在2~5 km/s的微米级尘埃探测的DTS结构参数;并得出DTS的电荷噪声比(charge to noise ratio, QNR)≥6.25才能从背景噪声中探测到尘埃信号。Xie等[5]同样利用COULOMB软件,发现探测器的QNR越小电荷量探测误差分布范围越大,当QNR=10时,尘埃电荷量探测误差为-1.5%~5%。为提高仪器对尘埃电荷量的探测灵敏度,Li等[16]提出了分段式低电容金属丝阵列结构,使尘埃电荷量探测阈值降至0.2 fC。Voronov等[17]介绍了一种基于网格面和金属丝阵列面交替的DTS设计,不仅可以测量带电微流星体的电荷量、轨迹,还能计算微流星体的质量。

      宇宙尘埃大都分布在亚微米到厘米量级,速度在每秒几千米到几十千米不等,DTS的探测精度也因尘埃性质不同而显示差异。本文基于典型的DTS结构,针对尘埃轨迹和电荷量的测量建立DTS数值仿真模型,研究从DTS金属丝信号提取尘埃粒子速度矢量和电荷量的方法,重点分析尘埃速度、带电量和入射方向对算法精度的影响,有助于探索适用于DTS探测的尘埃范围。

    • DTS一般由中间4个金属丝阵列面和分别位于顶部和底部的2个屏蔽栅网组成,每个金属丝阵列面包含7~16条相互平行的金属丝,两相邻平面的金属丝方向正交。如图1[5]所示,两屏蔽栅网接地;每条金属丝都独自与一个电荷灵敏放大器(charge sensitive amplifier, CSA)连接,CSA通道的输出信号被瞬态记录器以一定频率采集。当带电尘埃粒子穿过DTS时,会产生诱导电荷。根据同一平面金属丝信号特性可确定粒子穿过当前平面的一维坐标(x1y2x3y4),如图2[5]所示,进而可推导出尘埃粒子的速度矢量和电荷量[1, 5]

      图  1  DTS结构示意[5]

      Figure 1.  Structure schematic of DTS[5]

      图  2  带电尘埃粒子在DTS中的穿越轨迹示意[5]

      Figure 2.  Trajectory schematic of charged particles travelling in DTS[5]

      本文研究所用DTS模型含4个平行于xy平面的金属丝阵列面,每个平面包含7条金属丝,相邻平面间隔40 mm,金属丝直径0.4 mm,长140 mm,间隔20 mm; DTS的整体尺寸为160 mm×160 mm×200 mm(长×宽×高),如图3所示。平面编号从上到下依次增加,平面1、3的金属丝平行于y轴,平面2、4的金属丝平行于x轴。模型坐标原点位于结构的几何中心,图中以红点标记。

      图  3  DTS数值模型

      Figure 3.  Numerical model of DTS

      探测器金属丝的诱导电荷量与尘埃粒子的入射位置、角度、探测器尺寸和构型有关。Auer等[1]给出了当带电量为Q的球形尘埃颗粒穿过DTS时,金属丝i的诱导电荷qi的近似表达式

      $$ {q}_{i}=\frac{{\text{-}}(Q/{r}_{i})}{\sum _{j}(1/{r}_{j})} \text{,} $$ (1)

      式中:ri为尘埃到金属丝i的垂直距离;$ \sum _{j}(1/{r}_{j}) $代表求和遍历了DTS所有的金属丝。尘埃的位置(xp, yp, zp)随时间变化,因此尘埃到金属丝的距离和金属丝的诱导电荷均为时间的函数。尘埃电荷量可由所有金属丝诱导电荷相加近似得到。

      Auer等[1]在地面试验中采集到的随时间变化的DTS输出信号如图4所示,图中信号来自两金属丝阵列平面,每阵列包含16条金属丝。输出信号的形状反映金属丝与尘埃粒子的距离,信号强度反映金属丝带电量,信号峰值表明尘埃粒子近距离地从某金属丝旁穿过该金属丝所在平面。

      图  4  Auer在尘埃实验室记录到的DTS信号[1]

      Figure 4.  DTS signals recorded in the dust laboratory by Auer[1]

      假设高速尘埃粒子在DTS内的轨迹为一条直线,由尘埃粒子进、出DTS的时间(tintout)和位置坐标((xin, yin, zin)、(xout, yout, zout))可以确定尘埃粒子的速度和方向。通过式(2)可由4个一维坐标确定尘埃粒子进、出DTS的坐标(xin, yin)、(xout, yout) [5]

      $$ \left\{\begin{aligned} & {x}_{\mathrm{i}\mathrm{n}}=\frac{3}{2}{x}_{1}-\frac{1}{2}{x}_{3}\\ & {x}_{\mathrm{o}\mathrm{u}\mathrm{t}}=2{x}_{3}-{x}_{1}\\ & {y}_{\mathrm{i}\mathrm{n}}=2{y}_{2}-{y}_{4}\\ & {y}_{\mathrm{o}\mathrm{u}\mathrm{t}}=\frac{3}{2}{y}_{4}-\frac{1}{2}{y}_{2}\end{aligned}\right.。 $$ (2)

      粒子速度的大小为[5]

      $$ v=\frac{\sqrt{{({x}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{x}_{\mathrm{i}\mathrm{n}})}^{2}+{({y}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{y}_{\mathrm{i}\mathrm{n}})}^{2}+{({\textit{z}}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{\textit{z}}_{\mathrm{i}\mathrm{n}})}^{2}}}{{t}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{t}_{\mathrm{i}\mathrm{n}}} 。 $$ (3)

      粒子入射角为[5]

      $$ \left\{\begin{aligned} & {\theta }_{x}=\arctan\left(\frac{{x}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{x}_{\mathrm{i}\mathrm{n}}}{{\textit{z}}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{\textit{z}}_{\mathrm{i}\mathrm{n}}}\right)\\ & {\theta }_{y}=\arctan\left(\frac{{y}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{y}_{\mathrm{i}\mathrm{n}}}{{{\textit{z}}}_{\mathrm{o}\mathrm{u}\mathrm{t}}-{\textit{z}}_{\mathrm{i}\mathrm{n}}}\right)\end{aligned}\right., $$ (4)

      其中,θxθy分别为尘埃轨迹在xz、yz平面的投影与z轴的夹角。

    • 为提取尘埃粒子的速度矢量和电荷量,在已知探测器几何构型的前提下,表征尘埃轨迹和电荷量需要7个独立参量:尘埃电荷量(Q);尘埃进、出DTS的时刻(tintout);尘埃穿过4个金属丝阵列平面的一维坐标(x1y2x3y4),故设参数矢量为P=(Q, tin, tout, x1, y2, x3, y4)。文献[5]提出了一种基于探测数据与仿真数据拟合的参数提取方法,其反演流程如图5所示。

      图  5  参数提取反演流程

      Figure 5.  Inversion process for parameter extraction

      参数提取的具体步骤为:

      1)初始化参数矢量P0。尘埃初始电荷量Q可设置为DTS所有金属丝电荷量之和;x1y2x3y4根据式(1)估算得到。文献[5]中并未介绍tintout的取值方法,本文设DTS的采样频率为r(单位为MS/s,即106 sample/s),对应的信号采样个数为n,并将第1个采样信号的时刻t1设置为0作为参考,则tin∈[-1/r, 0]、tout∈[(n–1)/r, n/r],单位为μs,tintout的初值在各自范围内均匀随机取值。

      2)初始化参数的误差范围ΔP0。ΔQ为±10%,Δt为±1/r,Δx或Δy为±10 mm。

      3)数据拟合。从P0±ΔP0的范围中随机选取一组 P=(Q, tin, tout, x1, y2, x3, y4),根据式(1)以及DTS采样频率仿真计算DTS各金属丝的采样信号,然后计算仿真信号相对实测信号的误差,

      $$ {\chi }^{2}=\sum _{N}{({D}_{m}-{S}_{m})}^{2} 。 $$ (5)

      式中$ : $χ2为仿真信号误差;N为数据总个数;Dm为第m个实测数据;Sm为第m个仿真数据。

      重复上述随机选取参数Pχ2计算103次,χ2值随P变化,χ2最小值对应的参数集为P1, min

      4)优化参数矢量P1。由于χ2也是单一参数的函数,依次改变P1,min中第ll=1, 2,$ \cdots $, 7)个参数P(l)的值,使χ2最小的参数值成为该参数的新估计值P1(l)。

      5)缩小参数的误差范围ΔP1。ΔP1取2ΔP0/3和3$ \left|{\boldsymbol{P}}_{1}-{\boldsymbol{P}}_{0}\right| $中的大者。

      6)重复步骤3)~5),直到χ2的变化小于0.1%。

    • 为验证2.1节尘埃参数提取方法的有效性,针对典型算例进行计算分析。尘埃预设轨迹为入射点坐标(13, -22, 100)、出射点坐标(-7, 38, -100),入射角θx=-5.7°、θy=16.7°,z方向速度分量vz=5 km/s,尘埃穿过各金属丝平面的一维坐标为x1=9 mm,y2=2 mm,x3=1 mm,y4=26 mm,涵盖了尘埃与金属丝之间的近距离、中等距离和远距离。另外,金属丝输出信号用诱导电荷与尘埃电荷之比表示,设尘埃电荷为1 C。设置QNR=10,采样频率为r=10 MS/s,尘埃穿过DTS的时间为t=40 μs,因此每个金属丝有n=400个采样信号。设尘埃进入DTS的时间tin=0,则tout=40 μs,参数真值Ptrue=(1, 0, 40, 9, 2, 1, 26),本文用下标“opt”和“true”分别表示参数的最优估计值和真值。离尘埃近的金属丝将聚集大部分诱导电荷,因此在信号分析中,为节省计算时间只对每个平面中离尘埃最近的2条金属丝的信号进行拟合,这样用于拟合的信号一共有3200个。实际的观测信号包含了电荷灵敏放大器产生的随机噪声,为模拟真实信号,在式(1)计算出的电荷信号中添加正态分布的随机噪声,其均值为μ=0,标准差σ=1/QNR。根据上述仿真条件,随尘埃位置zp变化的无噪声金属丝模拟信号如图6(a)所示,图中虚线表示各平面对应的z轴坐标,顶端数字为平面序号。从图中可看出,当尘埃穿过平面时与金属丝距离很近(如y2=2 mm,x3=1 mm),该金属丝信号幅值骤升,随后下降;与此同时,尘埃另一侧相距较远的金属丝信号出现下凹。当尘埃穿过平面时与两侧金属丝距离相差不大(如x1=9 mm),两金属丝信号形状相似。图6(b)所示是叠加噪声后的金属丝信号,用于模拟DTS的实测信号。

      图  6  尘埃粒子从(13, -22, 100)到(-7, 38, -100)穿越DTS的模拟信号

      Figure 6.  Simulated signals generated by dust particles travelling from (13, -22, 100) to (-7, 38, -100) through DTS

      按上述尘埃轨迹与电荷信号的分析方法,经过6次迭代得到各参数的最优估计Popt=(0.99, -0.05, 41.27, 9.04, 2.06, 0.92, 26.07),χ2随迭代次数的变化如图7所示。参数估计偏差ΔP=PoptPtrue=(-0.01, -0.05, 1.27, 0.04, 0.06, -0.08, 0.07),尘埃电荷和轨迹的反演结果如表1所示。其中,速度的反演结果误差最大,被低估了7.63%,其原因可能是在独立参数的估计中算法能获得精度较高的尘埃位置,但尘埃穿过DTS的时间估值偏差较大(1.32 μs),导致反演的尘埃速度减小。

      图  7  χ2随迭代次数的变化

      Figure 7.  Variation of χ2 with iteration steps

      表 1  尘埃电荷与轨迹的反演结果

      Table 1.  Inversion results of the dust charge and trajectory

      参数真值反演值绝对误差相对误差/%
      Q/C1.000.99-0.01-1.00
      v/(km·s-1)5.244.84-0.40-7.63
      θx/(°)5.715.800.091.58
      θy/(°)-16.70-16.71-0.01-0.06

      参数初值P0与最优估计Popt对应的金属丝信号分别如图8(a)和(b)所示。

      图  8  不同参数估计对应的金属丝信号

      Figure 8.  Metal wire signals corresponding to different parameter estimates

      图中蓝色表示模拟的实际观测数据,红色表示根据P0Popt反演出的金属丝信号。可以看到,经过6次迭代后,Popt能复现实际的观测信号,说明时间估计的误差只影响尘埃速度大小的反演结果,对尘埃速度方向及DTS金属丝信号的反演结果没有影响。

    • 分别改变用于仿真的尘埃速度分量vz、入射角度θy、尘埃电荷量Q和探测器电荷噪声比QNR,分析它们对反演精度的影响。各参数取值如表2所示,当改变其中1个参量时,其他3个参量按表2中粗体数字取值。根据采用的DTS结构,尘埃穿过DTS的最大入射角为38.6°,仿真设置θx=0°,θy在0°~30°范围内变化,表3列出了与表2θy对应的预设尘埃轨迹(进、出DTS的位置坐标)。

      表 2  各控制变量的取值

      Table 2.  Values for control variables

      控制变量 取值
      vz/(km·s-1) 2.5、5.0、10.0、20.0、40.0、80.0
      θy/(°) 0、10、20、30
      Q/C 1、2、3、4、5
      QNR 10、20、30、40、50、60、70、80

      表 3  不同θy对应的尘埃轨迹

      Table 3.  Dust trajectories for different values of θy

      θy/(°) xin, yin, zin/mm xout, yout, zout/mm
      0 10.0, 16.0, 100.0 10.0, 16.0, -100.0
      10 10.0, 30.1, 100.0 10.0, -5.2, -100.0
      20 10.0, 45.1, 100.0 10.0, -27.7, -100.0
      30 10.0, 62.2, 100.0 10.0, -53.3, -100.0

      表2依次改变各控制变量数值,针对每一组控制变量(vzθyQ、QNR),重复图5的反演流程100次,每次的模拟信号叠加不同的噪声,会反演得到100组不同的尘埃电荷量(Q)及轨迹参数(vzθyθx)。尘埃入射角度的反演误差Δθ为反演最优值与真值之差,即

      $$ \mathrm{\Delta }\theta ={\theta }_{\mathrm{o}\mathrm{p}\mathrm{t}}-{\theta }_{\mathrm{t}\mathrm{r}\mathrm{u}\mathrm{e}} \text{;} $$ (6)

      电荷量和速度大小的反演误差Δa以相对误差表征,即

      $$ \mathrm{\Delta }a=\frac{{a}_{\mathrm{o}\mathrm{p}\mathrm{t}}-{a}_{\mathrm{t}\mathrm{r}\mathrm{u}\mathrm{e}}}{{a}_{\mathrm{t}\mathrm{r}\mathrm{u}\mathrm{e}}}\times 100{\text{%}} \text{,} $$ (7)

      其中符号a代表参数Qv

      计算得到100组反演结果的反演误差后,分析这100组反演误差的均值及误差绝对值的均值随表2中各控制变量的变化,结果如图9图12所示。其中,误差绝对值的均值(图中蓝线)反映反演结果相对 “真”值的整体偏离幅度,误差均值(图中粉线)反映反演结果相对“真”值的整体偏离方向(被高估或被低估)。

      图  9  反演误差均值及其绝对值均值随vz的变化

      Figure 9.  Variations of mean inversion error and mean absolute inversion error with vz

      图  10  反演误差均值及其绝对值均值随θy的变化

      Figure 10.  Variations of mean inversion error and mean absolute inversion error with θy

      图  11  反演误差均值及其绝对值均值随Q的变化

      Figure 11.  Variations of mean inversion error and mean absolute inversion error with Q

      图  12  反演误差均值及其绝对值均值随QNR的变化

      Figure 12.  Variations of mean inversion error and mean absolute inversion error with QNR

    • 反演误差随vz的变化如图9所示。

      图9可以看出,尘埃电荷量和轨迹的反演误差随尘埃速度的增大而增大(蓝线),这是由于随着尘埃速度的增大,尘埃穿过DTS的时间缩短,导致探测器的采样信号数量减少。当尘埃速度小于20 km/s时,误差随速度的变化率较大;当尘埃速度超过20 km/s时,误差的增长变缓。随着尘埃速度从2.5 km/s增加到80 km/s,$ \overline{\left|\mathrm{\Delta }Q\right|} $从0.25%增大到1.25%,$ \overline{\left|\mathrm{\Delta }v\right|} $从0.12%增大到0.74%,$ \overline{\left|\mathrm{\Delta }{\theta }_{x}\right|} $$ \overline{\left|\mathrm{\Delta }{\theta }_{y}\right|} $分别从0.1°增大到0.52°和0.56°。对于同一尘埃速度,$ \overline{\left|\mathrm{\Delta }Q\right|} $均比$ \overline{\left|\mathrm{\Delta }v\right|} $大,$ \overline{\left|\mathrm{\Delta }{\theta }_{x}\right|} $$ \overline{\left|\mathrm{\Delta }{\theta }_{y}\right|} $大小相当。需要注意的是,虽然各参数的平均反演精度较高,但对于单次反演结果仍可能出现较大的误差。例如,表1中尘埃速度的反演误差达到-7.63%。对于参数整体的偏离方向,尘埃电荷量和速度大小整体呈被低估的趋势(粉线),且当尘埃速度<40 km/s时,速度越大被低估得越多,电荷量和速度大小整体被低估的量不超过0.4%。当尘埃速度≥10 km/s时,尘埃轨迹倾角整体明显呈现被高估的趋势,θx整体被高估的量不超过0.05°,θy的不超过0.12°。

    • 反演误差随θy的变化如图10所示。

      图10可以看出,除了$ \overline{\left|\mathrm{\Delta }{\theta }_{y}\right|} $本身随着θy的增大而增大,轨迹倾角的变化对其他参数的反演误差没有显著影响。$ \overline{\left|\mathrm{\Delta }Q\right|} $在0.3%~0.4%之间,$ \overline{\left|\mathrm{\Delta }v\right|} $在0.15%~0.25%之间,$ \overline{\left|\mathrm{\Delta }{\theta }_{x}\right|} $$ \overline{\left|\mathrm{\Delta }{\theta }_{y}\right|} $均在0.12°左右。在反演结果偏离方向上,同样是尘埃电荷量和速度大小被整体低估,低估量都在0.15%以内,而轨迹倾角相对于“真”值没有明显的偏移方向。

    • 反演误差随Q的变化如图11所示。

      图11可以看出,尘埃电荷量的变化对各参数反演误差没有显著影响。$ \overline{\left|\mathrm{\Delta }Q\right|} $ 保持在0.4%左右,$ \overline{\left|\mathrm{\Delta }v\right|} $在0.2%左右,$ \overline{\left|\mathrm{\Delta }{\theta }_{x}\right|} $$ \overline{\left|\mathrm{\Delta }{\theta }_{y}\right|} $接近0.15°。在反演结果偏离方向上,尘埃电荷量和速度大小整体被低估0.1%左右,轨迹倾角相对于“真”值同样没有明显的偏移方向。

    • 反演误差随QNR的变化如图12所示。

      图12可以看出,尘埃电荷量和轨迹参数的反演误差随探测器QNR的增大而减小并逐渐趋于平缓,$ \overline{\left|\mathrm{\Delta }Q\right|} $小于0.4%,$ \overline{\left|\mathrm{\Delta }v\right|} $小于0.2%,$ \overline{\left|\mathrm{\Delta }{\theta }_{x}\right|} $$ \overline{\left|\mathrm{\Delta }{\theta }_{y}\right|} $均小于0.15°。在反演结果偏离方向上,同样是尘埃电荷量和速度大小被低估,但趋势随QNR的增大而减弱,低估量逐渐从0.1%附近变为接近于0,而轨迹倾角相对“真”值没有明显的偏移方向。

      通过上述分析可知,本文所提出的方法在其考虑的参数范围内,对尘埃电荷量的平均反演精度优于1.5%,速度大小的平均反演精度优于0.8%,轨迹倾角的平均反演精度优于0.6°。

    • 尘埃穿过DTS会诱导金属丝产生电信号,通过对观测信号和仿真信号进行拟合可提取尘埃的电荷量与轨迹信息。本文对DTS探测尘埃速度矢量进行了数值仿真,结果表明:

      1)尘埃电荷量和轨迹的反演误差随尘埃速度vz的增大而增大,随探测器QNR的增大而减小,基本不受轨迹倾角变化和尘埃带电量的影响。

      2)尘埃电荷量反演精度优于1.5%,速度大小反演精度优于0.8%,速度方向反演精度优于0.6°。

      本研究通过DTS探测获得了高精度的尘埃电荷量和轨迹信息,并对探测精度影响因素进行了仿真分析,可为DTS研制与科学探测提供参考。需要说明的是:由于反演过程需要7个独立参数,本文算法中通过随机计算103次来确定χ2的最小值很大可能是局部最小而不是全局最小;本文用于仿真的DTS结构为较优结构,并未针对其他结构类型(如不同的栅网平面间距、金属丝间距等)进行研究;本文对尘埃颗粒形状进行了简化(以球形代替),且用于计算探测器信号的公式与实际情况间存在一定偏差。后续将重点开展金属丝诱导电荷信号的理论和试验研究,依据更符合实际情况的信号修正本文的计算方法,进一步开展反演精度分析。

参考文献 (17)

目录

    /

    返回文章
    返回