CN105699973A - 利用等离子体线级联结构反演强扰动电离层参量的方法 - Google Patents
利用等离子体线级联结构反演强扰动电离层参量的方法 Download PDFInfo
- Publication number
- CN105699973A CN105699973A CN201610147359.0A CN201610147359A CN105699973A CN 105699973 A CN105699973 A CN 105699973A CN 201610147359 A CN201610147359 A CN 201610147359A CN 105699973 A CN105699973 A CN 105699973A
- Authority
- CN
- China
- Prior art keywords
- frequency
- ion
- omega
- cascade structure
- plasma
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Plasma Technology (AREA)
Abstract
一种利用等离子体线级联结构反演强扰动电离层参量的方法,包括步骤:(1)由非相干散射雷达观测数据,提取电离层加热实验引起的强扰动区域具有级联结构的异常非相干散射等离子体线功率谱;(2)对具有级联结构的异常非相干散射等离子体线功率谱中的相关峰值频差进行统计;(3)确定电离层加热过程中出现的低频波频率;(4)将步骤(2)中统计的频差与步骤(3)中计算得到的低频波频率进行对比,确定该级联结构是由电离层加热中参量衰变不稳定中哪种低频波产生的;(5)根据步骤(4)的结果,根据fpi的计算公式电子密度和电子温度。该方法能够很大程度提高反演精度,是获取强扰动电离层参量的一种新方法。
Description
技术领域
本发明涉及一种基于异常非相干散射等离子体线功率谱的级联结构反演强扰动电离层参量的方法,采用异常非相干散射等离子体线功率谱的级联结构和参量衰变不稳定性理论反演强扰动电离层参量的方法,该方法能够很大程度提高反演精度,是获取强扰动电离层参量的一种新方法。
背景技术
扰动电离层区域状态的探测是电离层加热异常特征实验和理论研究的重要部分之一。非相干散射雷达是电离层地面探测最为强大的手段之一。目前基于非相干散射雷达数据反演电离层参量使用的是平衡态非相干散射理论,但该方法并不能够很好地分析加热扰动电离层参量,尤其对于加热效应明显的强扰动区域。为了分析加热强扰动条件下的电离层状态建立了一些包括人工场向不规则体的非相干散射模型和碰撞等离子体非相干散射模型等理论,但在剧烈扰动区域,非相干散射功率谱的增强幅度多达数个量级的情况下,这些修正的非相干散射理论使用也会受限,使强扰动区域电离层参量的反演遇到巨大困难。
发明内容
本发明提供一种利用等离子体线级联结构反演强扰动电离层参量的方法,提取强扰动区域相应电离层高度的功率谱级联结构中相邻谱线频差,分析频差具有的规律性,根据参量不稳定性理论中低频振荡波特性获得电子温度和电子密度等参量。解决现有技术存在的在剧烈扰动区域,非相干散射功率谱的增强幅度多达数个量级的情况下,修正的非相干散射理论使用会受限,使强扰动区域电离层参量的反演遇到巨大困难得问题。
本发明采用的技术方案是:一种利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,包括以下步骤:
(1)由非相干散射雷达观测数据,提取电离层加热实验引起的强扰动区域具有级联结构的异常非相干散射等离子体线功率谱;
(2)对具有级联结构的异常非相干散射等离子体线功率谱中的相关峰值频差进行统计;
(3)确定电离层加热过程中出现的低频波频率,具体包括离子的等离子体频率fpi,离子的磁旋频率fci和离子声波频率fIA;
(4)将步骤(2)中统计的频差与步骤(3)中计算得到的低频波频率进行对比,确定该级联结构是由电离层加热中参量衰变不稳定中哪种低频波产生的;
(5)根据步骤(4)的结果,如果该级联结构频差与离子的等离子体频率一致,则根据fpi的计算公式可以获得电子密度;如果该级联结构的频差与离子声波频率fIA相等,则由离子声波的色散方程可以获得电子密度和电子温度。
本发明的优点是:该方法能够很大程度提高反演精度,是获取强扰动电离层参量的一种新方法。
具体实施方式
本发明一种基于非相干散射等离子体线功率谱的级联结构反演强扰动电离层参量的方法,包括以下步骤:
1)由非相干散射雷达观测数据,提取电离层加热实验引起的强扰动区域具有级联结构的异常非相干散射等离子体线功率谱。具体的实施方法是:
根据参量衰变不稳定性激发需要的电场阈值判断本次实验是否发生了参量衰变不稳定性。电场阈值计算公式为:
其中υe是电子碰撞频率,kb是玻尔兹曼常熟,ne是电子密度,Bmax为与电子温度和离子温度比值相关的系数,ε0为真空中的介电常数,ωpe是电子等离子体频率,Ti是离子温度。
通过上述判断,确定具有级联结构的异常非相干散射等离子体线功率谱是由参量衰变不稳定性引起的。具体的判断方法是:
将相关参数的值代入阈值公式计算获得EPDI,根据(ERP是有效辐射功率,z是电离层高度)计算加热电波在反射高度处的电场强度,讲两者对比,如果E(z)≥|EPDI|,则说明参量衰变不稳定性被激发,也就是说该级联结构是由参量衰变不稳定性引起的。
2)对具有级联结构的异常非相干散射等离子体线功率谱中的相关峰值频差进行统计。具体的实施方法是:
电离层加热引起的异常非相干散射等离子体线功率谱的级联结构会同时具有多个峰值。首先记录谱线中多个峰值对应的频率,计算第二个峰值对应频率与第一个峰值对应频率的差,第三个峰值对应频率与第二个峰值对应频率的差,直到第n个峰值对应频率与第n-1个峰值对应频率的差值。然后将统计的相邻峰值频差列表比较,通常这组数据是近似相等的。
3)确定电离层加热过程中出现的低频波频率,具体包括离子的等离子体频率fpi,离子的磁旋频率fci和离子声波频率fIA。
具体的实施方法是:根据本次实验中加热参数和电离层参数,求解电离层加热过程中可能出现的低频波的频率。三种低频频率的计算方法为:
离子的等离子体频率fpi和离子的磁旋频率fci分别利用下面两式求解:
fci==-eB/2πmi
离子声波频率fIA=ωIA/2π由色散方程求得。
其中,ni为离子数密度,mi为离子质量,e0为真空中的介电常数,e是单个电子的电荷数,B是磁感应强度,cIA为离子声速,角θ为波传播方向与地磁场的夹角,kIA是波数,ωIA是离子声波角频率,ωci是离子的磁旋角频率。
4)将步骤2)中统计的频差与步骤3)中计算得到的低频波频率进行对比,确定该级联结构是由电离层加热中参量衰变不稳定中哪种低频波产生的。
具体的实施方法是:通常离子的等离子体频率为几十Hz,离子声波频率典型值在几kHz,离子磁旋频率典型值在几十Hz,将步骤2)中统计的相邻峰值之间的频差与步骤3)中计算得到的三种低频波频率值对比,确定该频差的量级与哪种低频波频率对应,则认为该级联结构是由相应的低频振荡波产生。
5)根据步骤4)的分析,如果该级联结构频差与离子的等离子体频率一致,则根据fpi的计算公式可以获得电子密度;如果该级联结构的频差与离子声波频率fIA相等,则由离子声波的色散方程可以获得电子密度和电子温度。
补充:由离子等离子体频率fpi(步骤3中已补充公式)获得电子密度ne的公式为:
由离子声波频率的色散方程获得电子温度Te的近似公式为:
γe和γi为系数,可以近似取γe=1,γi=3,c是光速。
参量衰变不稳定性发生的物理过程中,需要满足的匹配条件为:
其中,ω1=ω0-ω2,ω0是加热泵波的角频率,ω2是加热电离层中本身存在一个平行于地磁场方向的低频电子密度扰动的角频率,k1是ω1对应的波数,k2是ω2对应的波数。b为与电子温度和离子温度相关的系数,vi为离子的动力学温度。λD是德拜长度。等离子体中电子的德拜长度可表示为:
类似的离子的德拜长度可以表示为:
ωpe是电子等离子体频率,具体表示为:
将获得的电子温度Te代入参量衰变不稳定性发生的匹配条件,计算得到电子等离子体频率ωpe。
利用电子密度ne与电子等离子体频率ωpe之间的关系,即:
可以求得电子密度。
Claims (6)
1.一种利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,包括以下步骤:
(1)由非相干散射雷达观测数据,提取电离层加热实验引起的强扰动区域具有级联结构的异常非相干散射等离子体线功率谱;
(2)对具有级联结构的异常非相干散射等离子体线功率谱中的相关峰值频差进行统计;
(3)确定电离层加热过程中出现的低频波频率,具体包括离子的等离子体频率fpi,离子的磁旋频率fci和离子声波频率fIA;
(4)将步骤(2)中统计的频差与步骤(3)中计算得到的低频波频率进行对比,确定该级联结构是由电离层加热中参量衰变不稳定中哪种低频波产生的;
(5)根据步骤(4)的结果,如果该级联结构频差与离子的等离子体频率一致,则根据fpi的计算公式可以获得电子密度;如果该级联结构的频差与离子声波频率fIA相等,则由离子声波的色散方程可以获得电子密度和电子温度。
2.根据权利要求1所述的利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,所述的步骤(1)的具体的方法是:
根据参量衰变不稳定性激发需要的电场阈值判断本次实验是否发生了参量衰变不稳定性,电场阈值计算公式为:
式中:EPDI是电场阈值,υe是电子碰撞频率,kb是玻尔兹曼常熟,ne是电子密度,Bmax为与电子温度和离子温度比值相关的系数,ε0为真空中的介电常数,ωpe是电子等离子体频率,Ti是离子温度;
通过上述判断,确定具有级联结构的异常非相干散射等离子体线功率谱是由参量衰变不稳定性引起的。
3.根据权利要求2所述的利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,所述的步骤(2)的具体的方法是:电离层加热引起的异常非相干散射等离子体线功率谱的级联结构会同时具有多个峰值;首先记录谱线中多个峰值对应的频率,计算第二个峰值对应频率与第一个峰值对应频率的差,第三个峰值对应频率与第二个峰值对应频率的差,直到第n个峰值对应频率与第n-1个峰值对应频率的差值;然后将统计的相邻峰值频差列表比较。
4.根据权利要求1所述的利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,所述的步骤(3)的具体的方法是:根据加热参数和电离层参数,求解电离层加热过程中可能出现的低频波的频率;三种低频频率的计算方法为:
离子的等离子体频率fpi和离子的磁旋频率fci分别利用下面两式求解:
fci==-eB/2πmi
离子声波频率fIA=ωIA/2π由色散方程求得;
式中,ni为离子数密度,mi为离子质量,ε0为真空中的介电常数,e是单个电子的电荷数,B是磁感应强度,cIA为离子声速,角θ为波传播方向与地磁场的夹角,kIA是波数,ωIA是离子声波角频率,ωci是离子的磁旋角频率。
5.根据权利要求3或4所述的利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,所述的步骤(4)的具体的方法是:通常离子的等离子体频率为几十Hz,离子声波频率典型值在几kHz,离子磁旋频率典型值在几十Hz,将步骤(2)中统计的相邻峰值之间的频差与步骤(3)中计算得到的三种低频波频率值对比,确定该频差的量级与哪种低频波频率对应,则认为该级联结构是由相应的低频振荡波产生。
6.根据权利要求4所述的利用等离子体线级联结构反演强扰动电离层参量的方法,其特征在于,所述的步骤(5)的具体的方法是:由步骤(3)中的离子等离子体频率fpi获得电子密度ne的公式为:
由离子声波频率的色散方程获得电子温度Te的近似公式为:
γe和γi为系数,可以近似取γe=1,γi=3,c是光速;
参量衰变不稳定性发生的物理过程中,需要满足的匹配条件为:
其中,ω1=ω0-ω2,ω0是加热泵波的角频率,ω2是加热电离层中本身存在一个平行于地磁场方向的低频电子密度扰动的角频率,k1是ω1对应的波数,k2是ω2对应的波数;b为与电子温度和离子温度相关的系数,vi为离子的动力学温度。λD是德拜长度。等离子体中电子的德拜长度可表示为:
类似的离子的德拜长度可以表示为:
ωpe是电子等离子体频率,具体表示为:
将获得的电子温度Te代入参量衰变不稳定性发生的匹配条件,计算得到电子等离子体频率ωpe;
利用电子密度ne与电子等离子体频率ωpe之间的关系,即:
求得电子密度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610147359.0A CN105699973B (zh) | 2016-03-16 | 2016-03-16 | 利用等离子体线级联结构反演强扰动电离层参量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610147359.0A CN105699973B (zh) | 2016-03-16 | 2016-03-16 | 利用等离子体线级联结构反演强扰动电离层参量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105699973A true CN105699973A (zh) | 2016-06-22 |
CN105699973B CN105699973B (zh) | 2018-04-06 |
Family
ID=56220603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610147359.0A Active CN105699973B (zh) | 2016-03-16 | 2016-03-16 | 利用等离子体线级联结构反演强扰动电离层参量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105699973B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107255802A (zh) * | 2017-07-21 | 2017-10-17 | 西安电子科技大学 | 基于离子线谱三峰结构反演强扰动电离层参量的方法 |
CN111580061A (zh) * | 2020-05-20 | 2020-08-25 | 中国科学院地质与地球物理研究所 | 基于clean算法的电离层电子密度反演方法 |
CN113686456A (zh) * | 2021-07-19 | 2021-11-23 | 中国人民解放军国防科技大学 | 基于rtg的反演电离层电子温度的方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008016914A2 (en) * | 2006-07-31 | 2008-02-07 | University Corporation For Atmospheric Research | Method and system for demodulation of open-loop gps radio occultation signals |
CN105022045A (zh) * | 2015-07-14 | 2015-11-04 | 华东交通大学 | 一种基于多源数据融合的三维电离层层析方法 |
CN105182304A (zh) * | 2015-06-04 | 2015-12-23 | 南昌大学 | 电离层非相干散射雷达非相干散射谱建模方法 |
-
2016
- 2016-03-16 CN CN201610147359.0A patent/CN105699973B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008016914A2 (en) * | 2006-07-31 | 2008-02-07 | University Corporation For Atmospheric Research | Method and system for demodulation of open-loop gps radio occultation signals |
CN105182304A (zh) * | 2015-06-04 | 2015-12-23 | 南昌大学 | 电离层非相干散射雷达非相干散射谱建模方法 |
CN105022045A (zh) * | 2015-07-14 | 2015-11-04 | 华东交通大学 | 一种基于多源数据融合的三维电离层层析方法 |
Non-Patent Citations (2)
Title |
---|
MUSONG CHENG ET AL.: "Observation of VHF incoherent scatter spectra disturbed by HF heating", 《JOURNAL OF ATMOSPHERIC AND SOLAR-TERRESTRIAL PHYSICS》 * |
程木松 等: "基于加热激发等离子体线级联结构的电子密度反演", 《电波科学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107255802A (zh) * | 2017-07-21 | 2017-10-17 | 西安电子科技大学 | 基于离子线谱三峰结构反演强扰动电离层参量的方法 |
CN111580061A (zh) * | 2020-05-20 | 2020-08-25 | 中国科学院地质与地球物理研究所 | 基于clean算法的电离层电子密度反演方法 |
CN113686456A (zh) * | 2021-07-19 | 2021-11-23 | 中国人民解放军国防科技大学 | 基于rtg的反演电离层电子温度的方法及系统 |
CN113686456B (zh) * | 2021-07-19 | 2023-09-01 | 中国人民解放军国防科技大学 | 基于rtg的反演电离层电子温度的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105699973B (zh) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Rotational shear effects on edge harmonic oscillations in DIII-D quiescent H-mode discharges | |
CN105699973A (zh) | 利用等离子体线级联结构反演强扰动电离层参量的方法 | |
CN104749156A (zh) | 拉曼光谱检测方法 | |
Lousto et al. | Exploring the outer limits of numerical relativity | |
CN112083509A (zh) | 时频电磁法中激发极化异常的检测方法 | |
Fu et al. | Electron temperature inversion by stimulated Brillouin scattering during electron gyroharmonic heating at EISCAT | |
CN102866429A (zh) | 一种地下水分布的确定方法 | |
Shi et al. | On the internal resonant modes in marching-on-in-time solution of the time domain electric field integral equation | |
CN107255802B (zh) | 基于离子线谱三峰结构反演强扰动电离层参量的方法 | |
Scales et al. | Observations and theory of ion gyro-harmonic structures in the stimulated radiation spectrum during second electron gyro-harmonic heating | |
Najmi et al. | Simulations of ionospheric turbulence produced by HF heating near the upper hybrid layer | |
Hikishima et al. | Evaluation of waveform data processing in wave-particle interaction analyzer | |
Fujii et al. | Frequency analysis of chaotic flow in transition to turbulence in taylor-Couette system with small aspect ratio | |
Mahmoudian et al. | Investigation of incoherent scatter radar spectra features with stimulated electromagnetic emissions at EISCAT | |
Mahmoudian et al. | Narrowband stimulated electromagnetic emissions (SEE) spectra: A new ionospheric diagnostic technique | |
Weichman et al. | Spatial variations of a passive tracer in a random wave field | |
Affolter et al. | Space charge frequency shifts of the cyclotron modes in multi-species ion plasmas | |
Aspden et al. | Inertia-gravity-wave generation: a geometric-optics approach | |
Kuzelev | Boundary conditions for the electromagnetic field in a waveguide with a thin-walled annular plasma in the context of application to plasma microwave electronics | |
CN113687149B (zh) | 基于rtg的电离层电子密度反演方法及系统 | |
Tanjia et al. | Plasmonic breathing modes in C 60 molecules: A quantum hydrodynamic approach | |
Toporkov | Doppler Properties of Sea Backscatter Simulated with Creamer or Choppy Wave Surface Models | |
Crameri et al. | A benchmark comparison of numerical topography-what are suitable sticky air parameters | |
Tsurugi et al. | Study on high frequency cut‐off characteristics of ground motions for intra‐slab earthquakes occurred in Southwest Japan | |
CN105717543A (zh) | 一种噪声压制方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |