国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

低軌導(dǎo)航增強(qiáng)GNSS發(fā)展綜述

2019-09-27 06:44:34張小紅馬福建
測繪學(xué)報(bào) 2019年9期
關(guān)鍵詞:定軌星座軌道

張小紅,馬福建

1. 武漢大學(xué)測繪學(xué)院,湖北 武漢 430079; 2. 地球空間環(huán)境與大地測量教育部重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430079

全球?qū)Ш叫l(wèi)星系統(tǒng)(GNSS)是重要的時(shí)空信息基礎(chǔ)設(shè)施,在國民經(jīng)濟(jì)建設(shè)與國防安全領(lǐng)域發(fā)揮著十分重要的作用[1-2],已廣泛應(yīng)用于導(dǎo)航、定位和授時(shí)的眾多領(lǐng)域,其應(yīng)用“僅受限于人們的想象力”。近年來,隨著5G/6G、物聯(lián)網(wǎng)、人工智能和無人駕駛等技術(shù)的發(fā)展,社會生產(chǎn)和生活對精準(zhǔn)時(shí)空信息的需求達(dá)到了前所未有的高度,已從過去的粗略、事后、靜態(tài)和區(qū)域,發(fā)展到現(xiàn)在的精準(zhǔn)、實(shí)時(shí)、動態(tài)和全球[3]。以無人駕駛汽車為例,不僅需要實(shí)時(shí)車道級的導(dǎo)航精度,更需要全路況的連續(xù)可用。但是,衛(wèi)星導(dǎo)航系統(tǒng)自身還無法實(shí)現(xiàn)高精度的室內(nèi)外無縫連續(xù)可用。首先,GNSS基本導(dǎo)航服務(wù)能提供的定位精度通常只有10 m左右,無法滿足高精度用戶的需求;其次,微弱的GNSS信號,不足以穿透物理遮蔽,無法在室內(nèi)、森林和城市峽谷等信號遮蔽區(qū)域提供可靠連續(xù)的導(dǎo)航定位服務(wù);此外,GNSS信號功率很低,容易受到干擾和欺騙,存在一定的安全隱患[4-5]。如何提升衛(wèi)星導(dǎo)航系統(tǒng)的服務(wù)性能,一直是衛(wèi)星導(dǎo)航技術(shù)發(fā)展的驅(qū)動力之一。

1 GNSS增強(qiáng)系統(tǒng)

為了提高衛(wèi)星導(dǎo)航定位的精度、可用性和完好性,在基本導(dǎo)航系統(tǒng)的基礎(chǔ)上,提出并研發(fā)了多種衛(wèi)星導(dǎo)航增強(qiáng)系統(tǒng)。按照增強(qiáng)方式可分為信息增強(qiáng)和信號增強(qiáng)兩大類。

1.1 信息增強(qiáng)與信號增強(qiáng)

信息增強(qiáng)是指通過修正衛(wèi)星導(dǎo)航定位系統(tǒng)的誤差來提高導(dǎo)航定位精度和可靠性等的一種技術(shù)方式。信息增強(qiáng)不提供額外的距離觀測量,只提供消除GNSS誤差的修正信息和完好性信息。信息增強(qiáng)通常需要的是傳輸信道,能夠把增強(qiáng)信息播發(fā)給用戶。按照增強(qiáng)信息傳輸?shù)钠脚_方式,可以分為地基增強(qiáng)和星基增強(qiáng)。網(wǎng)絡(luò)RTK、星基差分等技術(shù)是典型的信息增強(qiáng)系統(tǒng)。

信號增強(qiáng)是指通過除導(dǎo)航衛(wèi)星以外的平臺來發(fā)射導(dǎo)航信號,用戶同時(shí)可以接收導(dǎo)航衛(wèi)星系統(tǒng)本身的導(dǎo)航信號以及其他額外的導(dǎo)航信號,進(jìn)而提高導(dǎo)航定位精度和可用性的方法。信號增強(qiáng)通常需要信號發(fā)射機(jī)來為用戶提供測量信息。信號增強(qiáng)系統(tǒng)提供觀測量,可以與GNSS聯(lián)合定位或者獨(dú)立定位。按照發(fā)射導(dǎo)航信號平臺的位置,信號增強(qiáng)也可以分為地基增強(qiáng)和星基增強(qiáng)。

信息增強(qiáng)可以改善GNSS定位性能,但是對于收不到GNSS信號的情況來講,如室內(nèi)或者遮擋環(huán)境下,信息增強(qiáng)就無能為力了。而信號增強(qiáng)能夠提供新的觀測量,對于GNSS無能為力的場景,如立交橋下、城市峽谷地帶等衛(wèi)星數(shù)目可能不夠的情況下,通過信號增強(qiáng)的方式來彌補(bǔ),甚至對于室內(nèi)這種GNSS信號無法到達(dá)的領(lǐng)域,信號增強(qiáng)也能提供解決方案。

1.2 地基增強(qiáng)與星基增強(qiáng)

從20世紀(jì)90年代初開始,國際上先后提出了以信息增強(qiáng)為主、信號增強(qiáng)為輔的GNSS增強(qiáng)系統(tǒng),具有代表性的是美國的WAAS[6]、歐洲的EGNOS[7]等星基增強(qiáng)系統(tǒng),以LAAS[8]為代表的地基增強(qiáng)系統(tǒng);還提出了單純信息增強(qiáng)的GNSS增強(qiáng)系統(tǒng),主要包括利用移動通信播發(fā)增強(qiáng)信息的地基增強(qiáng)系統(tǒng)(如網(wǎng)絡(luò)RTK)和利用通信衛(wèi)星播發(fā)增強(qiáng)信息的星基增強(qiáng)系統(tǒng)(如StarFire[9]、OmniSTAR[10]等星基差分定位系統(tǒng))。在開闊環(huán)境下,RTK技術(shù)能滿足實(shí)時(shí)高精度應(yīng)用需求,但是需要區(qū)域密集參考站網(wǎng)的支持,PPP技術(shù)雖然不依賴于密集的參考站,但其定位的首次初始化時(shí)間較長,且信號失鎖后的重新初始化時(shí)間與首次初始化時(shí)間幾乎一樣長[11],還難以滿足實(shí)時(shí)性的要求。特別是在高樓和立交橋日益增加的城市環(huán)境下,大部分衛(wèi)星信號被遮擋,RTK和PPP均無法提供連續(xù)可用的導(dǎo)航定位服務(wù)。

1.3 對現(xiàn)有GNSS增強(qiáng)系統(tǒng)的討論

不管是信息增強(qiáng)還是信號增強(qiáng),采用地基增強(qiáng)GNSS的系統(tǒng)的覆蓋范圍有限,無法提供全球無縫統(tǒng)一的高精度服務(wù)。傳統(tǒng)的以信息增強(qiáng)為主、信號增強(qiáng)為輔的GNSS星基增強(qiáng)系統(tǒng),除了日本的QZSS,均將GEO衛(wèi)星作為增強(qiáng)信息播發(fā)平臺,主要是利用其對地靜止和覆蓋面廣的優(yōu)勢,服務(wù)區(qū)內(nèi)用戶可以始終接收該衛(wèi)星的增強(qiáng)服務(wù),但是GEO軌道資源非常有限、轉(zhuǎn)發(fā)通信時(shí)延高,且僅能覆蓋南北緯72°以下區(qū)域[12],而實(shí)際上兩極地區(qū)對于SBAS服務(wù)有強(qiáng)烈的需求——北冰洋地區(qū)航海貿(mào)易隨著全球變暖導(dǎo)致的夏季海冰消融而日漸繁榮[13],極地航線航程短節(jié)約燃油但對完好性要求高[14];QZSS目前由3顆高軌IGSO(偏心率約0.075)和1顆GEO衛(wèi)星組成,通過將IGSO的遠(yuǎn)地點(diǎn)置于日本上空,可以保證該區(qū)域天頂附近始終可見1顆高傾角衛(wèi)星,從而進(jìn)行信息和信號增強(qiáng),QZSS IGSO能夠覆蓋極區(qū),為極地星基增強(qiáng)提供了解決方案之一[15]。

但是,這些已經(jīng)建立的星基增強(qiáng)系統(tǒng)因其能播發(fā)增強(qiáng)信號的衛(wèi)星數(shù)量有限,信號增強(qiáng)的能力較為有限。特別是隨著BDS、Galileo等全球系統(tǒng),以及QZSS、IRNSS等區(qū)域衛(wèi)星導(dǎo)航系統(tǒng)的建成,對多系統(tǒng)用戶來說,傳統(tǒng)以高軌衛(wèi)星為信號播發(fā)載體建成的星基增強(qiáng)系統(tǒng),其信號增強(qiáng)的貢獻(xiàn)比較有限,在信號遮蔽區(qū)域和室內(nèi)仍然無法提供連續(xù)可用的導(dǎo)航定位服務(wù)。為了克服GNSS的脆弱性和局限性,近年來,國內(nèi)外學(xué)者提出通過建立低軌導(dǎo)航星座,利用低軌衛(wèi)星播發(fā)導(dǎo)航信號來增強(qiáng)GNSS的低軌導(dǎo)航增強(qiáng)的構(gòu)想[16-19]。

2 低軌導(dǎo)航星座發(fā)展概況

1963年,美國海軍研制的第1代衛(wèi)星導(dǎo)航系統(tǒng)——子午衛(wèi)星系統(tǒng)(Transit)就是采用低軌星座,由5—10顆運(yùn)行在高度約1100 km圓形極軌道上的低軌衛(wèi)星組成[20]。相比于傳統(tǒng)的地基無線電導(dǎo)航系統(tǒng),Transit在定位精度方面有了明顯改善[21-22],但無法實(shí)現(xiàn)瞬時(shí)單歷元定位,一般需要10~16 min的完整觀測,才能達(dá)到必要的幾何圖形條件。同時(shí),受星座衛(wèi)星數(shù)量少的影響,用戶往往需要等待30~100 min,才能觀測到第2次衛(wèi)星入境[20]。1967年,蘇聯(lián)也開始部署自己的軍方導(dǎo)航通信系統(tǒng)(Parus/Tsikada),在軌道類型、信號頻率、定位方式等方面都類似于Transit[23]。

直到1974年第1臺星載原子鐘的問世[24],以及偽隨機(jī)噪聲碼擴(kuò)頻調(diào)制和多址接入技術(shù)的發(fā)展[25],美國和蘇聯(lián)相繼建設(shè)GPS和GLONASS系統(tǒng)。為了盡快實(shí)現(xiàn)全球任意地區(qū)四重以上覆蓋,綜合考慮系統(tǒng)的建設(shè)和維護(hù)成本,GPS和GLONASS均采用中軌道星座,而沒有采用低軌星座。

隨著通信業(yè)務(wù)的發(fā)展,在20世紀(jì)末,出現(xiàn)了用于移動通信的低軌衛(wèi)星星座,典型代表是美國的銥星(Iridium)和全球星(Globalstar)星座。第1代銥星系統(tǒng)于1998年5月建成,空間段由66顆低軌衛(wèi)星均勻分布在6個近圓極軌道上,每個軌道11顆,高度780 km,傾角86.4°,順行軌道面夾角31.6°,逆行軌道面夾角22°,順行軌道面間相鄰主星相位差16.4°(文獻(xiàn)[26])。銥星運(yùn)行速度7.5 km/s,軌道周期100 min,產(chǎn)生的多普勒頻移量為±40 kHz。每顆衛(wèi)星形成48個點(diǎn)波束。星間鏈路使用Ka頻段與周邊4顆衛(wèi)星進(jìn)行通信;地面上行和下行鏈路也使用Ka頻段;衛(wèi)星語音和數(shù)據(jù)信號發(fā)射使用L波段(1616~1 626.5 MHz),采用符號率25 ksym/s的正交相移鍵控調(diào)制方式,信號幀長度90 ms,目前由Satelles公司取得唯一授權(quán)發(fā)射專用STL(satellite time and location)脈沖信號,既可以提供獨(dú)立的導(dǎo)航、定位和授時(shí)(PNT)服務(wù),也能對GNSS系統(tǒng)進(jìn)行信號增強(qiáng)(http:∥www.satellesinc.com/wp-content/uploads/2016/06/Satelles-White-Paper-Final.pdf)。銥星采用恒溫晶體振蕩器(OCXO)產(chǎn)生通信信號和維持系統(tǒng)時(shí),短期內(nèi)性能很穩(wěn)定,100 s后鐘差和鐘漂穩(wěn)定性低于星載原子鐘,銥星鐘差穩(wěn)定度一般優(yōu)于10-11,1000 s前后分別表現(xiàn)為閃變噪聲和隨機(jī)游走噪聲,因而每天至少需要2次時(shí)頻調(diào)整以獲得時(shí)間同步[27]。目前,銥星系統(tǒng)已于2019年1月11日完成了全部衛(wèi)星的更新?lián)Q代,“銥星二代”空間段由66顆主星和9顆備份星構(gòu)成[28]。

2015年以來,國際知名企業(yè)如美國的OneWeb[29]、SpaceX[30,31]和Boeing[32],韓國的Samsung[33]、中國的航天科技集團(tuán)[34]和航天科工集團(tuán)[35]等,先后宣布發(fā)射和部署各自的商用低軌星座,衛(wèi)星數(shù)量由數(shù)十至上萬顆不等,初衷是為了從太空提供全球范圍內(nèi)、無縫穩(wěn)定的寬帶互聯(lián)網(wǎng)通信服務(wù)。低軌星座建設(shè)進(jìn)入蓬勃發(fā)展時(shí)期,如表1所示。文獻(xiàn)[17]從整個系統(tǒng)的架構(gòu),包括星座幾何圖形條件、空間信號測距誤差、星載原子鐘性能和定軌方法等方面,全面探索其擴(kuò)展成為導(dǎo)航增強(qiáng)星座的可能性,結(jié)果令人振奮。

航天科技的“鴻雁”星座和航天科工的“虹云”工程等均將搭載導(dǎo)航增強(qiáng)有效載荷,既能為我國北斗導(dǎo)航衛(wèi)星系統(tǒng)提供增強(qiáng)改正數(shù)和完好性信息,又能夠自主地播發(fā)導(dǎo)航測距信號,增強(qiáng)PNT服務(wù)性能。

2018年3月,中國科學(xué)院光電研究院依托天儀研究院研制的衛(wèi)星平臺,開展了低軌衛(wèi)星導(dǎo)航信號增強(qiáng)在軌試驗(yàn),旨在驗(yàn)證通信與導(dǎo)航增強(qiáng)在信號層面深度融合新體制的功能和性能,探索基于低軌衛(wèi)星導(dǎo)航信號增強(qiáng)的應(yīng)用模式。試驗(yàn)衛(wèi)星軌道高度537 km,導(dǎo)航增強(qiáng)信號采用時(shí)分多載波體制、S頻段播發(fā),信號發(fā)射功率33 dBmW,落地電平-129~-118.5 dBmW,衛(wèi)星過境時(shí)長約10 min(http:∥www.beidou.gov.cn/yw/xydt/201803/t20180330_14297.html)。

2018年6月2日,武漢大學(xué)研制的“珞珈一號”科學(xué)實(shí)驗(yàn)衛(wèi)星01星(Luojia-1A)成功發(fā)射,是我國開展低軌夜光遙感、導(dǎo)航增強(qiáng)衛(wèi)星研究以及集成化空間信息系統(tǒng)建設(shè)的有益探索。“珞珈一號”01星位于高度645 km的太陽同步軌道,同樣使用OCXO代替原子鐘。初步導(dǎo)航性能顯示,星歷內(nèi)符合精度為0.1 m,鐘差穩(wěn)定度為3×10-10,偽距和載波測量精度分別為2.6 m和0.013 m,高仰角條件下優(yōu)于1.5 m和1.7 mm,因而能夠聯(lián)合現(xiàn)有的GNSS增強(qiáng)實(shí)時(shí)導(dǎo)航服務(wù)[44-45]。

2018年下半年,“鴻雁”星座、“虹云”工程和北京未來導(dǎo)航科技有限公司的“微厘空間”,均發(fā)射了各自的首顆驗(yàn)證衛(wèi)星開展測試論證工作。預(yù)計(jì)2023年前后,全球?qū)⑿纬啥鄠€低軌星座增強(qiáng)Multi-GNSS的新局面。

表1 部分已部署或提出的商用低軌星座匯總

3 低軌與中高軌衛(wèi)星導(dǎo)航關(guān)鍵指標(biāo)對比

根據(jù)軌道的偏心率不同,可將其劃分為近圓軌道(偏心率<0.14)和橢圓軌道(偏心率≥0.14),GNSS衛(wèi)星除了Galileo的FOC FM1和FOC FM2衛(wèi)星錯誤入軌至橢圓軌道之外,其余均工作在近圓軌道上。根據(jù)衛(wèi)星的運(yùn)行高度不同,近圓軌道又可進(jìn)一步細(xì)分為低軌、中軌和高軌。低軌道高度通常為300~2000 km,衛(wèi)星運(yùn)行速度7.7~6.9 km/s,運(yùn)行周期90.4~127.1 min;中軌道高度通常為2000~35 700 km,運(yùn)行周期約2~24 h,受范艾倫輻射帶高能粒子輻射的影響[46],很少有航天器在該區(qū)域活動,但GNSS衛(wèi)星卻普遍位于20 000 km左右,運(yùn)行速度3.9 km/s,運(yùn)行周期約12 h;高軌是指高度在35 700 km以上的軌道,典型的高軌衛(wèi)星是地球同步軌道衛(wèi)星,如BDS GEO和IGSO。截至2018年11月30日,全球在軌衛(wèi)星數(shù)量共計(jì)1957顆[47],包括1232顆LEO、126顆MEO、558顆地球同步軌道和41顆橢圓軌道衛(wèi)星,如圖1所示。

圖1 空間輻射電離總劑量及在軌衛(wèi)星數(shù)量分布Fig.1 Distributions of the total ionizing dose of space radiation and operational satellites on orbit

與中高軌衛(wèi)星相比,低軌衛(wèi)星具有距離地面近、運(yùn)行速度快、受到的攝動影響大這3大特點(diǎn)。

3.1 低軌星座優(yōu)勢

首先,軌道越低,克服重力勢能做功到達(dá)指定高度所需的發(fā)射速度和能量越小,需要的燃料越少,運(yùn)載火箭易于小型化。此外,輕型化的低軌小衛(wèi)星批量制造成本低,易于攜帶,可采用一箭多星技術(shù)減少發(fā)射次數(shù)。2017年1月14日,SpaceX就采用一箭十星的方式將“銥星二代”(Iridium NEXT)首批10顆低軌衛(wèi)星送入太空(http:∥www.gpsworld.com/spacex-launches-first-batch-of-iridium-next-satellites/),大大降低了衛(wèi)星發(fā)射成本。由于信號傳播距離短,自由空間損耗更少,高度780 km的銥星地面接收信號強(qiáng)度比GPS高約30 dB(1000倍)[48],有助于改善信號受遮蔽環(huán)境下的定位效果,提升抗干擾、防欺騙性能;有利于高速衛(wèi)星通信、寬帶互聯(lián)網(wǎng)接入和數(shù)據(jù)傳輸;有利于對地觀測技術(shù)的發(fā)展,可以獲取高分辨率遙感影像、大氣監(jiān)測信息等。

其次,低軌衛(wèi)星運(yùn)行速度快,多普勒頻移現(xiàn)象明顯,有利于提高測速的精度和基于多普勒觀測值的載波相位周跳探測效果,多普勒信息也可用于定位,如Transit;通常低軌衛(wèi)星在測站上方的可見時(shí)長僅10~20 min,如圖2所示,相同時(shí)間內(nèi)低軌衛(wèi)星劃過的天空軌跡更長,幾何圖形變化快,使得定位過程中歷元間觀測方程的相關(guān)性減弱[49],參數(shù)的可估性大大增強(qiáng),有望從根本上解決載波相位模糊度參數(shù)收斂和固定慢的問題,進(jìn)而實(shí)現(xiàn)快速精密定位;短時(shí)間內(nèi)反射信號不再是靜態(tài)的,與直達(dá)信號更易分離,多路徑降相關(guān)時(shí)間短,對信號相關(guān)峰和正交性影響減小,抗多路徑效果好。

圖2 GPS和低軌衛(wèi)星在12 min內(nèi)劃過的天空軌跡對比Fig.2 Comparison of 12-min sky plots between GPS and LEO satellites

最后,低軌衛(wèi)星軌道低,對重力場敏感度高,有利于反演更高階數(shù)的地球重力場模型。

3.2 低軌星座劣勢

由于衛(wèi)星軌道低,單星地面覆蓋范圍小。圖3顯示高度1100 km的“鴻雁”低軌星座單顆衛(wèi)星覆蓋面積只有BDS GEO衛(wèi)星的1/10,因此,在截止高度角10°的情況下,經(jīng)計(jì)算至少需要54顆衛(wèi)星才能保證全球單重以上覆蓋[50],低緯度地區(qū)僅可見1顆衛(wèi)星,若要滿足連續(xù)四重覆蓋的定位要求,則需要近200顆衛(wèi)星,并且“鴻雁”衛(wèi)星的天底距58°遠(yuǎn)大于BDS GEO的9°,為保證覆蓋區(qū)的信號增益,應(yīng)采用多波束天線[51]。接收機(jī)信號捕獲要求高,增加了接收機(jī)的負(fù)擔(dān),捕獲過程中多普勒搜索范圍更大,捕獲速度降低[52]。低軌衛(wèi)星精密定軌與預(yù)報(bào)難度大。圖4給出了不同高度軌道攝動力的加速度量級[53]??梢钥闯?,低軌所受到的大氣阻力、地球非球形引力和廣義相對論作用均明顯高于中高軌,需要選用更精細(xì)的力模型參數(shù)和處理策略進(jìn)行精密定軌才能達(dá)到厘米級的定軌精度[54],預(yù)報(bào)軌道根數(shù)的擬合時(shí)長需要20 min甚至更短,擬合參數(shù)更多[55-56];由于低軌道的空氣阻力大,衛(wèi)星速度會逐漸降低,軌道高度也會降低,需要頻繁啟動進(jìn)行軌道維持,但攜帶的燃料卻是有限的,低軌衛(wèi)星壽命短。

圖3 “鴻雁”星座增強(qiáng)北斗三代Fig.3 Sketch of the Hongyan constellation-augmented BDS-3

圖4 軌道攝動隨高度的變化[53]Fig.4 Orbital perturbations vary with altitude[53]

4 低軌導(dǎo)航增強(qiáng)帶來的機(jī)遇

與GPS、GLONASS、BDS、Galileo等中高軌導(dǎo)航衛(wèi)星相比,低軌導(dǎo)航信號將給聯(lián)合定軌、快速精密定位、空間天氣監(jiān)測、室內(nèi)定位等方面帶來好處。

4.1 聯(lián)合定軌

傳統(tǒng)的GNSS衛(wèi)星精密定軌是利用全球均勻分布的大量地面監(jiān)測站,對導(dǎo)航衛(wèi)星進(jìn)行偽距和載波相位測量,再結(jié)合精確的軌道動力學(xué)模型和誤差改正模型進(jìn)行數(shù)據(jù)處理,確定GNSS衛(wèi)星的精密軌道。然而,我國BDS監(jiān)測站建設(shè)受地緣因素影響較大,難以實(shí)現(xiàn)全球均勻布站[57]。搭載星載GNSS接收機(jī)的低軌衛(wèi)星可以作為“星基監(jiān)測站”,參與高中低軌衛(wèi)星聯(lián)合定軌,能彌補(bǔ)地面站的不足,極大增強(qiáng)GNSS衛(wèi)星跟蹤網(wǎng)的圖形強(qiáng)度,使軌道和力模型參數(shù)估計(jì)得更準(zhǔn)確,實(shí)現(xiàn)區(qū)域監(jiān)測站條件下的導(dǎo)航衛(wèi)星精密定軌。

2004年,文獻(xiàn)[58]最早提出采用全球分布的40個地面站與3顆低軌衛(wèi)星(CHAMP、GRACE-A/B)聯(lián)合定軌, GPS衛(wèi)星軌道精度由傳統(tǒng)定軌方法的7 cm提高為聯(lián)合定軌的3.7 cm。文獻(xiàn)[59]在PANDA軟件基礎(chǔ)上,開展了類似的研究工作,采用全球分布的21個地面站與3顆低軌衛(wèi)星聯(lián)合定軌,相比于僅使用43個地面站的常規(guī)定軌,GPS軌道精度提升40%。文獻(xiàn)[34,57,60]基于仿真研究表明,在區(qū)域監(jiān)測站條件下,聯(lián)合定軌也能顯著改善BDS尤其是BDS GEO衛(wèi)星的定軌精度。文獻(xiàn)[61]基于仿真分析,提出了適用于大型低軌星座與GNSS星座聯(lián)合定軌的最佳方案,即首先優(yōu)選部分低軌衛(wèi)星,綜合使用星載GNSS數(shù)據(jù)、低軌衛(wèi)星和GNSS衛(wèi)星導(dǎo)航數(shù)據(jù)聯(lián)合定軌,確定參與聯(lián)合定軌的衛(wèi)星的精密軌道和鐘差,再使用常規(guī)方法確定剩下的低軌衛(wèi)星軌道和鐘差,結(jié)果如圖5和圖6所示。該方案可以在保證計(jì)算效率的前提下,獲得理想的定軌和估鐘精度,空間信號測距誤差達(dá)到厘米級,能滿足實(shí)時(shí)精密定位服務(wù)的要求??梢灶A(yù)見,未來整個空間信息網(wǎng)路建成以后,如圖7所示,各類星地鏈路、星間鏈路將會帶來更多的觀測數(shù)據(jù)和更強(qiáng)的幾何圖形條件,進(jìn)一步提升聯(lián)合定軌的性能。

圖5 BDS+GPS+LEO星座聯(lián)合定軌在切向、法向和經(jīng)向分量上的軌道誤差[61]Fig.5 Orbit errors in along, cross, and radial components for BDS+GPS+LEO combined precise orbit determination[61]

圖6 BDS+GPS+LEO星座聯(lián)合定軌的估鐘誤差[61]Fig.6 Clock errors for BDS+GPS+LEO combined POD[61]

圖7 低軌導(dǎo)航和星間鏈路增強(qiáng)的高中低軌衛(wèi)星聯(lián)合定軌Fig.7 Combined POD for high, medium, and low Earth orbit satellites augmented with LEO navigation and inter-satellite link data

4.2 快速精密定位

模糊度參數(shù)的快速解算是實(shí)現(xiàn)快速精密定位的關(guān)鍵。1998年,文獻(xiàn)[16]首次提出利用全球星(Globalstar)輔助GPS整周模糊度解算,設(shè)計(jì)了Globalstar/GPS雙系統(tǒng)接收機(jī),并采用蒙特卡羅模擬測試其性能。結(jié)果表明,低軌衛(wèi)星幾何圖形的快速變化,有助于模糊度的快速解算,1 km和5 km的基線分別能夠在120 s和240 s內(nèi),實(shí)現(xiàn)厘米級的快速差分定位。文獻(xiàn)[62]基于仿真數(shù)據(jù),從基線精度、模糊度精度與相關(guān)性、模糊度精度因子和固定成功率等方面,全面評估銥星增強(qiáng)短基線GPS RTK對模糊度解算效果的提升。文獻(xiàn)[49]采用類似方法,研究了銥星增強(qiáng)BDS的效果。文獻(xiàn)[63]仿真分析了不同衛(wèi)星數(shù)量規(guī)模的低軌星座增強(qiáng)中長基線GPS RTK定位性能。結(jié)果表明,低軌衛(wèi)星數(shù)量充足的情況下,長基線RTK模糊度首次固定時(shí)間平均由12 min縮短至2 min,固定成功率由76%提升至96%。

為了分析低軌導(dǎo)航增強(qiáng)星座對PPP的貢獻(xiàn),2015年,文獻(xiàn)[64]最早仿真比較單GPS、GPS+GLONASS、GPS+LEO PPP的定位精度和收斂時(shí)間。結(jié)果表明,采用GPS+LEO和GPS+GLONASS組合PPP收斂時(shí)間分別縮短51.31%和3.93%,說明低軌衛(wèi)星對PPP收斂速度提升貢獻(xiàn)更大。文獻(xiàn)[18—19,65]基于更精細(xì)仿真的觀測值,評估不同低軌星座增強(qiáng)BDS或Multi-GNSS PPP浮點(diǎn)解收斂性能。結(jié)果表明,PPP收斂時(shí)間可由單系統(tǒng)的30 min顯著縮短至數(shù)分鐘甚至1 min以內(nèi),且收斂速度增強(qiáng)效果主要取決于可見低軌衛(wèi)星數(shù),如圖8所示。通過加入包含60、96、192和288顆星的低軌星座增強(qiáng)Multi-GNSS PPP,中緯度NNOR站上分別可見2.4、3.1、6.3和9.5顆低軌衛(wèi)星,收斂時(shí)間由12.5 min縮短至6.2 min、3.9 min、2.8 min和1.9 min[19]。在此基礎(chǔ)上,文獻(xiàn)[66]進(jìn)一步探究了低軌星座增強(qiáng)Multi-GNSS PPP三頻固定解的性能??傮w而言,在不依賴區(qū)域密集地面參考站增強(qiáng)的條件下,基于低軌全球?qū)Ш叫亲鰪?qiáng)GNSS,有望從根本上解決PPP初始化時(shí)間長的難題,真正發(fā)揮PPP技術(shù)的優(yōu)越性[67]。

4.3 空間天氣監(jiān)測

高中低軌導(dǎo)航星座聯(lián)合,為大氣監(jiān)測提供了新的技術(shù)手段。其優(yōu)勢在于:更多的可用衛(wèi)星,可以提取出數(shù)量更多的傾斜路徑延遲;短時(shí)間內(nèi)能夠提供更多有效的觀測數(shù)據(jù),有利于實(shí)現(xiàn)快速的大氣建模;單位時(shí)間內(nèi)低軌衛(wèi)星劃過的軌跡長,高度角和方位角變化大,使得有效監(jiān)測范圍擴(kuò)大,如圖9所示。初步模擬結(jié)果表明,LEO星群新的電離層觀測可為電離層建模提供極為豐富的觀測資料數(shù)據(jù),特別是可以有效填補(bǔ)70%的海洋上空沒有地基GNSS電離層觀測的資料空白;有助于提取和監(jiān)測對流層梯度、電離層梯度、電離層閃爍等快變參數(shù);BDS GEO衛(wèi)星具有獨(dú)特的靜地特性,有利于特定地區(qū)長時(shí)間的大氣監(jiān)測[68];不同高度衛(wèi)星星座組合,有助于研究電離層的垂直分層結(jié)構(gòu);低軌星座星載GNSS數(shù)據(jù),可用于提取等離子體層的電子含量[69]。

4.4 室內(nèi)定位

低軌衛(wèi)星距離地球表面近,地面接收信號強(qiáng)度高,有利于改善信號受遮蔽環(huán)境下的定位性能。2016年,Satelles公司對某高層建筑物STL室內(nèi)服務(wù)性能進(jìn)行了測試,選用獨(dú)立溫控晶體振蕩器(TCXO)作為接收機(jī)鐘,如圖10和圖11所示。結(jié)果表明,對于GPS,只有在最頂層(13樓)靠窗處才能接收到1—2顆衛(wèi)星的信號,其余層則接收不到信號;而對于銥星,其信號強(qiáng)度高,能夠穿透多層鋼筋混凝土材料的阻隔,即使在第2層也能獲得35~55dB-Hz的載噪比,相當(dāng)于GPS在野外開闊環(huán)境下的信號功率水平,STL室內(nèi)定位精度為20 m、授時(shí)精度為亞μs級[70]。2018年,STL最新測試結(jié)果表明,若采用差分?jǐn)?shù)據(jù)并選用精度更高的OCXO作為接收機(jī)鐘,室內(nèi)授時(shí)精度可提升至160 ns(http:∥www.gpsworld.com/satelles-shows-improved-pnt-accuracy-from-leo-con-stellation/)。在野外開闊地帶,STL定位精度不如GNSS高,但由于其采用專用的加密信號,抗干擾、防欺騙能力更強(qiáng),能夠增強(qiáng)GNSS的完好性[71-72]。

綜上,低軌導(dǎo)航增強(qiáng)采用了低軌道,增強(qiáng)信息播發(fā)時(shí)延小、傳輸數(shù)據(jù)量大;增強(qiáng)信號功率強(qiáng),抗干擾、防欺騙性能好,能夠增強(qiáng)室內(nèi)等遮蔽區(qū)域服務(wù)性能;增強(qiáng)信號也能顯著加快精密定位模糊度收斂,為聯(lián)合定軌和空間天氣監(jiān)測等提供更多的有效數(shù)據(jù)源。在覆蓋性方面,雖然低軌衛(wèi)星單星覆蓋面小,但多顆衛(wèi)星構(gòu)成的星座,能夠?qū)⑻峁┌▋蓸O地區(qū)在內(nèi)的全球性信息和信號增強(qiáng)。低軌星座具有地面接收信號強(qiáng)、幾何圖形變化快的優(yōu)勢,能夠與中高軌GNSS星座形成互補(bǔ),有望實(shí)現(xiàn)對導(dǎo)航系統(tǒng)精度、完好性、連續(xù)性和可用性的全面增強(qiáng)。

5 低軌導(dǎo)航增強(qiáng)關(guān)鍵技術(shù)

低軌全球?qū)Ш皆鰪?qiáng)系統(tǒng)由空間段、地面段和用戶段3部分組成??臻g段由數(shù)十至上百顆搭載導(dǎo)航增強(qiáng)有效載荷的低軌衛(wèi)星構(gòu)成,主要任務(wù)是向各類用戶播發(fā)導(dǎo)航信號、高中低軌導(dǎo)航衛(wèi)星增強(qiáng)信息等,具備轉(zhuǎn)發(fā)星和導(dǎo)航星功能;地面段包括地面運(yùn)控系統(tǒng)和地面監(jiān)測站,共同完成在軌衛(wèi)星的運(yùn)行管理和控制;用戶段包括各類型用戶終端、模塊、芯片及配套設(shè)備。低軌導(dǎo)航增強(qiáng)從現(xiàn)在的概念階段到未來實(shí)際運(yùn)行的業(yè)務(wù)系統(tǒng),在空間段、地面段和用戶段均需要突破一系列的關(guān)鍵技術(shù)。

5.1 空間段

星座構(gòu)型的設(shè)計(jì)與優(yōu)化,是空間段建設(shè)必須解決的首要問題,直接決定了增強(qiáng)系統(tǒng)的成本、覆蓋性能、幾何圖形強(qiáng)度和服務(wù)能力,應(yīng)當(dāng)在盡可能節(jié)約成本的前提下,利用具有相似類型和功能的多顆低軌衛(wèi)星,分布在相似或互補(bǔ)的軌道上,共享控制,協(xié)同完成導(dǎo)航增強(qiáng)任務(wù)。

低軌衛(wèi)星數(shù)目眾多,容易產(chǎn)生多址干擾,而且由于多頻GNSS普遍采用L波段,導(dǎo)致該頻段上信號日趨擁擠,信號干擾問題日益嚴(yán)峻。因此,需要設(shè)計(jì)具有一定抗干擾性能的新信號體制,開發(fā)新的導(dǎo)航頻段(如S頻段、C頻段),研究適用低軌導(dǎo)航增強(qiáng)信號的擴(kuò)頻碼優(yōu)化、信號調(diào)制、捕獲跟蹤、信道編碼和多路復(fù)用等技術(shù)。

圖8 不同衛(wèi)星數(shù)量規(guī)模的低軌星座增強(qiáng)Multi-GNSS PPP[19]Fig.8 Different satellite numbers of LEO constellation-augmented Multi-GNSS PPP[19]

圖9 Multi-GNSS與192顆星的低軌星座在1 h內(nèi)產(chǎn)生的地基電離層穿刺點(diǎn)分布對比Fig.9 Comparison of the distributions of 1 h ground-based ionospheric pierce points between the multi-GNSS and 192 satellites LEO constellation

圖10 銥星STL測試場景[70]Fig.10 Test scenario of Iridium based STL[70]

如何最大限度地發(fā)揮低軌星座自身優(yōu)勢,實(shí)現(xiàn)導(dǎo)航與通信的融合也是關(guān)鍵。目前已經(jīng)可以利用移動通信網(wǎng)絡(luò)來播發(fā)GNSS增強(qiáng)信息??梢哉J(rèn)為,衛(wèi)星導(dǎo)航與移動通信已經(jīng)實(shí)現(xiàn)了部分融合。由于低軌導(dǎo)航增強(qiáng)今后主要將依托于互聯(lián)網(wǎng)星座進(jìn)行推廣,未來一定會進(jìn)一步加強(qiáng)導(dǎo)航和通信在信息層面甚至信號層面的深度融合,通過導(dǎo)航增強(qiáng)信號與移動通信信號的聯(lián)合設(shè)計(jì),將導(dǎo)航信號的所有功率全部分配給導(dǎo)頻分量,有利于接收機(jī)的跟蹤、捕獲和測距,而電文數(shù)據(jù)則可以通過低軌互聯(lián)網(wǎng)通信來快速播發(fā)[73]。

在衛(wèi)星有效載荷的配置方面,應(yīng)依托低軌互聯(lián)網(wǎng)星座平臺,優(yōu)選部分衛(wèi)星實(shí)現(xiàn)一星多用,通過搭載不同有效載荷,同時(shí)開展導(dǎo)航增強(qiáng)、遙感、衛(wèi)星重力、衛(wèi)星測高、掩星探測、星基GNSS-R等科學(xué)研究。由于低軌區(qū)域的空間電離輻射較小,因而可以使用商用現(xiàn)成品(COTS)器件和芯片級原子鐘(CSAC)以降低成本。

近地空間環(huán)境碎片監(jiān)測也是必須關(guān)注的問題。2009年2月10日,銥星33號曾與一顆俄羅斯報(bào)廢衛(wèi)星相撞形成無數(shù)太空碎片[74],因此,如何避免近地空間數(shù)以千計(jì)的衛(wèi)星之間發(fā)生碰撞以及衛(wèi)星與空間碎片間的碰撞極其關(guān)鍵,需要精確監(jiān)測衛(wèi)星周圍的環(huán)境,及時(shí)機(jī)動規(guī)避風(fēng)險(xiǎn),同時(shí)依賴國際合作,共同遵守《外層空間條約》[75-76]。

5.2 地面段

地面運(yùn)控部分同樣面臨著一些挑戰(zhàn)。由于低軌道衛(wèi)星特殊的動力學(xué)特征,在星歷和鐘差參數(shù)的擬合時(shí)長、更新頻度、參數(shù)個數(shù)等方面,與傳統(tǒng)的中高軌GNSS衛(wèi)星均有所差異;低軌導(dǎo)航增強(qiáng)衛(wèi)星還需要額外播發(fā)增強(qiáng)改正數(shù)和完好性信息,與基本導(dǎo)航數(shù)據(jù)相比,增強(qiáng)改正數(shù)有效時(shí)間短、時(shí)效性強(qiáng),分為快變改正參數(shù)和慢變改正參數(shù),必須分析增強(qiáng)信息的參數(shù)特性,擬合參數(shù)精度指標(biāo)與播發(fā)時(shí)延的關(guān)系,得出各增強(qiáng)參數(shù)的實(shí)時(shí)性要求。低軌導(dǎo)航電文信息及其編排模式、快變和慢變電文參數(shù)等信息需要根據(jù)預(yù)留的信息位重新定義設(shè)計(jì),主要取決于系統(tǒng)播發(fā)能力和參數(shù)播發(fā)的實(shí)時(shí)性需求兩者間的協(xié)調(diào)關(guān)系。

低軌衛(wèi)星運(yùn)行速度快、過境時(shí)間短、信息傳輸時(shí)間相應(yīng)縮短、地面覆蓋小,在電文上注時(shí)需要更多的地面注入站或依靠星間鏈路通信進(jìn)行電文信息傳遞。若采用星間通信,應(yīng)與中高軌衛(wèi)星共同構(gòu)成混合星座,那么低軌衛(wèi)星與中高軌衛(wèi)星間通信方式的最優(yōu)選擇是有待論證和解決的問題。

地面段一項(xiàng)重要任務(wù)是建立與維持低軌導(dǎo)航增強(qiáng)星座的時(shí)空基準(zhǔn)。GNSS與低軌增強(qiáng)數(shù)據(jù)的融合處理,需要在高精度、統(tǒng)一的時(shí)空框架下完成。對于時(shí)間系統(tǒng),需要給出低軌導(dǎo)航增強(qiáng)系統(tǒng)時(shí)定義,并進(jìn)行系統(tǒng)時(shí)間的建立與維持、系統(tǒng)內(nèi)部時(shí)間同步、系統(tǒng)時(shí)間溯源及時(shí)差預(yù)報(bào)工作;對于坐標(biāo)系統(tǒng),需要給出低軌導(dǎo)航增強(qiáng)系統(tǒng)坐標(biāo)系定義、實(shí)現(xiàn)、更新、維護(hù),及其與其他坐標(biāo)系轉(zhuǎn)換的方法。而時(shí)空基準(zhǔn)的統(tǒng)一,離不開GNSS和低軌星座精密定軌估鐘。

為了獲取所有低軌衛(wèi)星實(shí)時(shí)精密軌道,地面主控站可以采用多種定軌策略,包括區(qū)域監(jiān)測站條件下GNSS與低軌星座聯(lián)合精密定軌及軌道預(yù)報(bào)、僅使用全球跟蹤站網(wǎng)而不依賴于GNSS的低軌星座定軌、星間鏈路定軌等,由于衛(wèi)星數(shù)量多、計(jì)算量大,單站觀測弧段短,且在局部區(qū)域存在衛(wèi)星出入境不穩(wěn)定、不連續(xù)的情況,因此,需要分析論證不同觀測弧長對定軌精度的影響,同時(shí)開發(fā)高效率的分布式、并行處理算法。定軌和軌道預(yù)報(bào)中最難模型化的是太陽光壓模型,該模型主要受衛(wèi)星平臺和衛(wèi)星姿態(tài)等因素的影響,需要針對不同型號低軌衛(wèi)星求取相應(yīng)的光壓反射系數(shù),研究以太陽光壓反射系數(shù)為主的動力學(xué)模型參數(shù)變化特性,建立與衛(wèi)星姿態(tài)及太陽位置相關(guān)的高精度數(shù)學(xué)模型,提升精密定軌和中長期預(yù)報(bào)精度。

在星地時(shí)間維持和同步方面,由于高性能星載原子鐘的功耗、尺寸、重量和成本等因素的影響,導(dǎo)致其不適用于低軌衛(wèi)星低成本需求,一般使用CASC或OCXO作為替代進(jìn)行時(shí)間維持,為了保證其高精度和時(shí)頻穩(wěn)定性,需要實(shí)時(shí)或近實(shí)時(shí)解算低軌衛(wèi)星鐘差,采用鐘差馴服的方式維持低軌衛(wèi)星時(shí)間基準(zhǔn),對低軌衛(wèi)星上載荷的計(jì)算能力提出了挑戰(zhàn);或者部分衛(wèi)星搭載高性能原子鐘,其余衛(wèi)星通過星間、星地通信實(shí)現(xiàn)時(shí)間同步。

地面控制部分還要負(fù)責(zé)低軌導(dǎo)航增強(qiáng)星座管理和衛(wèi)星控制工作。具體包括:制定長期和中期的任務(wù)計(jì)劃,確保信號的全球覆蓋和連續(xù)性;衛(wèi)星總體星座的維持;按計(jì)劃補(bǔ)充足夠的衛(wèi)星;處理意外事件和失效衛(wèi)星,將它們對服務(wù)的影響降到最?。槐O(jiān)視和控制每顆衛(wèi)星各方面的狀態(tài),保證衛(wèi)星和有效荷載的正常運(yùn)行以及在意外事件發(fā)生時(shí)對一些突發(fā)性事件處理;計(jì)劃并執(zhí)行軌道調(diào)度、平臺維持等;支持在軌軟件維護(hù)[77]。

5.3 用戶段

低軌導(dǎo)航增強(qiáng)對地面用戶軟硬件設(shè)備和數(shù)據(jù)處理方式提出了新的要求。

在低軌衛(wèi)星高動態(tài)條件下,接收機(jī)信號捕獲更加困難,多普勒搜索范圍更大,捕獲速度降低。短時(shí)間內(nèi)站星連線方向的距離變化大,使得接收信號強(qiáng)度變化大,導(dǎo)致偽距噪聲變化不均勻,影響碼觀測值精度,接收機(jī)射頻前端也要改變,同時(shí),站星連線方向加速度也會變化更大,引起載波頻率和碼相位的大幅變化,容易產(chǎn)生信號失鎖,更大的多普勒預(yù)測不確定性和更短的相干積分時(shí)間要求,給信號的高精度穩(wěn)定跟蹤帶來挑戰(zhàn)[78]。低軌衛(wèi)星數(shù)量多,且觀測弧段短、衛(wèi)星切換頻繁,需要提升接收機(jī)的通道數(shù)量、存儲容量、微處理器計(jì)算能力等,優(yōu)化天線單元和接收單元。針對低軌導(dǎo)航可能采用S或C頻段增強(qiáng)信號,需要設(shè)計(jì)L+S或L+C一體化天線,相應(yīng)的射頻、模數(shù)轉(zhuǎn)換、基帶信號處理、捕獲、跟蹤、解調(diào)等方式均需要做出調(diào)整。

在數(shù)據(jù)處理方面,應(yīng)設(shè)計(jì)開發(fā)更加嚴(yán)苛的質(zhì)量控制算法,以診斷衛(wèi)星軌道和衛(wèi)星鐘差產(chǎn)品的異常等。信號失鎖和頻繁的衛(wèi)星切換可能會產(chǎn)生更多的周跳,而低軌衛(wèi)星運(yùn)行速度快,歷元間電離層變化大,傳統(tǒng)周跳探測算法如電離層殘差法將不再適用,需要研究更有效的算法。高中低軌衛(wèi)星融合定位數(shù)學(xué)模型需要進(jìn)一步精化,因?yàn)樾碌膶?dǎo)航增強(qiáng)信號將引入更多與系統(tǒng)、軌道類型、碼類型和頻率等有關(guān)的偏差項(xiàng),需要對各類偏差的時(shí)域和空域特性進(jìn)行仔細(xì)分析和建模,在后續(xù)的數(shù)據(jù)處理中予以改正或估計(jì)。此外,多源異構(gòu)星座觀測值的合理定權(quán),也是需要關(guān)注的問題。

6 結(jié)論與展望

總之,低軌星座具有地面接收信號強(qiáng)度高、幾何圖形變化快等優(yōu)勢,能夠在全球范圍內(nèi)對GNSS形成信息增強(qiáng)和信號增強(qiáng),全面提升衛(wèi)星導(dǎo)航系統(tǒng)的精度、完好性、連續(xù)性和可用性,為高中低聯(lián)合精密定軌、快速精密定位、空間天氣監(jiān)測和室內(nèi)定位等實(shí)際應(yīng)用和科學(xué)研究帶來新的發(fā)展機(jī)遇,有望解決當(dāng)前增強(qiáng)系統(tǒng)在全球覆蓋、低落地功率和廣域精密定位初始化時(shí)間長、可靠性偏低等方面的問題,服務(wù)于未來以電網(wǎng)、銀行、證券、軍事等高價(jià)值安全用戶,以及以自動駕駛、無人機(jī)為代表的實(shí)時(shí)精密定位用戶。隨著5G/6G技術(shù)的發(fā)展,空、天、地、海泛在移動通信網(wǎng)絡(luò)的建立,智能手機(jī)等移動終端處理能力的提升,低軌導(dǎo)航增強(qiáng)最終有望走進(jìn)千家萬戶,實(shí)現(xiàn)大眾應(yīng)用。低軌導(dǎo)航增強(qiáng)也將是我國綜合PNT體系的重要組成部分。當(dāng)然,也如前文所述,低軌導(dǎo)航增強(qiáng)在空間段、地面段和用戶段也還存在許多關(guān)鍵技術(shù)有待突破和研究。

為了滿足智能時(shí)代人們對于精準(zhǔn)時(shí)空信息的需求,在星基與地基增強(qiáng)一體化[3]、“通導(dǎo)”或“通導(dǎo)遙”一體化的基礎(chǔ)上[79],突破衛(wèi)星導(dǎo)航與慣性、激光雷達(dá)、視覺等多源傳感器的自適應(yīng)融合導(dǎo)航的關(guān)鍵性技術(shù)[4,5],也是導(dǎo)航領(lǐng)域的發(fā)展方向。

猜你喜歡
定軌星座軌道
基于單純形法的TLE軌道確定
CryoSat提升軌道高度與ICESat-2同步運(yùn)行
朝美重回“相互羞辱軌道”?
星座
健康女性(2017年3期)2017-04-27 21:34:31
12星座之我愛洗澡
星座
都市麗人(2016年12期)2016-12-14 09:54:15
星座
都市麗人(2016年7期)2016-08-03 10:06:05
導(dǎo)航星座自主定軌抗差濾波算法
太陽軌道器
太空探索(2014年4期)2014-07-19 10:08:58
偽隨機(jī)脈沖在北斗衛(wèi)星精密定軌中的應(yīng)用
廉江市| 翼城县| 宁波市| 休宁县| 介休市| 博湖县| 上虞市| 合川市| 哈巴河县| 海阳市| 卢龙县| 鹿邑县| 南溪县| 镇巴县| 太白县| 姚安县| 芦山县| 凤冈县| 曲松县| 泾源县| 汉中市| 庄浪县| 湘乡市| 绥芬河市| 五峰| 商水县| 盐山县| 蒙山县| 肥乡县| 错那县| 乡宁县| 莒南县| 曲水县| 平乡县| 横峰县| 林芝县| 平安县| 涞源县| 澳门| 双柏县| 鞍山市|