Multi-GNSS精密定軌及精密定位服務(wù)

1、第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京Multi-GNSS衛(wèi)星精密定軌及精密定位服務(wù)第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)”2016年年09月月23-25日于南京日于南京iGMAS武漢大學(xué)分析中心武漢大學(xué)衛(wèi)星導(dǎo)航定位技術(shù)研究中心第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京 Multi-GNSS衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展現(xiàn)狀衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展現(xiàn)狀1 1 內(nèi) 容 Multi-GNSS精密

2、定位服務(wù)系統(tǒng)技術(shù)2 2 Multi-GNSS衛(wèi)星精密定軌3 總結(jié) 4 第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京Galileo系系統(tǒng)統(tǒng)星星座座 Multi-GNSS衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展現(xiàn)狀衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展現(xiàn)狀GPSGPS GalilGalileoeoCOMPASSCOMPASS GLONASSGLONASS實(shí)時(shí)單點(diǎn)定位,精度:4-20米第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京 Multi-GNSS衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展現(xiàn)狀1

3、 1 內(nèi) 容 Multi-GNSS精密定位服務(wù)系統(tǒng)技術(shù)精密定位服務(wù)系統(tǒng)技術(shù)2 2 Multi-GNSS衛(wèi)星精密定軌3 總結(jié) 4 第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京Multi-GNSS精密定位服務(wù)系統(tǒng)精密定位服務(wù)系統(tǒng)美國NavCom STAR FIRE國家和地區(qū)星基增強(qiáng)系統(tǒng)荷蘭Fugro star美國Trimble RTX第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京6Ge et al 2012PPP-RTK (S

4、SR)第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京秒級(jí)更新，形成廣域和局域相秒級(jí)更新，形成廣域和局域相結(jié)合的綜合誤差模型結(jié)合的綜合誤差模型鐘差秒級(jí)，電離層鐘差秒級(jí)，電離層FCBFCB分鐘級(jí)更新分鐘級(jí)更新 軌道每軌道每h更新，坐標(biāo)每周更新更新，坐標(biāo)每周更新數(shù)據(jù)管理數(shù)據(jù)管理 區(qū)域電離層、對(duì)流層模糊度參數(shù)廣域廣域 實(shí)時(shí)鐘實(shí)時(shí)鐘差、電離層、差、電離層、FCB GNSS軌道及軌道及地面坐標(biāo)地面坐標(biāo)引入廣域空間基準(zhǔn)引入廣域空間基準(zhǔn)引入廣域時(shí)間基準(zhǔn)引入廣域時(shí)間基準(zhǔn)分分布布式式單單站站數(shù)數(shù)據(jù)據(jù)預(yù)預(yù)處處理理數(shù)據(jù)流數(shù)據(jù)流干凈數(shù)

5、據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)數(shù)據(jù)流數(shù)據(jù)流數(shù)據(jù)流數(shù)據(jù)流數(shù)據(jù)流數(shù)據(jù)流iGMAS武漢大學(xué)分析中心實(shí)時(shí)數(shù)據(jù)處理武漢大學(xué)分析中心實(shí)時(shí)數(shù)據(jù)處理第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京實(shí)現(xiàn)大規(guī)模用戶信息實(shí)現(xiàn)大規(guī)模用戶信息實(shí)時(shí)播發(fā)實(shí)時(shí)播發(fā)大規(guī)模數(shù)大規(guī)模數(shù)據(jù)接入與據(jù)接入與播發(fā)播發(fā)增強(qiáng)信息編碼播發(fā) 用戶類型及用戶類型及通信鏈路通信鏈路 各類誤差模各類誤差模型型用用戶戶分分類類管管理理用戶用戶干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)干凈數(shù)據(jù)用戶用戶用戶用戶用戶用戶每每 秒更新秒更新生成改正信息生成改正信息iGM

6、AS武漢大學(xué)分析中心實(shí)時(shí)定位服務(wù)武漢大學(xué)分析中心實(shí)時(shí)定位服務(wù)實(shí)時(shí)產(chǎn)品實(shí)時(shí)產(chǎn)品數(shù)據(jù)庫數(shù)據(jù)庫9Real-time orbit & clock correction estimatorReal-time UPD/FCB estimatorReal-time trop. & iono. correction estimatorTrop.Ion.UPD/FCBOrbitClock實(shí)時(shí)應(yīng)用實(shí)時(shí)應(yīng)用服務(wù)器服務(wù)器RTDRTKPPPPPP-ARPPP-RTKGlobal & Regional real-time positioning service實(shí)時(shí)應(yīng)用實(shí)時(shí)應(yīng)用服務(wù)器服務(wù)器實(shí)時(shí)應(yīng)用實(shí)時(shí)應(yīng)用服務(wù)器服務(wù)器第二屆

7、第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京Initial Results of BDS PPP Integer Ambiguity Resolution Positioning Results AnalysisSystemEastcmNorthcmUpcmFloatFixedFloatFixedFloatFixedBDSGPSBDS/GPSBDS/GPS/GLONASS

8、Fig.17 RMS of Float and Fixed Solutions of BDS, GPS, BDS/GPS and BDS/GPS/GLONASS PPPTable.1 Average Values of RMS of Float and Fixed Solutions of BDS, GPS, BDS/GPS and BDS/GPS/GLONASS PPP第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京Initial Results of BDS PPP Integer Ambiguity Re

9、solution Positioning Results AnalysisFig.18 Time Series of Float and Fixed Solutions of BDS-only, GPS-only, BDS/GPS and BDS/GPS/GLONASS PPP at WUHN第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京實(shí)時(shí)實(shí)時(shí)PPP/RTD定位服務(wù)器定位服務(wù)器IP：65端口：2101第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)”

10、 2016年年09月月23-25日于南京日于南京 Multi-GNSS衛(wèi)星導(dǎo)航系統(tǒng)發(fā)展現(xiàn)狀1 1 內(nèi) 容 Multi-GNSS精密定位服務(wù)系統(tǒng)技術(shù)2 2 Multi-GNSS衛(wèi)星精密定軌衛(wèi)星精密定軌3 總結(jié) 4 第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京iGMAS分析中心分析中心Multi-GNSS精密定軌精密定軌iGMAS提供四大衛(wèi)星導(dǎo)航系統(tǒng)BDS/GPS/ GLONASS/Galileo的超快速/快速/最終的精密軌道和鐘差。實(shí)時(shí)應(yīng)用的是6小時(shí)更新的超快速軌道第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)

11、術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京MGEX的的Multi-GNSS精密定軌精密定軌IGS Multi-GNSS Experiment (MGEX): To collate and analyze all available GNSS signals; to provide Multi-GNSS products (orbits, clocks, biases, etc.)Analysis Centers: Center for Orbit Determination in Europe; Deutsches GeoForschungsZen

12、trum Potsdam; CNES/CLS; Japan Aerospace Exploration Agency; Technische Universitat Munchen; Wuhan University第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京iGMAS武漢大學(xué)分析中心提供的軌道產(chǎn)品武漢大學(xué)分析中心提供的軌道產(chǎn)品ftp:/ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京超快速軌道計(jì)算超快速軌道計(jì)算p更新頻率更新頻率 ： 每

13、6小時(shí)p軌道弧長(zhǎng)軌道弧長(zhǎng) ： 48小時(shí)，其中24小時(shí)計(jì)算部分、24小時(shí)預(yù)報(bào)部分p更新時(shí)延更新時(shí)延 ： 滯后2小時(shí)2424小時(shí)小時(shí)2424小時(shí)小時(shí)2 2小時(shí)小時(shí)2424小時(shí)小時(shí)2424小時(shí)小時(shí)2 2小時(shí)小時(shí)00061218000612180606121800061218002424小時(shí)小時(shí)2424小時(shí)小時(shí)2 2小時(shí)小時(shí)1212180006121800062424小時(shí)小時(shí)2424小時(shí)小時(shí)1818000612180006122 2小時(shí)小時(shí)6小時(shí)更新軌道小時(shí)更新軌道第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-8

14、2-8小時(shí)軌道精度小時(shí)軌道精度-GPS-GPS衛(wèi)星衛(wèi)星p統(tǒng)計(jì)GPS軌道預(yù)報(bào)軌道2-8小時(shí)精度第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-82-8小時(shí)軌道精度小時(shí)軌道精度-GLONASS-GLONASS衛(wèi)星衛(wèi)星p統(tǒng)計(jì)GLONASS軌道預(yù)報(bào)軌道2-8小時(shí)精度第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-82-8小時(shí)軌道精度小時(shí)軌道精度-Galileo-Galileo衛(wèi)星衛(wèi)星p統(tǒng)計(jì)Galileo軌道預(yù)報(bào)軌道2-8小時(shí)精度

15、第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-82-8小時(shí)軌道精度小時(shí)軌道精度BDSBDS的的IGSO/MEOIGSO/MEOp統(tǒng)計(jì)BDS軌道預(yù)報(bào)軌道2-8小時(shí)精度IGSO/MEO第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-82-8小時(shí)軌道精度小時(shí)軌道精度 BDSBDS的的GEOGEOp統(tǒng)計(jì)BDS軌道預(yù)報(bào)軌道2-8小時(shí)精度 GEO第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)”

16、2016年年09月月23-25日于南京日于南京2-82-8小時(shí)預(yù)報(bào)精度小時(shí)預(yù)報(bào)精度RMS(mm)RMS(mm)PRNACRG01G0176 76 47 47 10 10 G02G0279 79 47 47 15 15 G03G03263 263 61 61 45 45 G05G0587 87 45 45 13 13 G06G0675 75 46 46 11 11 G07G0786 86 55 55 16 16 G08G0887 87 43 43 13 13 G09G0982 82 46 46 14 14 G10G10367 367 66 66 65 65 G11G1191 91 45 45 1

17、5 15 G12G12117 117 49 49 20 20 G13G1373 73 49 49 14 14 G14G1477 77 60 60 14 14 G15G1584 84 46 46 14 14 G16G1691 91 52 52 15 15 G17G1765 65 48 48 14 14 G18G1887 87 41 41 15 15 G19G1977 77 46 46 14 14 G20G20115 115 50 50 15 15 G21G2188 88 48 48 14 14 PRNACRG22G2294 94 46 46 18 18 G23G2389 89 49 49 13

18、13 G24G2487 87 46 46 16 16 G25G2584 84 48 48 16 16 G26G2678 78 48 48 11 11 G27G2792 92 48 48 15 15 G28G2896 96 41 41 15 15 G29G2983 83 41 41 16 16 G30G3072 72 53 53 14 14 G31G3168 68 44 44 10 10 G32G3280 80 58 58 12 12 PRNACRR01R01133 133 70 70 29 29 R02R02100 100 57 57 16 16 R03R03113 113 62 62 17

19、17 R04R0482 82 70 70 19 19 R05R05107 107 57 57 22 22 R06R0693 93 54 54 26 26 R07R07137 137 68 68 32 32 R08R08114 114 66 66 25 25 R09R09149 149 70 70 33 33 R10R10125 125 79 79 26 26 R11R11114 114 69 69 28 28 R12R12115 115 70 70 20 20 R13R13136 136 63 63 21 21 R14R1494 94 55 55 23 23 R15R1595 95 70 70

20、 26 26 R16R16124 124 63 63 29 29 R17R17125 125 69 69 24 24 R18R1882 82 62 62 33 33 R19R1995 95 50 50 18 18 R20R2096 96 112 112 18 18 R21R2126 26 21 21 17 17 R22R2278 78 38 38 14 14 第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-82-8小時(shí)預(yù)報(bào)精度小時(shí)預(yù)報(bào)精度RMS(mm)RMS(mm)PRNACRE08E08207 207 6

21、7 67 45 45 E09E09213 213 60 60 46 46 E11E11234 234 147 147 69 69 E12E12229 229 93 93 38 38 E14E14203 203 73 73 50 50 E18E18196 196 69 69 37 37 E19E19163 163 94 94 35 35 E22E22141 141 72 72 36 36 E24E24145 145 86 86 33 33 E26E26328 328 90 90 54 54 E30E30134 134 55 55 28 28 PRNACRC01C01846 846 86 86 1

22、68 168 C02C02953 953 70 70 333 333 C03C03762 762 77 77 190 190 C04C041289 1289 88 88 268 268 C05C05994 994 69 69 257 257 C06C06357 357 223 223 79 79 C07C07268 268 161 161 108 108 C08C08381 381 127 127 266 266 C09C09174 174 175 175 68 68 C10C10273 273 123 123 126 126 C15C15282 282 167 167 171 171 C11

23、C11127 127 73 73 29 29 C12C12250 250 86 86 37 37 C14C14161 161 60 60 26 26 第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京超快速軌道計(jì)算超快速軌道計(jì)算p更新頻率更新頻率 ： 每3小時(shí)p軌道弧長(zhǎng)軌道弧長(zhǎng) ： 48小時(shí)，其中24小時(shí)計(jì)算部分、24小時(shí)預(yù)報(bào)部分p更新時(shí)延更新時(shí)延 ： 滯后2小時(shí)2424小時(shí)小時(shí)2424小時(shí)小時(shí)2 2小時(shí)小時(shí)2 2小時(shí)小時(shí)00061218000612180606121800061218002 2小時(shí)小時(shí)12121

24、80006121800062424小時(shí)小時(shí)1818000612180006122 2小時(shí)小時(shí)0303091521030915212 2小時(shí)小時(shí)0909152103091521032 2小時(shí)小時(shí)15152103091521092 2小時(shí)小時(shí)2121030915210309152 2小時(shí)小時(shí)3小時(shí)更新軌道小時(shí)更新軌道第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-52-5小時(shí)軌道精度小時(shí)軌道精度p統(tǒng)計(jì)BDS軌道預(yù)報(bào)軌道2-5小時(shí)精度第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大

25、會(huì)” 2016年年09月月23-25日于南京日于南京2-52-5小時(shí)軌道精度小時(shí)軌道精度p統(tǒng)計(jì)BDS軌道預(yù)報(bào)軌道2-5小時(shí)精度第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京2-52-5小時(shí)軌道跳躍小時(shí)軌道跳躍RMS(mm)RMS(mm)PRNACRC01C011482 1482 146 146 352 352 C02C02911 911 75 75 153 153 C03C031218 1218 71 71 226 226 C04C041727 1727 169 169 300 300 C05C051008

26、1008 98 98 248 248 C06C06391 391 186 186 104 104 C07C07453 453 292 292 85 85 C08C08487 487 151 151 157 157 C09C09315 315 147 147 84 84 C10C10366 366 257 257 123 123 C15C15302 302 208 208 88 88 C11C11181 181 101 101 43 43 C12C12282 282 81 81 35 35 C14C14184 184 72 72 33 33 PRNACRC01C011454 1454 123 1

27、23 227 227 C02C021353 1353 92 92 294 294 C03C031051 1051 96 96 263 263 C04C042367 2367 145 145 405 405 C05C051810 1810 81 81 238 238 C06C06542 542 338 338 142 142 C07C07346 346 210 210 101 101 C08C08486 486 125 125 295 295 C09C09280 280 199 199 130 130 C10C10373 373 90 90 177 177 C15C15297 297 311 3

28、11 102 102 C11C11180 180 120 120 50 50 C12C12346 346 150 150 36 36 C14C14207 207 98 98 38 38 3小時(shí)軌道更新跳躍RMS(mm)6小時(shí)軌道更新跳躍RMS(mm)第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京實(shí)時(shí)可用軌道24_060-2小時(shí)：38mm2-5小時(shí)：59mm24小時(shí)觀測(cè)數(shù)據(jù)軌道確定部分與預(yù)報(bào)部分精度，紅線后為實(shí)時(shí)軌道可用部分。第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)”

29、2016年年09月月23-25日于南京日于南京不同軌道預(yù)報(bào)策略 軌道預(yù)報(bào)24小時(shí)弧段預(yù)報(bào)47小時(shí)弧段預(yù)報(bào)毫毫米米0-22-5014360024057034362044155053852063859ave4057毫毫米米0-22-5013337023234033336042931053541064255ave3439第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京iGMAS武漢大學(xué)分析中心實(shí)時(shí)定軌軟件武漢大學(xué)分析中心實(shí)時(shí)定軌軟件武漢大學(xué)實(shí)時(shí)精密定軌的研究樓益棟（2010）滑動(dòng)窗口法方程疊加法；戴曉蕾（2016）實(shí)

30、時(shí)濾波法；徐小龍，馬宏陽（2016）分塊遞推最小二乘配置法1）兼容以上兩種方法；2）及時(shí)釋放內(nèi)存3）動(dòng)力學(xué)和幾何待估參數(shù)的鏈接和松弛4）模糊度的固定5）數(shù)據(jù)質(zhì)量控制6）超大規(guī)模的數(shù)據(jù)處理（低軌和高軌星群聯(lián)合實(shí)時(shí)處理）第二屆第二屆“ “中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)中國大地測(cè)量與地球物理學(xué)學(xué)術(shù)大會(huì)” 2016年年09月月23-25日于南京日于南京PRNACRC01C01223 223 22 22 122122C02C021506 1506 2727126126C03C033303302525133133C04C0451051062 62 144144C05C0563163135 35 166 166 C06C0610010032327575C07C0715615659599696C08C08808079798383C09C091421421071079696C10C10111111484874 74 C15C151701701501508181C11C1145 45 34347 7C12C1277 77 31312222C14C14636320202929PRNACRC01C01228 228 23 23 230 230 C02C02971 971 25 25 217 217 C03C0