CN113536211B - 一种用于east托卡马克磁信号快速提取的数值方法 - Google Patents
一种用于east托卡马克磁信号快速提取的数值方法 Download PDFInfo
- Publication number
- CN113536211B CN113536211B CN202110840131.0A CN202110840131A CN113536211B CN 113536211 B CN113536211 B CN 113536211B CN 202110840131 A CN202110840131 A CN 202110840131A CN 113536211 B CN113536211 B CN 113536211B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- signals
- response signal
- magnetic response
- weight coefficient
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000000605 extraction Methods 0.000 title claims abstract description 10
- 239000000523 sample Substances 0.000 claims abstract description 40
- 238000002474 experimental method Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000007547 defect Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000003491 array Methods 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000003556 assay Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/10—Nuclear fusion reactors
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Plasma & Fusion (AREA)
- Measuring Magnetic Variables (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
本发明提供了一种用于EAST托卡马克磁信号快速提取的数值方法。首先根据EAST托卡马克磁探针阵列在空间中的排布情况计算相应的权重系数并储存在控制系统中等待调用。在EAST实验放电过程中,通过磁探针阵列实时测量各处的磁响应信号,将相应位置的磁信号数值与权重系数相乘并通过计算公式提取特定环向分量的复数信号。最后将提取到的信号储存在系统中等待后续分析或传输给其他控制算法进行处理。本发明实现了快速提取EAST放电过程中磁响应信号的特定环向分量,该数值方法计算效率高,达到实时获取磁响应信号特定分量的目的,弥补了传统拟合方法因效率低导致只能通过实验后期处理信号的不足。
Description
技术领域
本发明涉及磁约束可控核聚变领域托卡马克装置放电实验数据的快速处理,特别涉及一种用于EAST托卡马克磁响应信号环向分量快速提取的数值方法。
背景技术
随着社会的进步,人类对能源的需求与日俱增,能源短缺问题逐渐凸显。目前人类所使用的能源主要为化石能源,化石能源储量有限且燃烧释放有害气体,因此亟需寻找化石能源的替代品。目前普遍认为核聚变能是解决能源问题的终极方案,其中磁约束法中的托卡马克装置是最有前景的实现方法。我国在此领域处于国际前列,其中中科院等离子体物理研究所的东方超环(EAST)托卡马克是我国重要装置之一。
为了监测以及实时控制EAST托卡马克放电过程中的等离子体状态,需要各种诊断设备来测量等离子体参数,其中磁探针阵列就是用来测量空间中磁响应信号的设备。通过对测量到的磁响应信号进行分析可以判断等离子体状态并采取相应的实时控制手段对等离子体进行控制来达到稳态运行的目的。某些实时控制算法需要磁响应信号的特定环向分量作为输入参数,但是传统的拟合手段效率较低,往往只能在实验后进行拟合分析,无法达到实时提取磁响应信号的特定环向分量的水平,这阻碍了许多实时控制算法的发展。本发明提出的方法具有较高的计算效率,足够实现实时提取磁响应信号特定分量的目的,同时保证了提取信号的准确性,弥补了传统拟合方法只能在实验完成后进行信号处理的不足,为发展各种以磁响应信号特定环向分量为输入的实时控制算法提供了基础,是一种高效且稳定可靠的数值方法。
发明内容
发明目的:实现实时提取EAST放电过程中磁响应信号的特定环向分量的目的,弥补传统拟合方法由于拟合效率较低而只能在实验后期对信号进行处理的不足,为发展多种以磁响应信号特定环向分量为输入参数的实时控制算法提供基础。
本发明的技术方案:
一种用于EAST托卡马克磁信号快速提取的数值方法,该数值方法计算效率高,可以达到实时获取磁响应信号特定分量的目的,同时可以保证提取信号的准确性,是一种高效且稳定可靠的数值方法,具体包括以下步骤:
步骤1:根据EAST托卡马克装置的磁探针阵列空间分布计算相应的权重系数W1、W2…Wn(n为环向磁探针数,通常为偶数);
步骤1.1:将n个磁探针分为两两一组 共/>组,并分别计算每组的权重系数;
步骤1.2:首先计算第一组磁探针所对应的权重系数W1和 上测量到的磁响应信号S1和/>满足如下方程组:
acosφ1-bsinφ1=S1
acosφn/2+1-bsinφn/2+1=Sn/2+1
其中a、b为待定系数(用于计算所需提取的磁响应信号),S为不同磁探针测量到的磁响应信号,φ为磁探针所在空间位置所对应的环向角大小。最终提取的磁响应信号环向分量为其中/>i为虚数符号;
求解如上方程组可得待定系数如下:
因此相应的权重系数W定义如下:
步骤1.3:根据步骤1.2中采用的方法计算剩余组磁探针所对应的权重系数;
步骤2:将步骤1中得到的权重系数W1、W2…Wn存储到控制系统中等待调用;
步骤3:在EAST托卡马克装置实验放电过程中,通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…Sn;
步骤4:调用步骤2中存储的权重系数并根据如下公式计算该时刻的磁响应信号环向分量Ct:
步骤5:将该磁响应信号环向分量Ct存储在系统中留待实验后分析或实时发送到需要该信号的实时控制算法作后续处理。
步骤6:重复进行步骤3-5即可实现实时提取磁响应信号特定环向分量的目的。
本发明的有益效果:本发明实现了实时提取EAST托卡马克装置放电过程中磁响应信号的特定环向分量的目的,弥补了传统拟合方法效率低的不足,计算效率高且同时保证信号准确性,是一种稳定高效的数值方法。
附图说明
图1本发明所适用的EAST托卡马克装置的横截面示意图。
图2本发明所适用的EAST托卡马克装置磁探针阵列环向空间分布示意图(俯视图)。
图3本发明快速提取的磁响应信号与传统拟合得到的磁响应信号的对比。
图4本发明实现实时提取EAST实验过程中磁响应信号特定环向分量的主要流程图。
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
EAST托卡马克装置具有如附图1所示截面位形,其中芯部位置为实验放电过程中包含等离子体的区域。等离子体外围装配多组磁探针阵列用于测量实验过程中产生的磁响应信号。每组磁探针阵列包含多个磁探针,磁探针沿着环向均匀分布,如附图2所示(以其中一组磁探针阵列为例)。环向均匀分布8个磁探针(Sen1、Sen2…Sen8),相邻磁探针环向相距磁探针测量的信号通常无法直接使用,需要对特定环向分量进行提取。本实施例通过计算各磁探针的权重系数来提高提取效率,可以实现快速提取特定环向分量的目的。首先将8个磁探针分成4组(Sen1-Sen5、Sen2-Sen6、Sen3-Sen7、Sen4-Sen8),根据下面具体实施步骤1中的方法计算相应的权重系数可得:
其中φ1-φ8分别为0、π、/>
在实验之前将权重系数存储到控制系统中等待调用。实验放电开始后,通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…S8,并采用上面计算得到的权重系数通过具体实施步骤4中的方法计算该时刻的磁响应信号环向分量Ct,然后将提取信号存储或发送到其他实时控制算法。重复测量、提取和存储发送等步骤多次直到实验放电结束即可达到实时获取磁响应信号特定环向分量的目的。与传统拟合方法(单次计算时间为秒量级)相比,该方法(单次计算时间为毫秒量级)极大提高了提取效率,并且保证了提取结果的准确性,提取结果与传统拟合结果十分接近,如附图3所示。
具体实施步骤如下:
步骤1:根据EAST托卡马克装置的磁探针阵列空间分布计算相应的权重系数W1、W2…Wn(n为环向磁探针数,通常为偶数);
步骤1.1:将n个磁探针分为两两一组 共/>组,并分别计算每组的权重系数;
步骤1.2:首先计算第一组磁探针所对应的权重系数W1和 上测量到的磁响应信号S1和/>满足如下方程组:
acosφ1-bsinφ1=S1
acosφn/2+1-bsinφn/2+1=Sn/2+1
其中a、b为待定系数(用于计算所需提取的磁响应信号),S为不同磁探针测量到的磁响应信号,φ为磁探针所在空间位置所对应的环向角大小。最终提取的磁响应信号环向分量为其中/>i为虚数符号;
求解如上方程组可得待定系数如下:
因此相应的权重系数W定义如下:
步骤1.3:根据步骤1.2中采用的方法计算剩余组磁探针所对应的权重系数;
步骤2:将步骤1中得到的权重系数W1、W2…Wn存储到控制系统中等待调用;
步骤3:在EAST托卡马克装置实验放电过程中,通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…Sn;
步骤4:调用步骤2中存储的权重系数并根据如下公式计算该时刻的磁响应信号环向分量Ct:
步骤5:将该磁响应信号环向分量Ct存储在系统中留待实验后分析或实时发送到需要该信号的实时控制算法作后续处理。
步骤6:重复进行步骤3-5即可实现实时提取磁响应信号特定环向分量的目的。
以上内容是结合优选技术方案对本发明所做的进一步详细说明,不能认定发明的具体实施仅限于这些说明。对本发明所属技术领域的普通技术人员来说,在不脱离本发明的构思的前提下,还可以做出简单的推演及替换,都应当视为本发明的保护范围。
Claims (1)
1.一种用于EAST托卡马克磁信号快速提取的数值方法,其特征在于,该数值方法实现了快速提取EAST托卡马克实验放电过程中磁探针阵列测量到的磁信号的特定分量,具体包括以下步骤:
步骤1:根据EAST托卡马克装置的磁探针阵列空间分布计算相应的权重系数W1、W2…Wn,n为环向磁探针数,为偶数;
步骤1.1:将n个磁探针分为两两一组共/>组,并分别计算每组的权重系数;
步骤1.2:首先计算第一组磁探针所对应的权重系数W1和 上测量到的磁响应信号S1和/>满足如下方程组:
acosφ1-bsinφ1=S1
acosφn/2+1-bsinφn/2+1=Sn/2+1
其中,a、b为待定系数,用于计算所需提取的磁响应信号;S为不同磁探针测量到的磁响应信号,φ为磁探针所在空间位置所对应的环向角大小;最终提取的磁响应信号环向分量为其中/>i为虚数符号;
求解如上方程组得待定系数如下:
因此相应的权重系数W定义如下:
步骤1.3:根据步骤1.2中采用的方法计算剩余组磁探针所对应的权重系数;
步骤2:将步骤1中得到的权重系数W1、W2…Wn存储到控制系统中,等待调用;
步骤3:在EAST托卡马克装置实验放电过程中,通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…Sn;
步骤4:调用步骤2中存储的权重系数并根据如下公式计算该时刻的磁响应信号环向分量Ct:
步骤5:将该磁响应信号环向分量Ct存储在系统中留待实验后分析或实时发送到需要该信号的实时控制算法作后续处理;
步骤6:重复进行步骤3-5,即实现实时提取磁响应信号特定环向分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110840131.0A CN113536211B (zh) | 2021-07-24 | 2021-07-24 | 一种用于east托卡马克磁信号快速提取的数值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110840131.0A CN113536211B (zh) | 2021-07-24 | 2021-07-24 | 一种用于east托卡马克磁信号快速提取的数值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113536211A CN113536211A (zh) | 2021-10-22 |
CN113536211B true CN113536211B (zh) | 2024-02-06 |
Family
ID=78088878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110840131.0A Active CN113536211B (zh) | 2021-07-24 | 2021-07-24 | 一种用于east托卡马克磁信号快速提取的数值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113536211B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116679863B (zh) * | 2023-07-28 | 2023-10-24 | 中国科学院合肥物质科学研究院 | 一种适用于east辐射偏滤器充气系统的数据处理系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5804965A (en) * | 1995-09-25 | 1998-09-08 | The United States Of America As Represented By The United States Department Of Energy | Method and apparatus for steady-state magnetic measurement of poloidal magnetic field near a tokamak plasma |
RU2191410C2 (ru) * | 2000-10-16 | 2002-10-20 | Российский научный центр "Курчатовский институт" | Система регистрации рентгеновского излучения в установках типа токамак |
CN106650951A (zh) * | 2016-12-29 | 2017-05-10 | 中国科学技术大学 | 基于时域全局相似度的多通道测量数据自动筛选方法 |
CN106991271A (zh) * | 2017-03-07 | 2017-07-28 | 中国科学院合肥物质科学研究院 | 一种适用于east偏滤器探针诊断数据处理的软件系统 |
CN107896412A (zh) * | 2017-09-27 | 2018-04-10 | 中国科学院合肥物质科学研究院 | 一种用于测量射频波特性的高频磁探针诊断系统 |
CN110232205A (zh) * | 2019-04-28 | 2019-09-13 | 大连理工大学 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
CN112930015A (zh) * | 2021-03-17 | 2021-06-08 | 核工业西南物理研究院 | 一种等离子体电磁特性诊断用电磁复合探针 |
CN112992385A (zh) * | 2021-05-06 | 2021-06-18 | 西南交通大学 | 一种准环对称仿星器磁场位形的设计方法 |
-
2021
- 2021-07-24 CN CN202110840131.0A patent/CN113536211B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5804965A (en) * | 1995-09-25 | 1998-09-08 | The United States Of America As Represented By The United States Department Of Energy | Method and apparatus for steady-state magnetic measurement of poloidal magnetic field near a tokamak plasma |
RU2191410C2 (ru) * | 2000-10-16 | 2002-10-20 | Российский научный центр "Курчатовский институт" | Система регистрации рентгеновского излучения в установках типа токамак |
CN106650951A (zh) * | 2016-12-29 | 2017-05-10 | 中国科学技术大学 | 基于时域全局相似度的多通道测量数据自动筛选方法 |
CN106991271A (zh) * | 2017-03-07 | 2017-07-28 | 中国科学院合肥物质科学研究院 | 一种适用于east偏滤器探针诊断数据处理的软件系统 |
CN107896412A (zh) * | 2017-09-27 | 2018-04-10 | 中国科学院合肥物质科学研究院 | 一种用于测量射频波特性的高频磁探针诊断系统 |
CN110232205A (zh) * | 2019-04-28 | 2019-09-13 | 大连理工大学 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
CN112930015A (zh) * | 2021-03-17 | 2021-06-08 | 核工业西南物理研究院 | 一种等离子体电磁特性诊断用电磁复合探针 |
CN112992385A (zh) * | 2021-05-06 | 2021-06-18 | 西南交通大学 | 一种准环对称仿星器磁场位形的设计方法 |
Non-Patent Citations (3)
Title |
---|
EAST上极向场线圈位置测量方法的研究与设计;李, 罗家融, 王华忠, 黄勤超;计算机测量与控制(08);全文 * |
多螺旋新经典撕裂模大规模模拟研究;王正汹等;《第十八届全国等离子体科学技术会议摘要集》;全文 * |
应用高速CCD诊断研究HT-7装置边界等离子体旋转;李凤娟;朱翔;杨建华;高翔;;核聚变与等离子体物理(03);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113536211A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102401217B1 (ko) | 타겟 음성 검출 방법 및 장치 | |
CN113536211B (zh) | 一种用于east托卡马克磁信号快速提取的数值方法 | |
Haskey et al. | Active spectroscopy measurements of the deuterium temperature, rotation, and density from the core to scrape off layer on the DIII-D tokamak | |
CN109934222A (zh) | 一种基于迁移学习的绝缘子串自爆识别方法 | |
WO2018041046A1 (zh) | 负荷开关电磁机构输出特性计算方法及容差分配方法 | |
CN110739091B (zh) | 一种核电厂次临界度直接测量系统及硬件平台和测量方法 | |
EP4009230A1 (en) | Current waveform fluctuation starting point identification method, electronic device and readable storage medium | |
CN106500991B (zh) | 基于自适应多尺度AVG-Hat变换的轴承故障信号特征提取方法 | |
CN107396392B (zh) | 无线接入点的射频指标测试方法、装置、设备及存储介质 | |
Sanyal et al. | Fracture behavior of thin-walled Zircaloy fuel clad tubes of Indian pressurized heavy water reactor | |
CN109669414B (zh) | 一种基于自相关特征分解的动态过程监测方法 | |
CN113963821B (zh) | 核反应堆泄漏及缺陷声发射信号的有效性判定方法和装置 | |
CN113163424B (zh) | 用于检测设备的nr小区pss搜索方法与检测设备 | |
CN110866346B (zh) | 一种用于获取干式空心电抗器固有振动特性的方法及系统 | |
CN115828088A (zh) | 基于深度学习的高压并联电抗器振动异常故障诊断方法 | |
CN113868031A (zh) | 一种硬盘性能一致性的测试方法、装置、系统 | |
CN111079267B (zh) | 一种基于微观残余应力突度的装载机底盘可靠性检测方法 | |
CN112651150B (zh) | 一种Al2O3/316L不锈钢阻氚系统表面多重裂纹行为的扩展有限元模拟分析方法 | |
CN109961788A (zh) | 风电中央监控系统的语音智能助手装置及语音控制系统 | |
CN110873672A (zh) | 一种wic镶嵌线铜比测量的方法 | |
CN110032758A (zh) | 计算电信号的能量的方法、装置和计算机存储介质 | |
CN111400869B (zh) | 一种堆芯中子通量时空演变预测方法、装置、介质及设备 | |
Li et al. | Research on online monitoring and fault diagnosis system based on multivariate empirical mode decomposition | |
CN114741658A (zh) | 一种本征正交基归一化方法、装置及存储介质 | |
CN117195743A (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 |