0 引言

电离层作为无线电信号的重要传播介质，其背景状态对信号的传输效果影响巨大.早在20世纪30年代，研究人员就提出了利用大功率高频电波与电离层相互作用产生人工光源照亮夜空的想法.20世纪70年代，美国、前苏联等相继建立大功率高频电离层加热装置，使人工改变电离层并影响无线电信号传输成为可能.最早，Utlaut等(1970)利用Platteville加热装置观测到了人工扩展F现象，并根据其特性推测在电离层中产生了人工不均匀体；随后，Thome和Blood(1974)利用位于新墨西哥州白沙导弹靶场的高频(High Frequency, HF)和甚高频(Very High Frequency, VHF)雷达，首次观测到了来自于人工不均匀体的雷达回波，证实了人工不均匀体沿地磁场分布的结构特性，因此被称为人工沿场不均匀体(Artificial Field Aligned Irregularities, AFAI)；在其进一步实验中，获得的雷达散射截面积(Radar Cross Section, RCS)在VHF频段最高达到108m2量级(Fialer, 1974)，基于此在多条链路上成功进行了不同速率报文、话音以及数据的远距离传输(Barry, 1974)，然而，此后公开文献中再无类似报道.21纪初，美国在阿拉斯加州建成了目前为止世界上功能最强的HAARP站，在其实验中，观测到了GPS接收信号扰动(Milikh et al., 2008)、人工扩展F(Kuo and Snyder, 2010; Kuo, 2015)等现象，进一步确认了AFAI的产生，但利用AFAI进行无线电信号传输的实验从未报道过，考虑到HAARP所处纬度高于Platteville，分析认为和AFAI散射特性随纬度的变化有关.

国外研究人员针对AFAI的理论研究大多集中在其产生机制方面，在AFAI的散射理论研究方面，Booker和Gordon(1950)曾经给出了无线电波在包含介电常数不均匀体的大气中的散射理论.虽然Rao和Thome(1974)曾经建立了人工散射体的RCS解析模型，但其电子密度扰动波数谱是利用Platteville站实验中的雷达数据拟合得到，无法给出AFAI散射特性随纬度的变化.

近年来，我国学者在电离层加热引起的电子温度和密度扰动(黄文耿等，2004；郝书吉等，2013)、气辉辐射(郝书吉等，2017)、VLF辐射与传播VLF传播(杨巨涛等，2017)以及人工沿场散射(闫玉波等，2008)等方面开展了大量理论研究.针对电离层加热产生人工不均匀体，国内学者主要是从电离层不稳定性的激发(赵正予和魏寒颖，2004，2005)和对电离层探测结果的影响(吕立斌等，2017)等方面开展研究，而关于AFAI散射特性方面的研究未见公开报道.分析中低纬度地区产生AFAI的散射特性及其和高纬度地区的差异，对我国开展相关应用研究具有重要的指导意义.

本文首先从各向异性介质的散射理论出发，通过推导电子密度扰动的波数谱表达式，建立AFAI的散射系数理论模型；然后利用国外实验数据对模型进行验证；在此基础上，重点利用模型仿真计算散射系数和方向特性随纬度的变化；最后讨论在不同纬度地区利用AFAI进行超短波远距离通信的优缺点，为在中低纬度开展AFAI试验及其散射特性应用研究提供理论基础.

1 理论模型

1.1 AFAI散射系数模型

在电离层加热实验中，采用适当的加热参数能够在电离层的反射高度(高频电波在电离层中反射的高度)附近产生大量沿地磁场分布的圆柱状电子密度不均匀体，其水平尺度通常可达百公里量级，被称为人工散射体.其中单个圆柱状电子密度不均匀体尺度通常在几米到几十米量级，由于其沿地磁场的结构特征，被称为人工沿场不均匀体(AFAI).

图1 AFAI几何模型 Fig.1 Geometry of AFAI

AFAI具有沿场结构特性，能够对入射的无线电波产生散射作用.为从理论上分析其散射特性，本文以AFAI为对象，建立其对无线电波的散射几何模型，如图1所示，假定从T站点发出功率为P0的无线电波照射人工沿场不均匀体v，根据Booker等的散射理论(Booker and Gordon, 1950)，在v中、和T距离为rT的点M处场强为

(1)

其中k和ζ分别为波矢量和特征阻抗.

若AFAI的电子密度扰动导致M点处的介电常数ε增大Δε，则距离M为rR的R处产生的极化势为

(2)

在v中选择参考点N，其与T和R的距离分别为和由于和远大于AFAI的线性尺寸，因此于是有

[4]张晖明、周岚岚、伍茜溪：《论税收手段对企业自主创新的激励作用——基于技术进步路径的视角》，《复旦大学学报》(社会科学版)2017年第6期。

(3)

其中

(4)

若是散射方向和电场E0方向的夹角，则R处的散射场为

(ls-li)2+(ms-mi)2=

(5)

散射系数可用单位体积的散射功率和入射功率之比表示，令(li,mi,ni)和(ls,ms,ns)分别为入射方向和散射方向的单位矢量，则散射系数可表示为

×P{k(ls-li),k(ms-mi),k(ns-ni)}，

(6)

其傅里叶变换

将(11)、(12)代入(10)中，可得到AFAI的波数谱

(7)

其中是电子密度扰动的均方差，用于衡量电子密度扰动强度.AFAI中的电子密度扰动近似服从高斯分布，可采用高斯自相关函数表示，在直角坐标系中形式为

1309 Long non-coding RNA-8439 up-regulates expression of pluripotency factor nanog in hepatoma carcinoma cells

(8)

其中P(l,m,n)是的自相关函数，即电子密度扰动波数谱.对于电离介质有

(9)

对于圆柱状AFAI，有a=b=T,c=L.其中，T是垂直对称轴的相关距离，L是沿对称轴的相关距离.则有

P{k(ls-li),k(ms-mi),k(ns-ni)}

[T2(ls-li)2

+T2(ms-mi)2+L2(ns-ni)2]}.

(10)

P{k(ls-li),k(ms-mi),k(ns-ni)}

BDBO主要针对畸形的顶点或者旋转轴的中心位于椎间盘区域的患者，或者重度矢状位畸形需要截骨的范围超过PSO所能达到的最大截骨角度时也可应用BDBO。BDBO不适用于椎间融合及脊柱前柱钙化或AS患者。

（2）从中国的出口商品结构看，改革开放40年来，中国出口商品的结构由原料和工业半成品为主逐渐转为机电设备、通信技术设备等科技产品为主，而此类产品具有低成本高利润的特点，无形中对美国同类产品形成竞争。因此，美国的厂商想方设法利用“337条款”夺回失去的市场份额。

(11)

和

E=k2sin

(12)

“唾弃你？不，你该珍惜自己的骨骼，不要憎恨自己，你的身躯将永远圣洁如莲。停尸房根本没有闹鬼，我会保护你免受他的骚扰，和你一起看守停尸房。”

在图1中以M为原点建立三维直角坐标系，z轴沿地磁场方向，y轴垂直向外，x轴垂直于地磁场方向.假设入射波矢量和x,y,z轴的夹角分别为αix,αiy,αiz，散射波矢量和x,y,z的夹角分别为αsx,αsy,αsz，则(li,mi,ni)和(ls,ms,ns)分别可表示为(cosαix,cosαiy,cosαiz)和(cosαsx,cosαsy,cosαsz).若散射角为φs，则有

(13)

其中k=2πf/c，f为信号频率，c为光速.将(7)和(13)代入(6)中，当水平极化波入射时，=90°，得到AFAI的散射系数

(14)

根据式(14)可以看出，影响AFAI散射系数的因素包括加热导致的电子密度扰动强度δn、背景等离子体波长λN、空间相关距离T和L、信号频率f、入射波和地磁场方向的夹角αiz、散射波和地磁场方向的夹角αsz以及散射角φs.其中，δn理论上可利用电离层F层加热模型进行计算，或是基于非相干散射雷达测量数据反演得到；λN等于大功率高频加热波的波长；T和L取决于AFAI的空间尺度，通常利用实验观测得到；αiz、αsz和φs取决于收发位置和AFAI的空间几何关系.根据这些参数，即可对不同纬度不同条件下产生AFAI的散射系数进行计算，并对其散射强度特性进行分析.

1.2 AFAI散射约束条件

由于AFAI沿地磁场的结构特征，能够将入射的无线电波沿以地磁场方向为轴的锥面进行散射，这种散射形象的被称为人工沿场散射(Artificial Field Aligned Scattering，AFAS)，如图2所示.由于AFAS具有很强的方向性，只有当入射波和散射波与地磁场方向的夹角相等时，才能获得有效的散射，因此，基于AFAI的无线电信号传输必须满足两个条件：第一，从地面发射的无线电信号能够到达AFAI；第二，从AFAI散射的无线电信号能够到达地面.

图2 AFAS及其约束条件 Fig.2 AFAS and the confine condition

地面上与AFAI的连线和地面相切的位置，即为发射站的最远位置，此时发射站到AFAI的距离为Dmax.令发射站到AFAI的距离为DTS，则基于AFAI进行无线电信号传输的第一个必要条件可表示为

染色体研究已经为蔬菜遗传育种提供了理论依据和实践指导。根据染色体形态进行染色体核型分析，还存在许多局限性，如对形态相似、大小相近、染色体较小等不好区分。因此，在蔬菜植物的核型分析中可以应用荧光原位杂交技术进行基因定位、染色体变异及种间与种内杂交研究，还可以开展多倍体研究来指导蔬菜育种。由于多倍体的生长优势、适应能力、群体产量及抗病力等各方面都明显比二倍体具有优势，多倍体研究对蔬菜品种选育具有重要意义。

DTS≤Dmax，

(15)

若切线和水平面的夹角为θ，地磁倾角为I，入射波和地磁场方向的夹角为αiz，则基于AFAI进行无线电信号传输的第二个必要条件可表示为

αiz<180-(θ+I).

(16)

(15)和(16)式共同构成了基于AFAI的人工沿场散射约束条件，满足约束条件的地面范围即为AFAS地面覆盖区.利用(15)和(16)式可对基于AFAI散射的方向性进行分析.

2 数值模拟

根据第1节建立的AFAI散射系数模型和AFAS约束条件，本节首先利用Platteville站的实验结果对模型进行验证，重点从散射强度和方向性两个方面，对AFAI的散射特性进行理论分析和研究.

2.1 模型验证

以Platteville站的典型加热实验(Minkoff, 1974)为例，其主要实验参数为：有效辐射功率100 MW，加热波波长54 m.实验测量的电子密度扰动强度δn约为1%，利用VHF雷达观测的人工散射体半宽度、半厚度和高度分别约为50 km、7.5 km和250 km，T和L的典型值分别为0.5 m和1 km.利用本文建立的散射系数模型对实验结果进行计算，20～80 MHz的散射系数约为-19.6～-13.9 dB.根据人工散射体的尺度，对计算得到的散射系数进行体积分，获得人工散射体RCS约为74.3～80 dBm2，和试验结果基本一致，验证了模型的正确性，如图3所示.

图3 Platteville站RCS理论算与实验结果对比 Fig.3 Comparision of RCS for Platteville

2.2 散射系数模拟结果

不同纬度地区地磁场存在差异，根据加热原理，在同等条件下加热产生的电子密度扰动强度不同，导致AFAI的散射强度(散射系数)存在差异.以110°E经线上10°N—80°N(对应的地磁倾角从6°变化到87°)为例；计算中背景电离层参数由国际参考电离层模型IRI2016给出，如图4所示，其中和NO+分别为电子、NO+离子离子和O+离子的密度，Te和Ti分别为电子和离子的温度；不同纬度地区的地磁场参数由国际参考地磁场模型IGRF-11给出；加热频率fheat=0.95f0，f0为电离层临界频率；有效辐射功率选取目前加热站常用的功率等级，即100～1000 MW；电子密度扰动强度利用电离层F层加热模型进行计算；T取0.5 m，L取1 km，φs取180°.

利用上述计算参数，仿真计算了频率分别为20 MHz、40 MHz、60 MHz和80 MHz时，AFAI的散射系数随纬度和有效辐射功率的变化如图5所示.

本报讯10月25日，史丹利发布2018年三季报，公司2018年1-9月实现营业收入42.91亿元，同比增长9.03%；化学制品行业已披露三季报个股的平均营业收入增长率为4.70%；归属于上市公司股东的净利润1.97亿元，同比下降17.87%，化学制品行业已披露三季报个股的平均净利润增长率为41.17%；公司每股收益为0.17元。

2017年10月底在各密度林分中选取3个具有代表性的样地，每样地测20株样本，用轮尺测量胸径、勃鲁莱测高器测量树高、卷尺测量冠幅，优势木高度的测定采取6株最高木法[12]，通过胸高形数法计算立木单株材积得到蓄积量。

图5计算结果表明，散射系数随纬度的升高而增大，从10°N到70°N，散射系数增大了约15 dB；散射系数随有效辐射功率的增大而增大，有效辐射功率从100 MW增大到1000 MW，散射系数增大了12 dB，但变化率逐渐减小；散射系数随信号频率的增大而减小，但在20～80 MHz频段，仍能产生最高约108m2量级的散射截面积，这为超短波信号的远距离传输提供了一种全新的途径.

图6进一步给出了HAARP站(地磁倾角75.6°)、Platteville站(地磁倾角67.7°)、武汉(地磁倾角44.2°)和广州(地磁倾角31.4°)四个典型站点在信号频率为40 MHz时的散射系数计算结果，可以看出，HAARP站的散射系数略高于Platteville站，比中低纬度地区的武汉大5 dB，比广州大10 dB.分析认为，由于高电离层中沿地磁场方向的扩散系数和热传导系数要远远大于横跨地磁场方向，电子密度扰动强度主要受沿地磁场方向的加热能量影响，高纬度地区容易产生电子密度扰动，进而能够获得比中低纬度大的散射系数.但在武汉和广州地区，若分别采用250 MW和650 MW的加热有效辐射功率，理论上可产生和Platteville站相当的散射系数，进而实现同样的无线电传输效果.

图4 背景电离层参数 (a) 电子和离子密度; (b) 电子和离子温度. Fig.4 Background parameters (a) Density; (b) Temperature.

图5 散射系数随地磁倾角变化 (a) 20 MHz; (b) 40 MHz; (c) 60 MHz; (d) 80 MHz. Fig.5 Changes of scattering coefficient with geomagnetic inclination

2.3 散射覆盖区模拟结果

同样由于人工沿场散射(AFAS)对地磁场方向的敏感性，不同纬度地区由于地磁场的差异，导致AFAS的散射覆盖区不同.利用(15)和(16)所示的约束条件，计算了AFAI高度从150 km变化到300 km时Platteville、HAARP、武汉以及广州的地面散射覆盖区，图7给出了AFAI高度为200 km和300 km的典型仿真结果，其中(a)—(d)为AFAI高度200 km时的散射覆盖区，(e)—(h)为AFAI高度300 km时的散射覆盖区.从图中可以看出：基于AFAI的散射覆盖区呈扇形区域，关于磁子午面对称并且大部分位于加热站以南；随着AFAI高度的升高，扇形区域的半径变大，张角减小，散射覆盖区整体向南偏移，纬度越高偏移越明显.值得注意的是，散射覆盖区随加热站点纬度的升高而变小甚至消失.对于高纬度站点，由于其较大的地磁倾角，使得能够满足式(13)和(14)条件的地面范围很小，例如HAARP站在AFAI高度200 km时散射覆盖区很小，而在300 km时不存在散射覆盖区，散射信号无法被地面接收.而对于中低纬度地区的武汉和广州，由于地磁倾角较小，具有较大的地面散射覆盖区，东西向范围可达3000 km，南北向可达1500 km，在基于电离层加热的AFAI研究及应用中具有独特的地理条件优势.

图6 典型站点的散射系数 Fig.6 Scattering coefficient for several typical stations

图7 不同站点的散射覆盖区随高度的变化 (a)—(d) 200 km； (e)—(h) 300 km. Fig.7 Change of scattering cover zone with height of different location

3 结论与讨论

本文建立了基于电离层加热的AFAI散射系数理论模型，给出了AFAS约束条件，在此基础上，仿真分析了AFAI散射系数和方向特性随纬度的变化，主要结论如下：

本发明公开了一种制备过渡金属硫化物二维纳米片层分散液的方法，将质量比为0.01～100∶1的物质A和过渡金属硫化物二维纳米片层的分散液加入到反应容器内超声，温度0～80 ℃，超声功率1～750 W，超声反应0.5～72 h；超声结束进行离心洗涤，得到分散在不同体系中的功能化产品的分散液；其中，物质A为HS(CH2)nR，R为COOH、NH2、OH或SH，n=2、3、4。该方法操作简单，操作条件温和，所制备得到的功能化产品在不同的有机溶液中都有较好的分散性，可以广泛用于新型复合材料的制备、新型光学及导电材料等领域。

(1) 对于同样背景条件，散射系数随有效辐射功率增大而增大，有效辐射功率从100 MW增大到1000 MW，散射系数增大约12 dB；

(2) 在同等有效辐射功率和背景条件下，AFAI的散射系数随着纬度的升高而增大，但变化率逐渐减小.根据仿真计算结果，纬度从10°～70°，散射系数最大增加约15 dB.典型的，Platteville和HAARP站相当，比武汉和广州大5～10 dB，根据结论1)，通过适当增大加热有效辐射功率，我国中低纬度地区可以产生和高纬度地区相当的散射系数；

己巳，诏曰：“朕闻至治之世，元、凯共朝，不为朋党，君明臣哲，垂荣无极，何其德之盛也。……至于属文之人，类亡体要，诋斥前圣，放肆异言，以讪上为能，以行怪为美。自今委中书、门下、御史台采察以闻。”［1］3718

(3) AFAI的散射覆盖区为关于磁子午面对称的扇形区域，主要位于加热站以南并且随着AFAI高度的升高逐渐向南移动；

(4) 受地磁场和人工沿场散射约束条件限制，散射覆盖区随纬度的升高而变小甚至消失.对于高纬度地区的HAARP站，散射覆盖区很小，在AFAI高度300 km时完全消失；而对于武汉和广州等我国中低纬度地区，散射覆盖区东西范围可达3000 km，南北向范围可达1500 km，具有明显的地理条件优势.

对于高纬度地区，尽管在电离层加热过程中能够产生较大的散射系数，但由于其地磁倾角较大，对于HAARP站，散射覆盖区很小，难以开展相关实验；而在Platteville站，虽然散射覆盖区相对较大，但相当一部分区域位于其国土之外.这可能是Platteville站没有深入开展无线电信号传输实验的原因.而我国大部分地区位于中低纬度，具有得天独厚的地理条件优势，通过采用较大的有效辐射功率或高效的AFAI激发方式，可以产生较大的散射系数，特别是在超短波频段，可获得约108m2的散射截面积，并且由于AFAI位于高度较高电离层F层，完全可用于超短波频段的超远距离传输，形成一种新型的超短波通信方式，本文研究结果可为中低纬度相关实验的开展提供理论支持.

①针对项目质量与成本控制等责任和奖罚措施，作业队应予以明确且在执行过程中做到奖罚分明，以此提升工作人员积极性；②对于项目成本控制措施，作业队操作人员应在现场技术人员领导下对其制定，同时利用培训的方式强化作业人员成本控制理念，并以项目实况为依据，针对成本核算科学建立岗位责任制，以此对作业人员在核算工作中的地位、作用、责任以及考核与奖励办法明确规定。此外，定期检查成本控制情况，通过阶段性实际成本与计划成本的比较，找出差距及分析原因，并以此采取措施实施纠偏处理，实现项目成本的有效控制。

监控系统采用SQL网络型数据库并提供了设备实时数据的第三方接口，提供手机APP应用，可以与用户实际管理体制进行配合管理，真正做到无人值守并融入用户的管理方式。将以前不注重但却对电力安全运营拥有相当影响的端子箱，纳入管理系统中来，进一步提高电力运维的安全性。监控平台界面显示采用全中文界面，图形化设计，支持电子地图功能。界面的结构、层次清晰明了，能够实时直观地显示设备的运行数据和运行状态，场景仿真。

此外，(14)式中的空间相关距离T和L取决于AFAI的空间尺度，由于包含在指数项内，不同的取值将对计算结果产生较大影响.以典型的东西向链路为例，不同频率的散射系数随T和L的相对变化如图8所示.可以看出，对于不同频率，均存在T0和L0，能够使散射系数最大，并且随着频率的增大，T0和L0的值增大.由此可知，通过控制加热参数生成特定结构的AFAI，可有效提升其散射能力.

然而，国外试验结果表明，电离层加热中产生的不规则体具有较宽的谱结构，横向尺度从cm到km量级，和背景电离层以及加热参数有密切的联系，但限于尚无成熟的不规则体生成理论，本文采用了电离层加热研究中常用的欧姆加热模型来计算电子密度扰动强度，同时利用国外典型试验结论进行分析，导致模型仍有一定的局限性.后续研究中将对不规则体的谱结构和背景电离层以及加热参数的关联性进行分析，研究特定结构不规则体的生成控制方法.

3.2.2 内环境保持相对稳定的意义 内环境是体内细胞生存的直接环境，细胞与内环境、内环境与外界环境之间不断进行物质交换。因此，细胞的代谢活动以及外界环境的不断变化，必然会影响内环境的理化性质。那么内环境的理化性质会变吗？再次出示血液报告单，并追问： 内环境的理化性质会剧烈的变化吗？学生自学讨论后，教师总结： 动物机体内部环境相对稳定的状态叫做稳态，稳态是一种可变却又相对稳定的状态。再次提问： 稳态会被破坏吗？请学生分别举些稳态的例子和稳态失调的例子，并说明稳态被破坏的危害。

图8 散射系数随T和L的变化 (a) 20 MHz; (b) 40 MHz; (c) 60 MHz; (d) 80 MHz. Fig.5 Changes of scattering coefficient with T and L

References

Barry G H. 1974. HF-VHF communications experiment using man-made field-aligned ionospheric scatterers. Radio Science, 9(11): 1025-1032.

Booker H G, Gordon W E. 1950. A theory of radio scattering in the troposphere.Proc. IRE, 38(4): 401-412.

Booker H G. 1956. A theory of scattering by nonisotropic irregularities with application to radar reflections from the aurora. Journal of Atmospheric and Terrestrial Physics, 8(4-5): 204-221.

Fialer P A. 1974. Field-aligned scattering from a heated region of the ionosphere-Observations at HF and VHF.Radio Science, 9(11): 923-940.

Hao S J, Li Q L, Yang J T, et al. 2013. Theory and numerical modeling of under-dense heating with powerful X-mode pump waves. Chinese Journal of Geophysics (in Chinese), 56(8): 2503-2510, doi: 10.6038/cjg20130801.

Hao S J, Ma G L, Yang J T, et al. 2017. Research on electron energy distribution with airglow measurements during ionospheric modification experiments.Chinese Journal of Geophysics (in Chinese), 60(3): 879-886, doi: 10.6038/cjg20170302.

Huang W G, Gu S F, Gong J C. 2004. Ionospheric heating by powerful high-frequency radio waves. Chinese Journal of Radio Science (in Chinese), 19(3): 296-301.

Kuo S,Snyder A. 2010. Observation of artificial spread-F and large region ionization enhancement in an HF heating experiment at HAARP. Geophysical Research Letters, 37: L07101, doi: 10.1029/2010GL042656.

Kuo S P. 2015. Ionospheric modifications in high frequency heating experiments.Physics of Plasmas, 22: 012901, doi: 10.1063/1.4905519.

Lv L B, Li Q L, Hao S J, et al. 2017. Theoretical analysis of effects on high frequency vertical sounding by artificial field-aligned irregularities.Acta Phys. Sin., 66(5): 059401, doi: 10.7498/aps.66.059401.

Milikh G, Gurevich A, Zybin K, et al. 2008. Perturbations of GPS signals by the ionospheric irregularities generated due to HF-heating at triple of electron gyrofrequency.Geophysical Research Letters, 35: L22102, doi: 10.1029/2008GL035527.

Minkoff J. 1974. Radio frequency scattering from a heated ionospheric volume, 3, Cross section calculations.Radio Science, 9(11): 997-1004.

Rao P B, Thome G D. 1974. A model for RF scattering from field-aligned heater-induced irregularities.Radio Science, 9(11): 987-996.

Thome G D, Blood D W. 1974. First observations of RF backscatter from field-aligned irregularities produced by ionospheric heating.Radio Science, 9(11): 917-921.

Utlaut W F, Violette E J, Paul A K. 1970. Some ionosonde observations of ionosphere modification by very high power, high frequency ground-based transmission. Journal of Geophysical Research, 75(31): 6429-6435.

doi: 10.6038/cjg20170201.

Yan Y B, Lv L B, Li Q L. 2008. Design of communication circuits based on AFAS.Chinese Journal of Radio Science (in Chinese), 23(3): 520-524, 538.

Yang J T, Li Q L, Wang J G, et al. 2017. Theoretical analysis of effects of artificial ionospheric disturbance on VLF propagation. Chinese Journal of Geophysics (in Chinese), 60(2): 463-469,

Zhao Z Y, Wei H Y. 2004. General dispersion relation for the three-wave process of parametric excitation (Ⅰ) The preferred excitation of an electron Langmuir wave and an ion acoustic wave. Chinese Journal of Space Science (in Chinese), 24(6): 441-447.

Zhao Z Y, Wei H Y. 2005. General dispersion relation for the three-wave process of parametric excitation (Ⅱ) The Basic Formula, the Field Threshold and the Growing Rate. Chinese Journal of Space Science (in Chinese), 25(1): 17-22.

附中文参考文献

郝书吉, 李清亮, 杨巨涛等. 2013. 大功率高频X波欠密加热电离层的理论与数值模拟. 地球物理学报, 56(8): 2503-2510, doi: 10.6038/cjg20130801.

郝书吉, 马广林, 杨巨涛等. 2017. 利用气辉测量结果研究高频电波加热电离层实验中电子能量分布. 地球物理学报, 60(3): 879-886, doi: 10.6038/cjg20170302.

黄文耿, 古士芬, 龚建村. 2004. 大功率高频无线电波加热电离层. 电波科学学报, 19(3): 296-301.

吕立斌, 李清亮, 郝书吉等. 2017. 人工沿场不均匀体对短波垂直探测影响的理论分析. 物理学报, 66(5): 059401, doi: 10.7498/aps.66.059401.

闫玉波, 吕立斌, 李清亮. 2008. 基于AFAS的通信链路设计. 电波科学学报, 23(3): 520-524, 538.

杨巨涛, 李清亮, 王建国等. 2017. 基于人工电离层扰动对地-电离层波导传播的VLF波的影响. 地球物理学报, 60(2): 463-469, doi: 10.6038/cjg20170201.

赵正予, 魏寒颖. 2004. 参量激励过程中三波耦合的一般色散关系(Ⅰ)最容易激励参量不稳定性的频率和波矢条件. 空间科学学报, 24(6): 441-447.

赵正予, 魏寒颖. 2005. 参量激励过程中三波耦合的一般色散关系(Ⅱ)一般色散关系以及泵场阈值和增长率. 空间科学学报, 25(1): 17-22.