大型稀疏陣列天基雷達系統分析

潘 潔 李道京 周建衛 盧曉春

?

大型稀疏陣列天基雷達系統分析

潘 潔*①②李道京①周建衛①②盧曉春③

①(中國科學院電子學研究所 北京 100190)②(中國科學院大學 北京 100049)③(中國科學院國家授時中心 西安 710600)

天基雷達可以在全球范圍內實現運動目標探測和對地高分辨率成像。該文研究了天基雷達系統問題，基于大型稀疏陣列天線，論述了天基雷達的技術體制和關鍵技術，設計了基于稀疏陣列的MEO軌道X波段天基雷達系統參數并分析了其系統性能，為天基雷達的技術實現提供了一種新的途徑。

天基雷達；稀疏陣列；目標探測；SAR成像；雷達系統

1 引言

天基雷達(Space-Based Radar, SBR) 以衛星作為觀測平臺，可實現對地高分辨率成像和對空間、空中、地面、海面運動目標探測，能連續對陸海空天目標進行全方位的偵察監視，實現天基預警。天基雷達作為獲取信息的核心裝備，在國家戰略防御預警系統中具有重要地位，現已成為當今世界各航天軍事大國的研究重點[1]。

衛星軌道通常按照高度分為低軌(LEO)、中軌(MEO)、同步衛星軌道(GEO)。LEO軌道高度約為700~1200 km，全球覆蓋約需20~30顆衛星；MEO軌道高度約為8000~13000 km，全球覆蓋約需15~25顆衛星；GEO軌道高度為35800 km，全球覆蓋約需3~4顆衛星。較高的衛星軌道，可以提高覆蓋范圍、減少觀測衛星數目、增加駐留時間用于運動目標探測，對天基雷達預警具有重要的意義。

由于探測距離遠，天基雷達需要大功率孔徑積、輕量化、可展開而且收攏體積小的有源相控陣天線[2]。美國空軍實驗室支持研制的X波段有源透鏡天線，發射前可收攏，發射后采用充氣硬化展開，但其全陣的陣面控制精度需要有特殊措施保證[3]。

考慮到星上供電能力有限，天基雷達必須采用大尺寸天線。當軌道高度大于10000 km時，用于運動目標探測的天基雷達，其功率孔徑積應大于65 dBwm2 [1,4]。當天線發射平均功率需控制在20~30 kw量級，其天線面積應在250 m2左右。在采用大面積天線保證功率孔徑積的條件下，進一步增大天線尺寸提高其空間分辨率有助于對運動目標探測定位精度的提升。這種大型天線采用傳統技術很難工程實現。

稀疏陣列天線及其孔徑綜合技術在射電天文望遠鏡、地基/星載/機載成像輻射計和地基對空觀測雷達中已獲得了廣泛的應用。近年來，中國科學院電子學研究所積極開展了稀疏微波成像理論和方法的研究[8]，深入研究了稀疏陣列天線雷達技術，已取得了一定的研究成果[9]。目前稀疏陣列天線應用于機載雷達[10]和艇載雷達[11,12]的研究工作已經展開，本文將相關概念引入天基雷達，論述了基于大型稀疏陣列的天基雷達的技術體制和關鍵技術，設計了基于稀疏陣列的MEO軌道X波段天基雷達系統參數并分析了其系統性能。

2 技術體制和工作原理

(1)概念： 基于大型稀疏陣列的天基雷達是指利用稀疏陣列天線孔徑綜合技術，通過空-時-頻信號處理，以較高的空間分辨率和較少的設備體積重量實現運動目標探測和對地成像的雷達系統。主要由大型稀疏陣列天線和中央電子設備兩大部分組成，其中稀疏陣列天線由多個可2維電掃描的有源相控子陣和天線展開機構組成。

(2)體制： 稀疏陣列天基雷達是各子陣采用正交信號實現寬波束發射，工作在寬發窄收方式下的多發多收體制雷達系統。其全陣高空間分辨率波束形成中的孔徑綜合和窄波束接收，可理解是基于子陣級接收DBF進行處理。

(3)稀疏陣列天線孔徑綜合原理： 利用在空間優化稀疏布局的多發多收結構天線系統，可獲得更小間隔的空間采樣點并產生新的天線相位中心，通過在空間重排處理，可使相位中心的數量和分布情況與滿陣天線相同，在獲取全陣高空間分辨率天線波束的同時，避免產生柵瓣和較高的旁瓣。

(4)大型陣列形變的測量和天線波束形成方法：

從原理上講，只要精確測量和估計陣列形變誤差并實施補償，就有可能形成理想的天線方向圖。由于陣列天線方向圖形成和合成孔徑雷達成像在本質上是一致的，合成孔徑雷達的成像可理解為利用不同時刻的空間采樣實現了大型陣列的瞬時空間采樣。在此基礎上，傳統的天線波束形成問題也可轉化為雷達成像處理問題，故合成孔徑雷達成像處理中的運動補償方法也可以用于大型陣列形變誤差的補償。

(5)特點： 采用稀疏陣列天線，在同樣的發射功率孔徑積條件下，不僅可以大幅度提高有源陣天線的空間分辨率，而且便于發射前收攏和星上展開。

3 關鍵技術和解決途徑

3.1 大型稀疏陣列天線設計

主要內容包括大型天線陣列優化布局技術、高效輕量有源相控子陣技術、接收和發射DBF技術、稀疏陣列天線收攏和展開機構技術等。

隨著技術的發展，目前高效輕量有源相控陣天線(AESA)技術已有長足的發展。AESA的關鍵部件T/R組件，初期是磚塊Brick結構，體積大而重，很難滿足的安裝尺寸要求，生產需要許多手工參與，成本高，不易小型化設計。隨著MMIC, LTCC, HTCC, SOC等的技術進步，美國、德國等西方先進國家已經設計制造出帶輻射陣元立體集成的瓦片式Tile的T/R，并已成功裝備到F-15戰斗機上使用。德國EADS公司已提出星載TerraSAR-X的磚塊Brick式T/R到瓦片Tile式T/R的過渡計劃。

2006年，國外在X波段由多層LTCC+HTCC的帶輻射元的立體集成T/R組件尺寸為15 mm×15 mm×8 mm，重量小于4 g，輸出功率4 W，噪聲系數NF為3.5 dB。

由于瓦片式T/R組件，薄而輕，易小型化設計，便于星上使用。瓦片式T/R組件可以大大減少生產中的手工參與，便于自動化生產和各工序中的檢測，有利于降低成本，改善和提高可靠性。

基于小型芯片T/R組件和平面集成技術，目前國外研制的X波段輕質低功率密度有源相控陣天線，已用于平流層飛艇的驗證試驗[13]，其重量密度達到1.8 kg/m2，在接收狀態下功耗為4.7 W/m2。設計并選擇低峰值功率芯片T/R組件有利于天基雷達天線的小型化和散熱，是目前的發展趨勢[14]。

本文大型稀疏陣列天線，子陣擬設計成剛性結構，全陣擬工作在非剛性狀態下，其收攏和展開方式可和衛星太陽翼類似。

3.2 多子陣相位中心高精度測量和陣列誤差補償

主要內容包括基于分布式位置和姿態測量系統(POS)多子陣相位中心高精度測量技術、基于回波信號的陣列誤差估計技術、大型陣列誤差高精度補償技術等。

本文大型天線全陣擬工作在非剛性狀態下，多子陣相位中心高精度測量和陣列誤差補償就成為一項關鍵技術，擬采用安裝在各子陣上的分布式POS形成多節點測量系統，完成陣列形變誤差測量并引入SAR成像處理中的運動補償方法實施誤差補償，同時考慮使用雷達成像常用的自聚焦技術實施陣列誤差估計和補償。

陣列形變誤差測量精度的要求約為1/16個波長，在C和X波段對應的相對空間位置測量精度分別約為3.0 mm和1.5 mm。文獻[15]分析了分布式POS的測量精度和用于陣列形變誤差測量的可行性，文獻[12]分析了艇載陣列形變誤差對成像的影響并實施了陣列形變誤差補償。

3.3 運動目標探測和對地高分辨率成像

主要內容包括低信噪比運動目標探測技術、高速運動目標探測技術、地面雜波抑制技術、多模式對地高分辨率成像技術等。

遠距離觀測和對低RCS運動目標探測，要求天基雷達具有對低信噪比微弱信號探測能力。雷達對遠距離目標信號可采用長時間相干積累，對近距離目標信號可采用短時間相干積累，并自動實現一定的靈敏度時間控制(STC)功能。通過設置不同的波駐時間，雷達可在作用距離、波束掃描范圍和數據率等參數間進行折中選擇。采用多頻信號和雙頻共軛處理解除回波信號的多普勒模糊，利用雷達成像處理方法完成高速目標大范圍距離徙動校正，解決低信噪比高速運動目標探測問題[16]。為減少波駐時間，可使用稀疏陣列中的密集相鄰子陣做發射DBF，減少瞬時波束覆蓋范圍，提高探測距離。

由于軌道較高，天基雷達的探測對象可包括陸海空天運動目標，應同時考慮對空間目標的成像問題，故在系統設計上對運動目標探測，應考慮設置遠距離窄帶探測和中距離寬帶高分辨率成像模式。

關于對地高分辨率SAR成像模式，由于本文天基雷達天線主要工作在子陣寬波束發射，全陣接收DBF處理方式下，客觀上為多模式對地成像提供了有利條件。

3.4 天基信號高速采集和高時效處理技術

主要內容包括多通道信號高速采集技術、稀疏陣列天線孔徑綜合和接收DBF處理技術、高時效信號處理技術等。

由于天線子陣數量和通道數較多，數據量較大，為保證信號質量和通道一致性，信號采集擬使用中頻采樣，海量信號存儲擬使用高密度固態盤，信號采集設備并入天線子陣設備中。考慮到天基雷達稀疏陣列天線孔徑綜合、接收DBF處理和對運動目標探測的實際應用需求，星上應有足夠的實時和高時效信號處理能力，相關的信號處理硬件結構可采用FPGA+DSP的方式實現。

4 系統參數設計

4.1 波段

波段是雷達的重要參數，使用較短波長可使雷達的空間分辨率和目標探測定位精度較高。使用較長波長(如S波段)可使雷達T/R組件的效率較高，由于波長較長，也可降低對多子陣相位中心測量和陣列誤差補償的精度要求。目前天基雷達具有明顯的向高頻段發展的趨勢，且在高頻段高密度小功率芯片T/R組件更容易工程實現，本文雷達系統參數分析主要圍繞X波段(中心頻率10 GHz)展開。

4.2 稀疏陣列

本文X波段大型稀疏陣列由20個有源相控子陣組成，形成一個長為238 m的稀疏線陣。稀疏線陣在交軌方向展開，在獲得交軌向高空間分辨率的同時，可考慮利用天基雷達和運動目標存在的相對運動，根據SAR/ISAR成像原理提高順軌方向的空間分辨率和對運動目標的探測定位精度。

本文稀疏陣選為主動冗余線列陣，采用模擬退火算法得到的稀疏線陣布陣方式[5,9]為：[1 2 4 5 6 9 15 21 27 33 36 42 48 54 60 63 64 65 67 68]，圖1為其孔徑綜合的結果。

陣列可獲得的等效相位中心信號數量為400，第1次換為滿陣的孔徑綜合使用的相位中心信號數量為135，第2次近似換為滿陣的孔徑綜合使用的相位中心信號數量為127，可以近似得到兩次滿陣孔徑綜合的結果[5]。圖1橫坐標和縱坐標表示陣列位置，圖中第2條橫虛線表示第1次孔徑綜合處理結果，第1條橫虛線表示第2次孔徑綜合處理結果。稀疏陣列天線的參數如表1所示。

圖1 稀疏線陣孔徑綜合結果

表1 稀疏陣列天線的參數

子陣參數數值全陣參數數值 信號波長(m)0.03波束展寬因子1.2~1.4 子陣尺寸(m2)3.5×3.5全陣尺寸(m2)238×3.50 子陣波束寬度(°)0.6×0.6全陣波束寬度(°)0.01×0.60 波束2維掃描范圍(°)±25全陣面積(m2)245 子陣發射平均功率(kW)1全陣發射平均功率(kW)20 子陣數量(個)20功率孔徑積(dBWm2)66.9

對本文子陣T/R規模和功率分析如下：

4.3 信號波形

為實現多發多收，本文天基雷達需采用正交信號波形，擬使用頻分正交信號和碼分正交信號，當距離分辨率為1 m時，其總的信號帶寬應優于200 MHz。

在對地SAR成像模式時，需對稀疏陣做孔徑綜合，使綜合后的相位中心數量和滿陣天線相同，主要使用頻分正交信號。

在運動目標探測模式時，使用碼分正交信號時，通過瞬時孔徑綜合將稀疏陣換成滿陣利用STAP實現運動目標探測；使用頻分正交信號時，由于運動目標信號具有稀疏特性，可將壓縮感知(CS)理論引入稀疏陣列信號處理過程，實施運動目標探測和圖像重構[17]，因此頻分正交信號將是本文天基雷達的首選信號。

圖2為碼分信號和頻分信號與子陣的收發關系示意圖。

圖2 正交信號和子陣的收發關系示意圖

5 系統性能分析

5.1 信噪比和作用距離

本文信噪比分析引入孔徑綜合處理，主要比對了密集陣窄發窄收、密集陣寬發窄收和稀疏陣寬發窄收3種情況。由于運動目標探測模式所需的功率孔徑積遠高于對地SAR成像模式，本文的信噪比分析主要針對多目標搜索模式進行。雷達系統參數如表2所示。

表2 雷達系統參數

參數數值參數數值 噪聲系數Fn(dB)3軌道高度(km)10000 噪聲溫度(K)290目標RCS3 波駐時間tot(s)5地物散射系數-15 系統損耗L(dB)7子陣面積A(m2)3.5×3.5 最大作用距離R(km)12000天線效率50 SAR斜距分辨率1衛星速度(km/s)5

用單脈沖信噪比表示的雷達方程如式(1)：

(1)

其中，為子陣平均發射功率，為波爾茲曼常數。

(1) 密集布陣：

參考條件1：1個子陣自發自收，工作在寬發寬收狀態下，信噪比為：。

參考條件2(A1)：1個子陣發射，20個子陣密集排布，接收做DBF處理，實現寬發窄收，相對于參考條件1，接收增益提高了20倍，信噪比為：。

參考條件3：20個子陣密集排布，發射碼分正交信號，接收做DBF處理，實現寬發窄收，相對于參考條件2，發射功率提高了20倍，對20組碼的信號做相參積累，信噪比可提高20倍，為：。

參考條件4：20個子陣密集排布，實現窄發窄收，適用于單目標跟蹤模式，相對于參考條件1，收/發天線增益和發射功率均增加了20倍，信噪比可提高(功率)×20(發射面積)×20(接收面積)=26.8 dB。

(2)稀疏布陣： 20個子陣稀疏排布，陣列布局如圖1所示，實現多子陣寬發，接收DBF窄收，信噪比推導過程如下：

對于稀疏陣，由于可以得到兩次近似滿陣孔徑綜合的結果，可獲得增益3 dB，則稀疏陣做兩次孔徑綜合的信噪比為：。

總結上述分析結果，當采用碼分信號時，密集陣的信噪比僅比稀疏陣的信噪比高1.61 dB。在上述參數下，各種收發方式的信噪比分析結果如表3所示。

5.2 功率孔徑積

在功率孔徑積分析時考慮的兩種雷達工作模式[18]如下：

(1)多目標搜索模式下，

(3)

當檢測信噪比SNR為13 dB時，由系統參數可知多目標搜索模式下的功率口徑積需PA=63.1 dB W·m2。

(2)SAR成像模式下，

當成像信噪比SNR為15 dB時，由系統參數可知SAR成像模式下功率口徑積需PA2=84.2 dB W·m4。

本文設計的天基雷達的功率孔徑積基本符合MEO軌道使用要求。

5.3 波束覆蓋范圍和空間分辨率

天基雷達主要用于運動目標探測，當軌道高度為10000 km，衛星地速約為5 km/s。當波束擦地角為45°時，在最大作用距離12000 km處，由子陣波束寬度決定的瞬時覆蓋面積約22500 km2(180 km×125 km)。

為減少波駐時間，可考慮利用稀疏陣列中的密集相鄰子陣做發射DBF，減少瞬時波束覆蓋范圍。若要進一步提高作用距離，需考慮再增加波駐時間。通過有源相控子陣波束的2維電掃可擴大探測范圍。

當本文大型稀疏線陣天線在交軌方向展開，雷達工作在MEO軌道10000 km高度時，全陣波束寬度0.01°×0.60°形成星下點空間分辨率為1.75 km(交軌向)×104.70 km(順軌向)。對地面和空中運動目標，利用SAR/ISAR成像處理方法進一步提高順軌方向的空間分辨率到1 km量級是可能的。

6 結束語

研究大功率孔徑積的天基雷達系統具有重要意義，目前中高軌天基雷達[19]已成為新的研究熱點。基于大型稀疏陣列天線，本文論述了天基雷達的技術體制和關鍵技術，設計了基于MEO軌道X波段稀疏陣列天基雷達的系統參數并分析了其系統性能。

采用稀疏陣列天線，在同樣發射功率孔徑積情況下，可將有源陣天線的空間分辨率提高3倍，對提高運動目標遠距離探測定位精度具有重要意義。本文論述的大型稀疏陣列天線為中高軌天基雷達的技術實現提供了一種新的途徑。

表3 各種收發方式的信噪比分析結果

項目密集陣窄發窄收密集陣寬發窄收稀疏陣寬發窄收 信噪比(dB)26.813.812.2 接收分辨率(°)0.036×0.60.036×0.60.01×0.6 技術特點適用于單目標跟蹤模式，密集大陣工程實現困難適用于多目標搜索模式空間分辨率提高3倍，目標探測定位精度高，陣列稀疏布局容易工程實施，信噪比略下降1.6 dB，適用于多目標搜索模式

[1] 賁德, 林幼權. 天基監視雷達[J]. 現代雷達, 2005, 27(4): 1-4. doi: 10.3969/j.issn.1004-7859.2005.04.001.

BEN De and LIN Youquan. Space-based surveillance radar[J]., 2005, 27(4): 1-4. doi: 10.3969/j.issn.1004- 7859.2005.04.001.

[2] 陳升友. 天基雷達大型可展開相控陣天線及其關鍵技術[J]. 現代雷達, 2008, 30(1): 5-8. doi: 10.3969/j.issn.1004-7859. 2008.01.002.

CHEN Shengyou. Large deployable phased array antenna for space-based radar and its key techniques[J]., 2008, 30(1): 5-8. doi: 10.3969/j.issn.1004-7859.2008.01.002.

[3] GUERCI J and JASKA E. ISAT-Innovative space-based- radar antenna technology[C]. 2003 IEEE International Symposium on Phased Array Systems and Technology, Boston, Massachusetts, USA, 2003: 45-51.

[4] 羅倩, 閆鴻慧. 天基雷達系統設計相關問題分析[J]. 艦船電子工程, 2008, 28(5): 118-121. doi: 10.3969/j.issn.1627-9730. 2008.05.033.

LUO Qian and YAN Honghui. Analysis for space-based radar system design issues[J]., 2008, 28(5): 118-121. doi: 10.3969/j.issn.1627-9730.2008.05.033.

[5] MARKUS P, HELMUT S, STEPHAN D,. Imaging technologies and applications of microwave radiometry[C]. European Radar Conference 2004, Amsterdam, The Netherlands 2004: 356-359.

[6] CHEN D F, CHEN B X, and ZHANG S H. Muti-input muti-output radar and sparse array synthetic impulse and aperture radar[C]. International Conference on Radar, Shanghai, China, 2006: 28-31.

[7] MA C, YEO T S, TAN H S,Three-dimensional ISAR imaging using a two-dimensional sparse antenna array[J]., 2008, 5(3): 378-382. doi: 10.1109/LGRS.2008.916071.

[8] 吳一戎, 洪文, 張冰塵, 等. 稀疏微波成像研究進展[J]. 雷達學報, 2014, 3(4): 383-395. doi: 10.3724/SP.J.1300.2014. 14105.

WU Yirong, HONG Wen, ZHANG Bingchen,Current developments of sparse microwave imaging[J]., 2014, 3(4): 383-395. doi: 10.3724/SP.J.1300.2014. 14105.

[9] 李道京, 侯穎妮, 滕秀敏, 等. 稀疏陣列天線雷達技術及其應用[M]. 北京: 科學出版社, 2014: 152-169.

LI Daojing, HOU Yingni, TENG Xiumin,Radar Technology and Its Application in Sparse Array Antenna[M]. Beijing: Science Press, 2014: 152-169.

[10] 侯穎妮, 李道京, 洪文. 基于稀疏陣列和碼分信號的機載預警雷達STAP研究[J]. 航空學報, 2009, 30(4): 732-737. doi: 10.3321/j.issn:1000-6893.2009.04.025.

HOU Yingni, LI Daojing, and HONG Wen. Airborn early warning radar STAP based on thinned array and code division signal[J]., 2009, 30(4): 732-737. doi: 10.3321/j.issn:1000-6893.2009.04. 025.

[11] 滕秀敏, 李道京. 艇載共形稀疏陣列天線雷達成像研究[J]. 電波科學學報, 2012, 27(4): 644-649. doi: 10.13443/j.cjors.2012. 04.015.

TENG Xiumin and LI Daojing. Airship conformal sparse array radar imaging[J]., 2012, 27(4): 644-649. doi: 10.13443/j.cjors.2012.04.015.

[12] 李烈辰. 變換域稀疏壓縮感知雷達成像技術研究[D]. [博士論文], 中國科學院電子學研究所, 2015.

LI Liechen. Study on the technique of sparse compressed sensing radar imaging in transform domain[D]. [Ph.D. dissertation], Institute of Electrics, Chinese Academy of Sciences, 2015.

[13] CLARK T and JASKA E. Million element ISIS array[C]. 2010 IEEE International Symposium on Phased Array Systems and Technology (ARRAY), Boston, Massachusetts, USA, 2010: 29-36.

[14] VAN BEZOUWEN H, FELDLE H, and HOLPP W. Status and trends in AESA-based radar[C]. 2010 IEEE MTT-S International on Microwave Symposium Digest (MTT), California, USA, 2010: 526-529.

[15] 李道京, 滕秀敏, 潘舟浩. 分布式位置和姿態測量系統的概念與應用方向[J]. 雷達學報, 2013, 2(4): 400-405. doi: 10.3724/ SP.J.1300.2013.13086.

LI Daojing, TENG Xiumin, and PAN Zhouhao. The concept and application of distributed POS[J]., 2013, 2(4): 400-405. doi: 10.3724/SP.J.1300.2013.13086.

[16] 李道京, 等. 高分辨率雷達運動目標成像探測技術[M]. 北京: 國防工業出版社, 2014: 89-101.

LI Daojing,. High Resolution Radar Moving Target Imaging Detection Technology[M]. Beijing: National Defence Industry Press, 2014: 89-101.

[17] CETIN M, STOJANOVIC I, ONHON N,. Sparsity- Driven synthetic aperture radar imaging: reconstruction, autofocusing, moving targets, and compressed sensing[J].,2014, 31(4): 27-40. doi: 10.1109/MSP.2014.2312834.

[18] 王文生, 王小謨. 天基雷達的幾個基本問題[J]. 現代雷達, 2007, 29(4): 1-4. doi: 10.3969/j.issn.1004-7859.2007.04.001.

WANG Wensheng and WANG Xiaomo. Several basic questions for Space-based radar[J]., 2007, 29(4): 1-4. doi: 10.3969/j.issn.1004-7859.2007.04.001.

[19] MATAR J, LOPEZ-DEKKER P, and KRIEGER G. Potential and limitations of MEO SAR[C]. EUSAR 2016, Hamburg, Germany, 2016: 1035-1039.

潘 潔： 女，1977年生，高級工程師，博士生，研究方向為稀疏陣列雷達系統.

李道京： 男，1964年生，研究員，博士生導師，研究方向為雷達系統和雷達信號處理.

周建衛： 男，1989年生，碩士，研究方向為外輻射源雷達系統.

盧曉春： 女，1970年生，研究員，博士生導師，研究方向為精密時間信息傳輸與信息處理.

Space-based Radar System Analysis Based on Large Scale Sparse Array

PAN Jie①②LI Daojing①ZHOU Jianwei①②LU Xiaochun③

①(,,100190,)②(,100049,)③(,,,710600,)

Space-Based Radar (SBR) can realize moving target detection and the high resolution imaging on a global scale. Based on large scale sparse array antenna, this paper discusses the space-based radar technology system and key technology. The system parameters are designed and its performance is analyzed for an X band MEO space-based radar. The large scale sparse array antenna provides a new selection for space-based radar.

Space-Based Radar(SBR); Sparse array; Target detection; SAR imaging; Radar system

TN957

A

1009-5896(2016)12-3269-06

10.11999/JEIT160672

2016-06-27；改回日期：2016-11-09；

2016-12-02

潘潔 panj@mail.ie.ac.cn