CN113916717A - 一种基于低轨航天器掩星的平流层大气密度反演方法 - Google Patents
一种基于低轨航天器掩星的平流层大气密度反演方法 Download PDFInfo
- Publication number
- CN113916717A CN113916717A CN202111183403.0A CN202111183403A CN113916717A CN 113916717 A CN113916717 A CN 113916717A CN 202111183403 A CN202111183403 A CN 202111183403A CN 113916717 A CN113916717 A CN 113916717A
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- atmospheric density
- occultation
- inversion
- refraction
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
-
- 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
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种基于低轨航天器掩星的平流层大气密度反演方法,属于遥感技术领域。首先给出了平流层大气密度反演模型,此模型为本发明实施的基础;其次结合星光折射的几何关系以及航天器轨道的确定,给出了反演大气密度点的三维坐标以及掩星开始时航天器位置确定的方法。利用本发明可以获得全球三维、准实时的大气密度模型。弥补了直接测量以及地基遥感的空间覆盖率低的缺陷,并且相对无线电掩星方法成本更低,结构更加简单相对易于实现。
Description
技术领域
本发明涉及一种基于低轨航天器掩星的平流层大气密度反演方法,属于遥感技术领域。
背景技术
大气模型是对大气中的各种参数的时空分布的描述,大气参数的确定对于气象、天文以及航天等领域都有着重要的影响,所以很早就已经得到了广泛的研究。目前的大气探测主要包括直接量测、地基遥感以及天基遥感三类方法。
直接测量是指利用一系列仪器对载体所在位置大气参数的测量,例如温度计、测压计以及风速表等。在对空中气象进行测量时一般需要搭载探空气球、火箭以及卫星等航空航天飞行器。其中探空气球可以用来获取低层大气的垂直数据,其遇到的主要问题就是高度范围受限[1]和水平干扰[2]。卫星大气阻力反演则可以用于测量高层大气密度数据[3],但是卫星阻力反演的时空分辨率以及精度均存在一定的局限。
与直接量测相比,地基遥感在高度探测范围上有了质的突破。一般而言,地基遥感的可靠性、成本都很不错,但是这类方式同样会受到高山、沙漠、海洋等特殊地域的限制。根据中频、甚高频、流星、激光等不同雷达所选择的频率范围,其空间分辨率从1公里至数公里不等,时间分辨率从10分钟到60分钟不等。
天基遥感,即气象卫星出现的直接原因是气象的数值预报工作对于增加气象资料时间和空间分布密度的需求。直接观测以及地基遥感虽然可以获得比较精确的资料,但是其覆盖范围相对于全球而言仍只占很小一部分。美国发射的NOAA(National Oceanic andAtmospheric Administration)卫星通过不同垫面在可见光和红外波段下的特性分析温度、臭氧含量等成分,其一个瞬时观测区可以占据地球表面面积的2%。一个极轨卫星可以在半天时间内得到全球主要地区的气象资料。赤道上的四颗静止卫星,理论上可以覆盖全部赤道地区[4]。
无线电掩星是一种新颖的大气观测方式,不同于传统的气象卫星,无线电掩星不是向星下点接受信号,而是接受其它卫星的电磁信号,这一特点赋予了它更灵活的探测范围以及更多的测量信息。最早的无线电反演思想源于60年代,G.Fjeldbo等利用该方法分析了火星低层大气的折射率变化情况[5]。直至90年代以后,无线电掩星才开始被用于地球的大气探测[6]。目前,无线电掩星已经发展出几何光学反演、菲涅尔衍射反演等多种反演手段,几何光学反演对于星光折射反演具备明显的借鉴价值[7]。无线电掩星的问题在于信号源的成本,如果要实现大范围的大气观测,那么卫星的轨道设计、卫星间的协调、不同大气成分对电磁波信号的影响都是需要解决的困难。
参考文献:
[1]刘志欣,申延美,王蕊.提高探空气球探测高度的方法[J].黑龙江气象,2004,5(2):38-39.
[2]刘红亚,薛纪善,沈桐立,等.探空气球漂移及其对数值预报影响的研究[J].应用气象学报,2005,16(04):105-113.
[3]李若曦.利用低轨卫星数据反演热层大气密度[D]:中国科学技术大学,2017.
[4]G.Davis.History of the NOAA satellite program[J].Journal ofApplied Remote Sensing,2007,1(1):341-353.
[5]G.Fjeldbo,V.R.Eshleman.The atmosphere of mars analyzed by integralinversion ofthe Mariner IV occultation data[J].Planetary Space Science,1968,16(8):1035-1059.
[6]E.R.Kursinski,G.A.Hajj,W.I.Bertiger,et al.Initial Results of RadioOccultation Observations of Earth's Atmosphere Using the Global PositioningSystem[J].Science,1996,271(5252):1107-1110.
[7]王鑫,吕达仁.GPS无线电掩星技术反演大气参数方法对比[J].地球物理学报,2007,02):346-353.
发明内容
本发明技术解决问题:克服现有技术的不足,提供一种基于低轨航天器掩星的平流层大气密度反演方法,基于低轨航天器利用星光折射信息反演大气参数,获得全球平流层三维、准实时的大气模型参数,同时弥补了直接测量以及地基遥感的空间覆盖率低的缺陷,并且相对无线电掩星方法成本更低,结构更加简单相对易于实现。
本发明技术解决方案:一种基于低轨航天器掩星的平流层大气密度反演方法,基于星敏感器利用星光折射信息的大气密度反演方法,包括以下步骤:
(1)计算并记录掩星开始时刻航天器的位置;
(2)从航天器到达步骤(1)确定的位置开始,对星敏感器拍摄星图进行星图处理,记录折射前星光矢量方向和星光折射角,同时记录航天器导航***输出的位置信息,持续记录数据直到掩星过程结束;
(3)利用步骤(2)记录的航天器位置、折射前星光矢量方向以及折射角数据,根据几何关系计算视高度;
(4)利用步骤(2)记录的折射角数据以及步骤(3)计算的视高度数据,利用基于星光折射的大气密度反演公式计算每组数据对应的大气密度;
(5)利用步骤(2)记录的航天器位置、折射前星光矢量方向以及折射角数据,根据几何关系计算反演大气密度点的三维坐标;
所述步骤(4)中,基于星光折射的大气密度反演公式为:
式中,ρ代表反演大气密度点处的密度值,re为地球半径,ha,max和ha分别为掩星过程中的最大视高度以及反演大气密度点处的视高度,R(ha1)代表视高度为ha1时折射角大小,kG-D代表Gladstone-Dale常数。
所述步骤(5)中,反演大气密度点处的三维坐标计算方法如下:
式中,(xGP,yGP,zGP)代表反演大气密度点的三维坐标;(x,y,z)为航天器当前位置;(l,m,n)为折射前星光的方向矢量;航天器与反演大气密度点之间的距离Td的确定方法如下:
Td=rs(sinχ-cosχtan R/2) (3)
式中,χ表示航天器轨道平面与折射前星光相切平面的夹角,rs表示航天器地心距。
所述步骤(1)中,掩星开始时航天器所在位置的计算方法如下:
航天器轨道表示为:
ax+by+cz=0 (5)
其中(a,b,c)代表了航天器轨道平面的空间指向;
对于能观测到折射星的航天器位置来说,满足下式:
lx+my+nz+rs sinχ=0 (6)
根据式(4)~式(6)计算得到航天器的位置,(x,y,z)为航天器当前三维坐标;
其中:
α=bl-am
β=cl-an
γ=cm-bn
W=rs sinχ
由式(4)有:
另外,航天器位置还需要满足:
lx+my+nz<0 (10)
本发明与现有技术相比的优点在于:
(1)可以获得三维的、准实时的平流层大气密度模型,对于类似星光折射导航技术这种需要用到高精度全球大气密度模型的领域有不可或缺的作用;
(2)相对现有的大气密度模型获取方法,弥补了直接测量以及地基遥感的空间覆盖率低的缺陷;直接测量及地基遥感只能测量一个或多个点在不同高度剖面的密度值,而本发明可以做到覆盖全球,覆盖率与航天器轨道、星敏感器极限星等有关系,例如选择一颗近极地轨道卫星,极限星等为6,可以在63小时内覆盖全球96%的地区。
(3)无需考虑无线电掩星方法中信号源的成本及星间协同问题。无线电掩星的方法也可以做到全球覆盖,但是该方法需要考虑到星间协同的问题,问题的复杂度会大大提升,成本同样也会提升。而本发明观测的是自然天体无需考虑该问题。
附图说明
图1为本发明的技术要点图;
图2为折射星光与航天器的几何关系;
图3为本发明实施流程框图。
具体实施方式
下面结合附图及实施例对本发明进行详细说明。
如图1,3所示,本发明技术解决方案:基于星光折射的全球大气密度反演方法,步骤如下:
1、对于某一确定的航天器,其轨道方程已知,查找星表可以得到星光矢量方向,由技术要点3中表述的几何关系,计算出某星光掩星开始时航天器的位置p0;
2、从航天器到达步骤1确定的位置p0开始,对星敏感器拍摄星图进行星图处理,星图处理的流程如下:(1)对星敏感器拍摄的星图进行预处理减少噪声,提升信噪比;(2)使用种子生长法进行星图分割,分离星点及背景;(3)利用带阈值的质心法提取星点质心坐标;(4)使用三角形匹配算法进行星图识别,使拍摄星图与星表对应;(5)得到星光折射前后的星光矢量方向,折射角及折射前后星光矢量的夹角;记录折射前星光矢量方向μ0i和计算得到的星光折射角αi,同时记录航天器导航***输出的位置信息xi,持续记录数据直到掩星过程结束;
3、利用步骤2记录的航天器位置xi、折射前星光矢量方向μ0i以及折射角αi,根据几何关系计算视高度hai;
4、利用步骤2记录的折射角αi以及步骤3计算的视高度hai,利用技术要点1给出的基于星光折射的大气密度反演公式式(1),计算每组数据对应的大气密度ρ;根据以上信息代入式(1)可直接求解;
5、利用步骤2记录的航天器位置xi、折射前星光矢量方向μ0i以及折射角αi,根据技术要点2给出的方法计算反演大气密度点处的三维坐标(xGP,yGP,zGP)。
如图1所示,本发明的技术要点如下:
1、基于星光折射的大气密度反演模型
本发明是星光折射定位的逆应用,利用已知的航天器位置以及观测星光折射信息反演计算折射高度处的大气密度。基于星光折射的大气密度反演公式为:
式中,ρ代表反演大气密度点处的密度值,re为地球半径,ha,max和ha分别为掩星过程中的最大视高度以及反演大气密度点处的视高度,R(ha1)代表视高度为ha1时折射角大小,kG-D代表Gladstone-Dale常数。
由式(11)可知,大气密度反演计算需要已知掩星过程中的观测视高度及该视高度对应的折射角。折射角可由星敏感器观测得到,视高度的计算方法如下:
2、反演位置确定
根据技术要点1可以计算反演大气密度点处的密度值,为获得全球三维大气密度模型还需要得到反演大气密度点处的三维坐标。计算方法如下:
式中,(xGP,yGP,zGP)代表反演大气密度点的三维坐标;(x,y,z)为航天器当前位置;(l,m,n)为折射前星光的方向矢量;航天器与反演大气密度点之间的距离Td的计算方法如式(14):
Td=rs(sinχ-cosχtan R/2) (14)
式中,χ表示航天器轨道平面与折射前星光相切平面的夹角,具体表示方法见图2,rs表示航天器地心距,R表示星光折射角的大小。
3、掩星开始时航天器位置的确定
大气反演计算是一个积分过程,所以求解密度还需要知道积分开始和结束时航天器的位置。掩星结束根据拍摄的星图很容易辨别,当目标恒星无法被观测到时掩星过程结束;而掩星开始时由于折射很微弱,单纯靠拍摄星图难以确定掩星过程何时开始,所以需要依赖航天器和星光的几何关系进行解算。掩星开始时航天器所在位置的计算方法如下:
航天器轨道可以表示为:
x2+y2+z2=rs 2 (15)
ax+by+cz=0 (16)
其中(a,b,c)代表了航天器轨道平面的空间指向。
对于能观测到折射星的航天器位置来说,满足下式:
lx+my+nz+rs sinχ=0 (17)
根据式(15)~式(17)可以计算得到航天器的位置,计算结果为:
其中:
α=bl-am
β=cl-an
γ=cm-bn
W=rs sinχ
由式(15)有:
另外,航天器位置还需要满足:
lx+my+nz<0 (21)
以上虽然描述了本发明的具体实施方法,但是本领域的技术人员应当理解,这些仅是举例说明,在不背离本发明原理和实现的前提下,可以对这些实施方案做出多种变更或修改,因此,本发明的保护范围由所附权利要求书限定。
Claims (4)
1.一种基于低轨航天器掩星的平流层大气密度反演方法,其特征在于:包括以下步骤:
(1)计算并记录掩星开始时刻航天器的位置;
(2)从航天器到达步骤(1)确定的位置开始,对星敏感器拍摄星图进行星图处理,记录折射前星光矢量方向和星光折射角,同时记录航天器导航***输出的位置信息,持续记录数据直至掩星过程结束;
(3)利用航天器导航***输出的航天器位置、折射前星光矢量方向以及星光折射角的数据,基于几何关系计算视高度;
(4)利用步骤(2)记录的星光折射角数据以及步骤(3)计算的视高度数据,利用基于星光折射的大气密度反演公式计算每组数据对应的大气密度;
(5)利用步骤(2)记录的航天器位置、折射前星光矢量方向以及星光折射角的数据,基于几何关系计算反演大气密度点的三维坐标。
4.根据权利要求1所述的一种基于低轨航天器掩星的平流层大气密度反演方法,其特征在于:所述步骤(1)中,掩星开始时航天器所在位置的计算方法如下:
航天器轨道表示为:
x2+y2+z2=rs 2 (4)
ax+by+cz=0 (5)
其中(a,b,c)代表了航天器轨道平面的空间指向;
对于能观测到折射星的航天器来说,其位置满足下式:
lx+my+nz+rssinχ=0 (6)
根据式(4)~式(6)计算得到航天器的位置,式中(x,y,z)表示航天器三维坐标;
其中:
α=bl-am
β=cl-an
γ=cm-bn
W=rs sinχ
由式(4)有:
另外,航天器位置还需要满足:
lx+my+nz<0 (10)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111183403.0A CN113916717B (zh) | 2021-10-11 | 2021-10-11 | 一种基于低轨航天器掩星的平流层大气密度反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111183403.0A CN113916717B (zh) | 2021-10-11 | 2021-10-11 | 一种基于低轨航天器掩星的平流层大气密度反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113916717A true CN113916717A (zh) | 2022-01-11 |
CN113916717B CN113916717B (zh) | 2023-08-11 |
Family
ID=79239368
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111183403.0A Active CN113916717B (zh) | 2021-10-11 | 2021-10-11 | 一种基于低轨航天器掩星的平流层大气密度反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113916717B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080174487A1 (en) * | 2006-07-31 | 2008-07-24 | Sokolovskiy Sergey V | Method and system for demodulation of open-loop gps radio occultation signals |
CN101270993A (zh) * | 2007-12-12 | 2008-09-24 | 北京航空航天大学 | 一种远程高精度自主组合导航定位方法 |
CN105352500A (zh) * | 2015-10-21 | 2016-02-24 | 北京航空航天大学 | 带天体干扰的自适应选星方法及*** |
US20160327686A1 (en) * | 2015-05-08 | 2016-11-10 | European Space Agency | Interferometric radio occultation |
CN107782288A (zh) * | 2017-09-22 | 2018-03-09 | 北京空间机电研究所 | 基于背景纹影成像的大气扰动式光学监测飞行器的方法 |
JP2019203889A (ja) * | 2018-05-21 | 2019-11-28 | 修 三浦 | オゾン層の観測方法 |
CN111257497A (zh) * | 2020-02-17 | 2020-06-09 | 上海交通大学 | 一种火焰温度场测量装置及方法 |
CN112329354A (zh) * | 2020-10-22 | 2021-02-05 | 北京航空航天大学 | 一种基于进入圆球的火星高层大气风场及密度反演方法 |
CN112880669A (zh) * | 2020-12-14 | 2021-06-01 | 北京航空航天大学 | 一种航天器星光折射和单轴旋转调制惯性组合导航方法 |
CN113220739A (zh) * | 2021-04-13 | 2021-08-06 | 中国科学院国家空间科学中心 | 一种基于多源gnss掩星数据融合的气候数据反演方法及*** |
-
2021
- 2021-10-11 CN CN202111183403.0A patent/CN113916717B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080174487A1 (en) * | 2006-07-31 | 2008-07-24 | Sokolovskiy Sergey V | Method and system for demodulation of open-loop gps radio occultation signals |
CN101270993A (zh) * | 2007-12-12 | 2008-09-24 | 北京航空航天大学 | 一种远程高精度自主组合导航定位方法 |
US20160327686A1 (en) * | 2015-05-08 | 2016-11-10 | European Space Agency | Interferometric radio occultation |
CN105352500A (zh) * | 2015-10-21 | 2016-02-24 | 北京航空航天大学 | 带天体干扰的自适应选星方法及*** |
CN107782288A (zh) * | 2017-09-22 | 2018-03-09 | 北京空间机电研究所 | 基于背景纹影成像的大气扰动式光学监测飞行器的方法 |
JP2019203889A (ja) * | 2018-05-21 | 2019-11-28 | 修 三浦 | オゾン層の観測方法 |
CN111257497A (zh) * | 2020-02-17 | 2020-06-09 | 上海交通大学 | 一种火焰温度场测量装置及方法 |
CN112329354A (zh) * | 2020-10-22 | 2021-02-05 | 北京航空航天大学 | 一种基于进入圆球的火星高层大气风场及密度反演方法 |
CN112880669A (zh) * | 2020-12-14 | 2021-06-01 | 北京航空航天大学 | 一种航天器星光折射和单轴旋转调制惯性组合导航方法 |
CN113220739A (zh) * | 2021-04-13 | 2021-08-06 | 中国科学院国家空间科学中心 | 一种基于多源gnss掩星数据融合的气候数据反演方法及*** |
Non-Patent Citations (3)
Title |
---|
XIAOLIN NING ET AL.: "Autonomous satellite navigation using starlight refraction angle measurements", 《ADVANCES IN SPACE RESEARCH》, pages 1761 - 1772 * |
王鑫 等: "GPS无线电掩星技术反演大气参数方法对比", 《地球物理学报》, vol. 50, no. 2, pages 346 - 353 * |
胡雄 等: "火星电离层无线电掩星探测仿真研究", 《天文学报》, vol. 50, no. 3, pages 301 - 311 * |
Also Published As
Publication number | Publication date |
---|---|
CN113916717B (zh) | 2023-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hajj et al. | A technical description of atmospheric sounding by GPS occultation | |
Zou et al. | A ray‐tracing operator and its adjoint for the use of GPS/MET refraction angle measurements | |
Hajj et al. | COSMIC GPS ionospheric sensing and space weather | |
Hajj et al. | Ionospheric electron density profiles obtained with the Global Positioning System: Results from the GPS/MET experiment | |
Hinson et al. | Initial results from radio occultation measurements with Mars Global Surveyor | |
Yuan et al. | Sensing climate change using the global positioning system | |
Kursinski et al. | Observing Earth's atmosphere with radio occultation measurements using the Global Positioning System | |
Liao et al. | Preliminary validation of the refractivity from the new radio occultation sounder GNOS/FY-3C | |
Xie et al. | Profiling the atmosphere using the airborne GPS radio occultation technique: A sensitivity study | |
US8519886B2 (en) | System and method for using GPS interferometry to determine atmospheric conditions | |
Muradyan et al. | GPS/INS navigation precision and its effect on airborne radio occultation retrieval accuracy | |
Guerova et al. | Global navigation satellite system monitoring of the atmosphere | |
CN113916717B (zh) | 一种基于低轨航天器掩星的平流层大气密度反演方法 | |
Kolosov et al. | Radio occultation of the Venusian atmosphere and bistatic radiolocation of the surface of Venus using the Venera‐9 and Venera‐10 satellites | |
O'Sullivan et al. | Retrieval of water vapor profiles from GPS/MET radio occultations | |
Norman et al. | Simulating the impact of refractive transverse gradients resulting from a severe troposphere weather event on GPS signal propagation | |
Xie et al. | Sensitivity of airborne radio occultation to tropospheric properties over ocean and land | |
Ghoniem et al. | GNSS-RO LEO satellite orbit optimization for Egypt and the Middle East region | |
McCoy | Space weather comes of age: new sensors and models for ionospheric specification and forecast | |
Aoyama et al. | Down looking GPS occultation measurement on the top of Mt. Fuji | |
Webb | Kinematic GNSS tropospheric estimation and mitigation over a range of altitudes | |
Mousa et al. | Arabic Summary | |
Dodson et al. | Global positioning system as a passive integrated atmospheric water vapor sensing device | |
Mousa et al. | Retrieval of Key Climate Variables Using Occultation Geometry of a Mountain top GPS Receiver. | |
Kursinski et al. | Atmospheric profiles from radio occultation measurements of GPS satellites |
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 |