• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    顧及地磁影響的GNSS電離層層析不等像素間距算法*

    2020-12-07 09:25:48霍星亮劉昊杰
    國防科技大學學報 2020年6期
    關鍵詞:電子密度格網層析

    霍星亮,劉 琦,2,劉昊杰,2

    (1. 中國科學院精密測量科學與技術創(chuàng)新研究院, 大地測量與地球動力學國家重點實驗室, 湖北 武漢 430077;2. 中國科學院大學, 北京 100049)

    全球導航衛(wèi)星系統(tǒng)(Global Navigation Satellite System, GNSS)電離層層析反演技術以其獨特的空間大尺度電子密度三維監(jiān)測能力日益受到重視。但由于衛(wèi)星觀測射線視角有限,水平方向的觀測射線缺乏,同時地面監(jiān)測站布設不均勻且觀測稀疏等因素的多重影響,層析反演模型中的不適定問題是影響GNSS電離層層析反演精度及其推廣應用的重要因素。國內外學者已在GNSS電離層層析反演方法與技術方面開展了系列研究,克服了存在的不適定問題,取得了許多重要成果[1-7]。主要分為兩大類:第一類是以代數重構為代表的行迭代類重構算法,包括加法代數重構(Algebraic Reconstruction Technique,ART)[1]、乘法代數重構(Multiplicative Algebraic Reconstruction Technique,MART)[8-9]、同時迭代重構(Simultaneous Iterative Reconstruction Technique,SIRT)[10-11]等。該類算法要求提供精度相對較高的電離層電子密度初值,通過將GNSS射線在觀測方程組成的超平面內進行投影迭代,逐步縮小觀測值與投影重構值之間的差距,進而估算出最終的電子密度結果。該類方法的主要優(yōu)點是避開了對由于觀測射線不足在層析系統(tǒng)中形成的大型稀疏矩陣的求逆計算,直接利用觀測方程進行迭代反演計算,提高了計算效率和反演結果的穩(wěn)定性,但它對初值模型的依賴性較強,獲得的結果是近似的局部最優(yōu)解。第二類是非迭代算法,包括:正則化算法[4]、正交函數法[12]、奇異值分解[13]與廣義奇異值分解法[14]、混合重建算法[3]等,該類算法通常要求施加一定的約束條件改良或克服GNSS電離層層析系統(tǒng)的不適定性,如認為在電離層層析系統(tǒng)中相鄰“像素”網格內的電子密度具有較好的平滑性或假定電離層狀態(tài)具備線性變化特性等,進而將電離層電子密度反演轉化為附有約束條件的最優(yōu)化問題,從而使得電離層電子密度的反演值逼近實際值。本文將主要圍繞GNSS迭代類重構算法展開討論。

    傳統(tǒng)以像素類為基函數的GNSS層析技術假定待反演區(qū)域的電離層電子密度按照一定的空間格網間隔離散分布,每個像素格網內的電子密度均勻分布。在此基礎上,ART、SIRT、MART等迭代重構模型的共同特點都是以GNSS觀測射線和層析系統(tǒng)像素格網交叉形成的截距長度為權重比,對電離層電子總含量(Total Electron Content,TEC)層析投影重建值和GNSS TEC實測值之間的誤差進行分配。然而,在迭代重構模型層析反演過程中,投影重建的TEC貢獻來源于GNSS射線截距與格網像素內電子密度估算值之間的乘積,GNSS射線截距在迭代反演過程中保持不變,電子密度誤差是TEC實測值與其投影重構估算值差異的主要來源,也就是說,在GNSS層析系統(tǒng)迭代反演誤差中,電子密度誤差是決定因素,GNSS射線截距對誤差起放大作用。因此,傳統(tǒng)ART、SIRT、MART等GNSS層析迭代反演算法中僅以射線截距作為誤差分配準則的唯一要素是不合理的。為此,國內外部分學者考慮在GNSS層析迭代模型中引入電子密度參量[15-17],提高電子密度反演精度。文獻[15-16]在經典ART算法的基礎上通過引入電離層電子密度相對于GNSS觀測視線方向的最大電子密度的比值作為迭代重構時的權重,使得電離層電子密度的修正值成比例于電離層初值模型提供的電子密度;文獻[17]在GNSS每條觀測射線的每輪迭代修正中引入了電子密度參數作為迭代參量加速收斂;文獻[18]則認為該類修正算法的物理意義不太明確。實際上,此類改進方法的本質是確保反演得到的電離層電子密度垂直剖面結構相對于電離層初值背景模型保持穩(wěn)定,彌補水平方向觀測射線的不足,提高GNSS電離層層析系統(tǒng)反演計算效率與精度[19]。近年來,充分考慮電離層實際物理變化,文獻[20]提出以GNSS射線穿越層析像素格網形成的射線截距與對應的像素格網內電子密度乘積為組合自變量建立層析迭代表達式,合理分配TEC實測值與反演值之間的差距,并設計了一組新的松弛因子控制與削弱了噪聲對電子密度反演結果的影響。但需要說明的是,該方法的松弛因子表達式沒有全面反映出不同高度電離層電子密度變化,因此,構造合理有效的GNSS電離層層析反演模型,提高電離層電子密度反演精度,是本文討論的主要內容之一。

    需要指出的是,目前的GNSS電離層層析迭代反演方法都是基于地理坐標系構造,電離層電子密度像素格網按照地理經度、地理緯度與電離層高度劃分層析反演格網系統(tǒng)。地理坐標系下的電離層層析反演算法及其相關約束研究主要考慮了電離層受太陽輻射的影響。但相關研究表明,電離層電子密度結構不僅與太陽活動密切相關,同時也受到地球磁場及其相關的電動力學過程影響。因此,在現有的GNSS TEC建模及電離層高階項改正研究中,許多學者已注意到地球磁場對電離層變化活動的影響。諸如,國際GNSS服務組織下屬的電離層分析中心之一的歐洲定軌中心利用GNSS觀測數據并在日固地磁坐標系下采用球諧函數模型描述全球電離層TEC變化活動[21];著名的經驗電離層模型IRI也提供了地磁坐標系下的電離層TEC及電子密度變化等[22]。因此,基于地理坐標系建立電離層層析系統(tǒng)的像素格網約束,可能會對電子密度施加不符合實際變化的強約束,導致電子密度反演結果發(fā)生畸變。

    此外,許多電離層層析系統(tǒng)反演以等間距的方式劃分電子密度像素格網并假定電子密度在每個像素格網內均勻分布,忽略了電離層電子密度在垂直高度方向上變化較大的特點。這種等間距劃分的像素格網對于高度方向電子密度變化較大區(qū)域(尤其對F2區(qū)域的峰值電子密度)起到了一定平滑作用,容易降低層析技術估算的電離層峰值電子密度精度。文獻[23-24]研究了小范圍內多尺度電離層層析算法對電子密度反演結果的影響。但截至目前,尚未見到在地磁坐標系下以不等間距電子密度像素格網為基礎的GNSS電離層層析反演研究報道,這也是本文討論的主要內容。

    綜上所述,本文提出在地磁坐標系下沿地磁經度、地磁緯度及電離層高度劃分與建立GNSS電離層層析電子密度像素格網系統(tǒng),并考慮在不同的電離層高度方向上劃分不等間距的像素格網(如在電離層電子密度變化較大的F2層峰值高度區(qū)域采用10 km高度的空間格網間距,在電子密度變化較小的其他高度區(qū)域采用30 km或100 km的空間格網間距)。在此基礎上,通過調整GNSS電離層層析迭代算法中的松弛因子,實現對電離層電子密度重構精度和結果平滑度之間的調節(jié)。最后,利用IRI—2007模型采用模擬測試的方法以及利用GNSS數據反演電子密度并與電離層測高儀電子密度實測數據進行對比的方法,驗證了本文提出的顧及地磁影響的GNSS電離層不等像素間距層析方法反演電離層電子密度的有效性與可靠性。

    1 改進的GNSS電離層層析重構反演模型

    基于經典ART算法的GNSS電離層電子密度三維層析反演技術在迭代前利用經驗電離層模型給電離層區(qū)域內的每個電子密度像素賦予一個初值,然后采用迭代方式估計電子密度。每一步迭代估計對應于一條GNSS TEC觀測射線,完成全部的GNSS TEC測量迭代稱之為一輪迭代,并通過多輪迭代的方式逐步改善重構的電離層電子密度。但由于經典ART算法在每輪迭代過程中迭代松弛因子保持不變,且迭代模型僅與對電子密度誤差起放大作用的射線截距權重相關,使得層析重構反演結果精度不高。針對上述問題,文獻[20]提出一種顧及電離層變化的GNSS層析反演算法,引入在層析格網像素內的射線截距與電子密度的乘積為組合自變量(該組合變量表示GNSS TEC在觀測射線方向上對應層析像素格網內的電子總量分量),合理分配不同電子密度像素格網內TEC實測值與反演值之間的誤差。在此基礎上,構造了與電子密度變化相關的松弛因子表達式,抑制傳播噪聲對電子密度估值的影響,實現平衡與調節(jié)層析反演的電子密度精度與結果的平滑程度。上述GNSS電離層層析迭代模型如下:

    (1)

    其中

    (2)

    (3)

    圖1 GNSS射線穿越電離層層析像素格網示意圖Fig.1 Illustration of the GNSS ray passing through an example reconstruction grid of ionospheric tomography

    從上述GNSS電離層層析迭代模型可以看出,當d

    (4)

    綜合式(1)、式(2)和式(4)即為本文改進的GNSS電離層層析重構反演模型(Improved ART, IART)。

    2 地磁坐標系下劃分不等像素間距的電離層層析格網

    基于GNSS三維電離層層析過程中,為了簡化計算,通常需要將待研究的電離層空間進行離散化。GNSS電離層層析中的離散化工作通常是沿著地理經度、地理緯度、高度方向并按照一定的空間間隔距離建立電離層電子密度像素格網(如,將電離層空間沿經度和緯度方向分別按照5°和2.5°的間隔,高度方向以50 km的間隔將待反演區(qū)域離散化為一些小的空間像素格網),同時假定每個小像素格網內的電子密度均勻分布且在較短的反演時間內保持不變??紤]到電離層變化活動易受地球電磁場的影響,區(qū)別于常見的GNSS電離層層析方法在地理坐標系下劃分與建立電子密度像素格網系統(tǒng),本文提出在地磁坐標系下沿著地磁經度、地磁緯度、高度劃分與建立GNSS電離層層析電子密度像素格網,能夠更加合理地約束不同層析系統(tǒng)格網內的電子密度,并提高反演結果的可靠性。

    此外,以等間距方式劃分電子密度像素格網的GNSS電離層層析反演方法假定電子密度在層析像素格網內均勻分布,忽略了電離層電子密度在垂直高度方向上變化較大的特性,這對于電離層電子密度變化較大的峰值電子密度區(qū)域是不合適的。本文考慮電離層電子密度隨高度變化較大的特點,在電離層高度方向上將電子密度變化較大的區(qū)域劃分較小的電子密度像素格網,在電子密度變化較小的區(qū)域劃分較大的像素格網,盡可能在電離層變化較大的高度方向上提高像素格網的空間分辨率(如在電離層電子密度變化較大的峰值電子密度F2區(qū)域采用10 km高度的像素格網間距;在電子密度變化相對較小的高度區(qū)域采用100 km的電子密度像素格網間距)。本文以中國大陸上空的電離層區(qū)域為研究討論,具體劃分不等像素間距的電離層層析格網方法如下:

    1)考慮本文采用的GPS測站分布范圍,選定的經緯度反演區(qū)域范圍為地磁經度150°~200°(對應的地理經度取75° E~135° E)與地磁緯度7°~42°(對應的地理緯度取10° N~55° N)。在對選定的電離層區(qū)域進行空間離散化時,考慮到電離層沿經度方向上的變化活動比緯度方向上的變化活動較弱,在地磁經度和地磁緯度方向上分別取為5°與2.5°的電子密度像素格網間距。

    2)選定的電離層高度區(qū)域為90~1000 km,分別在距離地面90~210 km的高度取30 km的電子密度像素格網間距、距離地面210~400 km的高度取10 km的電子密度像素格網間距、距離地面400~700 km的高度取50 km的電子密度像素格網間距、距離地面700~1000 km高度取100 km的電子密度像素格網間距。

    3 實驗數據及處理策略

    本文采用2011年12月1日、12月2日、12月6日與12月7日的中國陸態(tài)網絡145個基準站的GPS觀測數據(12月3日—12月5日觀測數據較少)進行實驗討論。觀測數據以IGS標準采樣(30 s)獲取,衛(wèi)星截止高度角為15°。分別利用各個觀測站的雙頻載波相位平滑碼數據,并采用“兩步法”消除站星硬件延遲影響[26],計算對應各個GPS觀測站在站星方向上的斜距電離層TEC值。

    本文將分別采用模擬與實測數據對比的方法,驗證本文提出顧及地磁影響的GNSS電離層層析不等像素間距算法的有效性。在模擬比較過程中,利用IRI—2007電離層模型提供的電子密度作為“真值”,并通過增加電離層隨機誤差及電離層電子密度積分模擬不同的GPS站星斜距方向上TEC的“觀測數據”;在實測數據驗證對比過程中,分別利用北京(40.3° N, 116.2° E)和武漢(30.5° N, 114.3° E)電離層測高儀獲得的電子密度作為觀測“真值”。

    為評估不同GNSS電離層層析反演算法,在實驗討論過程中,分別比較了傳統(tǒng)地理坐標系下電子密度像素等間距格網ART算法(方案一)、地理坐標系下電子密度像素等間距格網IART算法(方案二)、地理坐標系下電子密度像素不等間距格網IART算法(方案三)和地磁坐標系下電子密度像素不等間距格網IART算法(方案四)反演的電子密度結果誤差和精度。需要補充說明的是,在方案一和方案二采用的電子密度像素等間距層析反演系統(tǒng)中,電離層反演高度范圍為90~990 km,電子密度像素格網間距為30 km。此外,為了克服由于觀測數據分布不均勻導致的在沒有觀測信息覆蓋的電子密度像素對初值的完全依賴,以及避免像素間電子密度值發(fā)生大的跳變,本文參照文獻[27]與文獻[28]對待反演區(qū)域的電離層電子密度進行了約束處理,提高了觀測信息稀少的電離層像素內電子密度反演精度。

    4 實驗結果和分析

    4.1 模擬比較與研究

    4.1.1 模擬實驗方法

    本文利用IRI—2007電離層模型提供的電子密度模擬實驗中的“真值”,并作為電離層變化趨勢項,在此基礎上引入電離層隨機誤差。根據文獻[29]電離層隨機項誤差的概率分布基本符合正態(tài)分布規(guī)律,并進一步參照文獻[30]和文獻[31]全球GIM TEC精度統(tǒng)計結果,本文研究作簡化處理,在中緯度地區(qū)(本文研究取地磁緯度大于20°)模擬電離層隨機誤差分布滿足期望值為0 TECu、標準差值為2 TECu的正態(tài)分布,在低緯度地區(qū)(地磁緯度小于20°)模擬電離層隨機誤差分布滿足期望值為0 TECu、標準差值為4 TECu的正態(tài)分布。在綜合IRI—2007模型提供的電子密度積分形成電離層TEC趨勢項和模擬的電離層隨機誤差項即可組合出各條GPS觀測射線上的電離層TEC模擬觀測值。在此基礎上,采用四種不同方案電離層層析算法反演電子密度結果,分別給出電離層在高度剖面和電子密度緯度與高度剖面的變化,并從定性分析和定量統(tǒng)計比較的方式,討論不同GNSS電離層層析反演算法的有效性與可靠性。模擬實驗計算過程如下:

    步驟1:構造系數陣:利用中國陸態(tài)網絡基準站和GPS衛(wèi)星坐標,計算反演時段內GPS信號傳播路徑在其所經過的電子密度像素格網內的截距,以此截距構造出電離層層析反演系數矩陣A。

    步驟2:模擬電子密度真值:利用IRI—2007模型模擬待反演時段各個像素內的電離層電子密度,并以此作為電子密度反演真值x真。

    步驟3:模擬電子密度反演所需要的電離層斜距TEC的真值:聯合步驟1與步驟2計算各條射線傳播路徑上的電離層斜距TEC值y真=A·x真。

    步驟4:構造帶誤差的TEC觀測值:顧及實際觀測中觀測噪聲和離散誤差的存在,在模擬計算中,向模擬電離層斜距TEC真值中加入服從正態(tài)分布的隨機誤差e,于是有:y=y真+e。

    步驟5:綜合步驟1~4,獲得模擬重構的電離層層析反演矩陣:y=A·x,其中,A為步驟1構造的觀測矩陣,y為步驟4)模擬的站星TEC值,x為對應的待反演層析系統(tǒng)像素格網內的電子密度。

    步驟6:在步驟5的基礎上,分別采用四種不同方案的GPS電離層層析算法反演電子密度結果,并通過與IRI—2007模型提供的電子密度“真值”進行比較,論證不同GPS電離層層析反演算法的有效性。

    4.1.2 模擬驗證結果

    1)電離層電子密度高度剖面反演結果比較

    考慮論文長度,本文在討論中僅給出2011年12月6日北京時間11:00、13:00和15:00的電離層電子密度剖面反演結果圖為例,其他時間段的結果見電子密度反演誤差與精度統(tǒng)計表。圖2定性地給出分別由傳統(tǒng)代數重構ART算法(方案一,圖中綠色虛線表示)、地理坐標系下電子密度像素等間距格網層析改進的迭代反演(Improved ART, IART)算法(方案二,圖中藍色實線與加號表示)、地理坐標系下電子密度像素不等間距格網層析IART反演算法(方案三,圖中青色實線與方塊表示)和地磁坐標系下電子密度像素不等間距格網層析IART反演算法(方案四,圖中紅色實線與三角形表示)反演的電子密度隨高度變化結果,及其相對于IRI—2007電離層模型提供的電子密度“真值”(圖中黑色實線表示)之間的比較情況。為了方便將本節(jié)模擬結果和本文后續(xù)GNSS實測數據結果進行對比討論,本節(jié)模擬給出北京(40.3°N,116.2°E)與武漢(30.5°N,114.3° E)上空電離層電子密度剖面結果。

    圖2展示了GNSS電離層層析算法四種不同反演方案估算的電離層電子密度隨高度變化剖面結果,可以看出,不同方案的電離層電子密度隨高度變化和電離層測高儀觀測得到的電子密度變化趨勢相同。但可以明顯看出,與IRI—2007模型提供的電離層電子密度“真值”相比較,方案一采用的傳統(tǒng)層析代數重構算法ART在不同時段和不同測站的電子密度反演結果最差,這與該算法對電離層電子密度反演重構值和實測值之間的誤差迭代分配時僅以觀測射線截距作為準則有關,使得電子密度重構結果誤差較大。

    (a) 2011年12月6日北京電離層電子密度模擬結果(a) Simulated ionospheric electron density profiles over Beijing on December 6, 2011

    (b) 2011年12月6日武漢電離層電子密度模擬結果(b) Simulated ionospheric electron density profiles over Wuhan on December 6, 2011圖2 不同GNSS電離層層析算法重構的電子密度與IRI—2007電子密度“真值”比較Fig.2 Comparison of the ionospheric electron density profiles estimated from different GNSS tomographic algorithms using the simulation observations with the ionospheric electron density profiles from the IRI—2007 model

    方案二IART算法通過引入層析格網電子密度像素內射線截距與電子密度乘積為組合變量,較為合理地分配不同電子密度像素格網內實測值與反演值之間的誤差,同時通過構造與電子密度變化相關的迭代松弛因子,抑制傳播噪聲對電子密度反演結果影響,降低了電子密度反演誤差。從圖2可以看出,方案二提供的以藍色實線與加號組合代表的電子密度結果(尤其是在高度200~300 km)更靠近IRI—2007模型提供的電子密度“真值”(黑色實線),優(yōu)于ART算法結果(綠色虛線)。在方案二的基礎上,方案三和方案四在電離層高度上采用了不等間距層析格網劃分方法,充分考慮了電子密度在垂直高度方向上變化較大的特性,在峰值電子密度區(qū)域劃分較小的像素格網而其他區(qū)域劃分較大像素格網,合理調節(jié)不同高度層析格網內射線截距對電子密度誤差的放大作用,提高電子密度反演結果精度??梢钥吹剑瑘D2中方案三提供的電子密度(青色實線與方塊表示)和方案四提供的電子密度(紅色實線與三角形表示)優(yōu)于前兩個方案的電子密度結果。不同于方案三在地理坐標系下劃分層析電子密度像素格網,顧及電離層受地球電磁場影響,方案四提出在地磁坐標系下劃分電子密度像素格網,像素類層析算法假定了每個像素格網內的電子密度均勻分布,同時兼顧實測數據分布稀疏,往往在實際反演計算中沿著經度與緯度對不同電子密度像素格網施加約束,因此在地磁坐標系下的電子密度約束反演結果將導致和在地理坐標系下的電子密度約束反演結果不同。圖2展示的結果反映出方案四在地磁坐標系下的電子密度結果優(yōu)于方案三地理坐標系下的電子密度結果,也是四種電離層層析算法反演方案中結果最優(yōu)的方案。

    2)電離層電子密度反演誤差與精度

    表1與表2定量給出12月1日、12月2日、12月6日、12月7日四種不同GNSS電離層層析反演方案計算的電子密度在高度上的剖面結果相對于IRI—2007電子密度“真值”的峰值電子密度絕對誤差、電子密度剖面平均絕對百分比誤差及均方根值(Root Mean Square, RMS)。其中,表1給出的是北京上空的平均統(tǒng)計結果,表2給出的是武漢上空的平均統(tǒng)計結果。需要補充說明的是,底部電離層電子密度值和頂部電離層電子密度值相對峰值電子密度較小,部分底部與頂部電子密度反演結果中的誤差變化可能會導致電子密度剖面誤差百分比的統(tǒng)計結果出現較大波動。為了更準確地表達不同GNSS電離層層析方法反演電子密度誤差的總體效果,本文在模擬實驗中僅統(tǒng)計電離層高度處于200~450 km處的電子密度結果百分比誤差。

    根據表1~2統(tǒng)計結果可以看出,相對于IRI—2007模型提供的電子密度“真值”,方案一ART算法反演的峰值電子密度誤差、電離層電子密度剖面結果的平均絕對百分比誤差及RMS值,其結果都最大;方案二電離層電子密度反演誤差與RMS值其次;方案三電離層電子密度反演誤差與RMS值總體優(yōu)于方案二結果,方案四的電離層電子密度反演誤差與RMS值最小。上述研究結果表明,本文通過改進ART算法、建立不等間距電子密度像素格網及在地磁坐標系下反演電子密度,都不同程度地提高了GNSS電離層層析反演電子密度結果的可靠性。另外,從表1~2的統(tǒng)計結果可以看出,不同方案反演的電離層峰值電子密度誤差和電子密度百分比統(tǒng)計誤差的結果差異最大,表明不同方案在靠近電子密度較大的峰值區(qū)域結果改進最明顯,這和圖2展示的結果一致;電子密度剖面結果的RMS誤差則反映了不同方案在電離層全部高度區(qū)域的整體改進效果。綜上所述,在四個不同的電離層層析反演方案中,方案四反演的電子密度結果均是最優(yōu)的。

    表1 不同GNSS電離層層析算法模擬重構的北京電子密度剖面相對于IRI—2007“真值”的誤差與精度統(tǒng)計Tab.1 Error and accuracy statistics of the ionospheric electron density derived from the different GNSS tomographic algorithms using the simulation observations in comparison with those provided by the IRI—2007 model over Beijing during the different days

    表2 不同GNSS電離層層析算法模擬重構的武漢電子密度剖面相對于IRI—2007“真值”的誤差與精度統(tǒng)計Tab.2 Error and accuracy statistics of the ionospheric electron density derived from the different GNSS tomographic algorithms using the simulation observations in comparison with those provided by the IRI—2007 model over Wuhan during the different days

    4.2 實測數據比較研究

    利用2011年12月1日、12月2日、12月6日、12月7日的中國陸態(tài)網絡145個GPS地面監(jiān)測站的實測數據開展電離層電子密度層析算法的實測數據試驗比較研究。和4.1小節(jié)模擬實驗一致,實測數據比較討論時采用了四種不同的電離層層析反演方案。同時,為比較評估GPS實測數據層析反演的電離層電子密度精度,實驗過程中采用了北京站和武漢站電離層測高儀觀測獲得的電子密度觀測結果作為“真值”進行比較論證。需要說明的是,電離層測高儀能夠通過直接觀測獲取峰值電子密度及其高度以下底部電子密度的準確結果,電子密度峰值高度以上的頂部電子密度卻是根據一定的算法獲得[32],因此,本小節(jié)在比較不同方法獲得的頂部電子密度時,僅進行定性的趨勢變化比較而不開展定量的結果差異比較,對底部電離層電子密度及峰值電子密度進行了定量分析與討論。

    1)電離層電子密度高度剖面反演結果比較

    圖3給出了北京和武漢上空在2011年12月6日北京時間11:00、13:00和15:00點電離層電子密度在高度方向上的剖面變化結果。其他時段結果類似,其電離層電子密度反演誤差與精度統(tǒng)計結果見表3~4。

    在圖3中電離層電子密度“真值”由電離層測高儀觀測提供(黑色實線表示),其他四種反演方案得到的電離層電子密度剖面結果表達與圖2保持一致。從圖3可以看出,方案四反演的電離層電子密度結果最接近于電離層測高儀觀測獲得的電子密度剖面,尤其是在靠近電子密度峰值區(qū)域的結果具有明顯改善效果,這與4.1.2小節(jié)利用IRI—2007模型開展的模擬實驗結果保持一致;但也可以看出方案四在部分時段的反演結果與方案三結果接近,這和在模擬實驗中方案四相對其他方案都具有明顯優(yōu)勢不同,這可能與實際電離層變化活動除受地球電磁場影響之外還受中性風等其他因素影響,從而導致其更為復雜多變有關。

    (a) 2011年12月6日北京電離層電子密度結果(a) Ionospheric electron density profiles over Beijing on December 6, 2011

    (b) 2011年12月6日武漢電離層電子密度結果(b) Ionospheric electron density over Wuhan on December 6, 2011圖3 不同GNSS層析算法重構的電子密度與測高儀電子密度“真值”比較結果Fig.3 Comparison of the ionospheric electron density profiles derived from the different GNSS tomographic algorithms using GNSS data with the ionospheric electron density profiles from the ionosonde measurements

    表3 不同GNSS電離層層析算法重構的北京電子密度相對于電離層測高儀觀測“真值”的誤差與精度統(tǒng)計Tab.3 Error and accuracy statistics of the ionospheric electron density derived from the different GNSS tomographic algorithms using GNSS data in comparison with those from the ionosonde measurements over Beijing

    表4 不同GNSS電離層層析算法重構的武漢電子密度相對于電離層測高儀觀測“真值”的誤差與精度統(tǒng)計

    圖3結果也顯示,GNSS電離層層析系統(tǒng)采用不等間距電子密度像素格網劃分方式后,方案三反演的電離層電子密度結果優(yōu)于方案二等間距電子密度像素格網算法反演的電子密度,傳統(tǒng)ART算法反演的電離層電子密度效果最差,這和模擬實驗結果保持一致。

    2)電離層電子密度反演誤差與精度

    進一步以電離層測高儀觀測提供的電子密度為“真值”,定量統(tǒng)計四種不同的GNSS電離層層析方案反演的峰值電子密度誤差的絕對值、底部電子密度剖面平均絕對百分比誤差及均方根值。同樣需要補充說明的是,為避免部分底部電離層電子密度反演結果變化可能導致的電子密度剖面誤差百分比統(tǒng)計出現較大波動,本文在實測數據比較研究中僅統(tǒng)計電離層高度處于200 km至峰值電子密度高度處結果的百分比誤差,表3與表4分別給出北京和武漢上空不同層析反演方案在12月1日、12月2日、12月6日、12月7日的電子密度結果誤差與精度的平均統(tǒng)計結果。

    從表3~4可以看到,相對于電離層測高儀觀測獲得的電離層電子密度“真值”,方案四反演的電離層峰值電子密度誤差、底部電離層電子密度剖面結果平均絕對百分比誤差以及RMS統(tǒng)計結果基本最?。环桨敢环囱蓦娮用芏冉Y果的誤差與精度最大;方案二和方案三反演的電子密度結果均優(yōu)于方案一的結果,方案三結果優(yōu)于方案二結果(12月7日武漢站結果除外)。從峰值電子密度誤差和平均百分比誤差看,不同方案反演的結果差異最為明顯。上述實測數據研究結果和模擬實驗結果基本保持一致。此外,從總體反演效果看,不同方案反演的武漢上空電離層電子密度誤差與精度大于北京上空的電子密度誤差與精度,這或許與我國低緯度區(qū)域電離層實際物理變化活動比中緯度區(qū)域更為復雜多變有關。

    最后需要補充說明的是,受太陽活動等多種因素影響,電離層空間環(huán)境變化復雜,存在逐日、月、半年、年甚至更長時間尺度(11年)的變化活動規(guī)律。因此,本文在地磁坐標系下建立不等像素間距的GNSS層析反演方法的討論與分析結果僅是初步的,更深入的分析與論證需要利用長時間尺度的實測數據進行研究。

    5 結論

    本文初步探討了GNSS電離層電子密度層析技術四種不同的反演方案,包括傳統(tǒng)的ART算法、IART算法、在電離層高度上建立電子密度像素不等間距格網的IART算法及在地磁坐標系下建立電子密度像素不等間距格網的IART算法。利用IRI—2007電離層模型、GNSS實測數據及電離層測高儀觀測數據,從模擬實驗和實測數據對比兩個方面,定性分析與定量討論了不同GNSS電離層層析方案的電子密度反演誤差與精度,主要結論如下:

    1)相對于ART電離層電子密度層析算法,本文在引入觀測射線截距與電子密度乘積組合為層析迭代自變量的基礎上,通過構造與電子密度變化相關的松弛因子實現平衡與調節(jié)層析反演的電子密度精度與結果平滑程度,降低電子密度反演誤差并提高反演精度;

    2)GNSS電離層層析方法等間距劃分電子密度像素格網,忽略了電離層電子密度在垂直高度方向上變化起伏較大的不均勻性,本文提出在接近電子密度變化較大的高度區(qū)域建立較小的電子密度像素格網而其他高度區(qū)域建立較大的電子密度像素格網,較為合理地提高電離層層析系統(tǒng)的空間分辨率,明顯改善靠近電子密度峰值區(qū)域反演效果;

    3)考慮電離層受地球電磁場影響較大,提出在地磁坐標系下建立電離層電子密度層析格網,使得能夠在電離層層析系統(tǒng)中沿著空間經緯度施加更為合理的電子密度平滑約束,提高電離層電子密度反演精度。

    猜你喜歡
    電子密度格網層析
    犬細小病毒量子點免疫層析試紙條的研制
    實時電離層格網數據精度評估
    不同GPS掩星電離層剖面產品相關性分析
    測繪通報(2019年11期)2019-12-03 01:47:34
    等離子體電子密度分布信息提取方法研究
    飛秒激光脈沖成絲的診斷及其應用概述
    一種適用于電離層電子密度重構的AMART算法
    測繪學報(2018年1期)2018-02-27 02:23:07
    基于空間信息格網與BP神經網絡的災損快速評估系統(tǒng)
    A族鏈球菌膠體金免疫層析試紙條的制備及應用
    新型B族鏈球菌膠體金免疫層析試紙條的臨床應用評價
    平均Helmert空間重力異常格網構制方法
    少妇高潮的动态图| 伦理电影免费视频| 高清视频免费观看一区二区| 日日爽夜夜爽网站| 狂野欧美激情性xxxx在线观看| 国产精品一区二区在线不卡| 亚洲国产欧美在线一区| 日本午夜av视频| 午夜福利网站1000一区二区三区| 婷婷色麻豆天堂久久| av.在线天堂| 精品人妻熟女毛片av久久网站| 91久久精品国产一区二区三区| 九草在线视频观看| 午夜免费观看性视频| 成年女人在线观看亚洲视频| 国产亚洲av片在线观看秒播厂| 99热这里只有是精品在线观看| 51国产日韩欧美| 国产精品99久久99久久久不卡 | 天堂俺去俺来也www色官网| 日韩视频在线欧美| 纵有疾风起免费观看全集完整版| 国国产精品蜜臀av免费| 欧美精品亚洲一区二区| 国产午夜精品一二区理论片| 人成视频在线观看免费观看| 国国产精品蜜臀av免费| 建设人人有责人人尽责人人享有的| 秋霞伦理黄片| av在线观看视频网站免费| 亚洲精品乱久久久久久| 久久综合国产亚洲精品| 91在线精品国自产拍蜜月| 只有这里有精品99| 国产 一区精品| 亚洲欧洲日产国产| 99久国产av精品国产电影| 热99国产精品久久久久久7| 超碰97精品在线观看| 亚洲美女搞黄在线观看| 久久久欧美国产精品| 亚洲中文av在线| 三级国产精品欧美在线观看| 日韩不卡一区二区三区视频在线| 有码 亚洲区| 性色avwww在线观看| 中文字幕亚洲精品专区| 亚洲,欧美,日韩| av网站免费在线观看视频| 成人黄色视频免费在线看| 草草在线视频免费看| 国产精品人妻久久久影院| 大香蕉久久成人网| 九九久久精品国产亚洲av麻豆| 亚洲国产av影院在线观看| 大码成人一级视频| 色哟哟·www| 欧美性感艳星| 亚洲色图综合在线观看| 亚洲国产精品国产精品| 美女脱内裤让男人舔精品视频| 久久久久久久久大av| 成人无遮挡网站| 内地一区二区视频在线| 日韩熟女老妇一区二区性免费视频| 26uuu在线亚洲综合色| 亚洲欧美日韩卡通动漫| 涩涩av久久男人的天堂| 亚洲av欧美aⅴ国产| 高清av免费在线| 美女国产高潮福利片在线看| 亚洲精华国产精华液的使用体验| 狠狠婷婷综合久久久久久88av| 精品一区二区免费观看| 色婷婷久久久亚洲欧美| 大又大粗又爽又黄少妇毛片口| 久久人妻熟女aⅴ| 国产精品国产三级国产av玫瑰| 国语对白做爰xxxⅹ性视频网站| 91成人精品电影| 九九在线视频观看精品| 日本猛色少妇xxxxx猛交久久| 美女国产高潮福利片在线看| 国产午夜精品一二区理论片| 欧美bdsm另类| 天堂8中文在线网| 啦啦啦中文免费视频观看日本| 亚洲av二区三区四区| 欧美日韩一区二区视频在线观看视频在线| 一级毛片黄色毛片免费观看视频| 交换朋友夫妻互换小说| 欧美日韩视频精品一区| 日日撸夜夜添| 精品一区在线观看国产| 黄色毛片三级朝国网站| 在线 av 中文字幕| 欧美 日韩 精品 国产| 国产高清三级在线| 国产一区二区三区综合在线观看 | 三级国产精品欧美在线观看| 亚洲av男天堂| 午夜免费鲁丝| 精品国产国语对白av| 女性被躁到高潮视频| 伊人久久国产一区二区| 国产在线视频一区二区| 亚洲精品乱码久久久久久按摩| 最近最新中文字幕免费大全7| 啦啦啦中文免费视频观看日本| 97精品久久久久久久久久精品| 一个人免费看片子| 国产欧美亚洲国产| 美女大奶头黄色视频| 国产成人91sexporn| 亚洲av电影在线观看一区二区三区| 在线观看免费高清a一片| 热99国产精品久久久久久7| 亚洲国产成人一精品久久久| 黄色欧美视频在线观看| 欧美人与善性xxx| 91久久精品电影网| 国产欧美日韩综合在线一区二区| 黑人猛操日本美女一级片| 午夜老司机福利剧场| 麻豆精品久久久久久蜜桃| 免费人成在线观看视频色| 精品国产一区二区三区久久久樱花| 欧美日韩在线观看h| 国产精品三级大全| 国产在视频线精品| 国产精品一区www在线观看| 亚洲国产精品国产精品| 免费大片黄手机在线观看| freevideosex欧美| 成人国语在线视频| 热re99久久国产66热| 啦啦啦啦在线视频资源| 国产精品国产av在线观看| 日本-黄色视频高清免费观看| 精品国产一区二区久久| 久久久久久久精品精品| 亚洲欧美日韩另类电影网站| 欧美一级a爱片免费观看看| 国产精品一区二区在线观看99| 男人操女人黄网站| 18禁在线播放成人免费| 国产精品一区www在线观看| 亚洲第一av免费看| 国产av码专区亚洲av| 又大又黄又爽视频免费| xxxhd国产人妻xxx| 国产精品一二三区在线看| 80岁老熟妇乱子伦牲交| 大话2 男鬼变身卡| 久久久久视频综合| 午夜老司机福利剧场| 亚洲美女搞黄在线观看| 岛国毛片在线播放| 成人黄色视频免费在线看| 免费大片18禁| 99热这里只有是精品在线观看| 插阴视频在线观看视频| 亚洲av福利一区| av一本久久久久| 男人操女人黄网站| 伦理电影免费视频| 亚洲人与动物交配视频| 午夜免费观看性视频| 中文字幕精品免费在线观看视频 | 国产乱来视频区| 一区二区日韩欧美中文字幕 | 亚洲中文av在线| 日韩伦理黄色片| 国产成人av激情在线播放 | 波野结衣二区三区在线| 日韩不卡一区二区三区视频在线| 天美传媒精品一区二区| 亚洲精品第二区| 免费黄频网站在线观看国产| 久久久久国产网址| 国产极品天堂在线| 在线观看三级黄色| 精品少妇内射三级| 久久精品久久久久久久性| 日韩欧美一区视频在线观看| 美女脱内裤让男人舔精品视频| 哪个播放器可以免费观看大片| 夫妻午夜视频| 亚洲欧美中文字幕日韩二区| 少妇精品久久久久久久| 久久鲁丝午夜福利片| 国产熟女午夜一区二区三区 | 亚洲av二区三区四区| 九九久久精品国产亚洲av麻豆| av国产久精品久网站免费入址| 欧美成人午夜免费资源| av视频免费观看在线观看| 久久99一区二区三区| 一级片'在线观看视频| 国内精品宾馆在线| 丝袜喷水一区| 你懂的网址亚洲精品在线观看| 十分钟在线观看高清视频www| 大片免费播放器 马上看| 欧美成人午夜免费资源| 国产精品一区二区在线观看99| 九九久久精品国产亚洲av麻豆| 美女内射精品一级片tv| 欧美 日韩 精品 国产| 人人妻人人爽人人添夜夜欢视频| 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 成人二区视频| 看十八女毛片水多多多| 日本vs欧美在线观看视频| 插阴视频在线观看视频| 天天影视国产精品| 两个人的视频大全免费| 一级毛片电影观看| 国产视频内射| 十八禁网站网址无遮挡| 国产一区二区在线观看日韩| 久久韩国三级中文字幕| 黑人巨大精品欧美一区二区蜜桃 | 这个男人来自地球电影免费观看 | 亚洲精品乱码久久久v下载方式| 国产在线免费精品| 国产精品秋霞免费鲁丝片| 美女内射精品一级片tv| 十八禁高潮呻吟视频| 免费观看av网站的网址| 成人国语在线视频| 久久国产亚洲av麻豆专区| 中文天堂在线官网| 人人妻人人澡人人看| 天美传媒精品一区二区| 99热6这里只有精品| 久久这里有精品视频免费| 免费观看在线日韩| 日韩不卡一区二区三区视频在线| 亚洲欧美日韩卡通动漫| 亚洲少妇的诱惑av| 我的老师免费观看完整版| 欧美少妇被猛烈插入视频| 日本黄色片子视频| 午夜福利视频在线观看免费| 最近中文字幕2019免费版| 大码成人一级视频| 国产黄色视频一区二区在线观看| 中国美白少妇内射xxxbb| 日韩中字成人| videossex国产| 日韩视频在线欧美| 2022亚洲国产成人精品| 国产黄片视频在线免费观看| 自线自在国产av| 国产精品久久久久久久电影| 国产成人精品婷婷| 一边摸一边做爽爽视频免费| 99国产精品免费福利视频| 日韩不卡一区二区三区视频在线| 亚洲精品国产av成人精品| 秋霞伦理黄片| 午夜av观看不卡| 91精品伊人久久大香线蕉| 丁香六月天网| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久 | 草草在线视频免费看| 久久久精品94久久精品| 妹子高潮喷水视频| 国产有黄有色有爽视频| 精品久久国产蜜桃| 嘟嘟电影网在线观看| videosex国产| 最近最新中文字幕免费大全7| 丝袜脚勾引网站| 高清视频免费观看一区二区| 满18在线观看网站| 亚洲精品一二三| 黄片无遮挡物在线观看| 久久ye,这里只有精品| 国产亚洲一区二区精品| 大陆偷拍与自拍| 18禁观看日本| 18禁在线播放成人免费| 全区人妻精品视频| 视频在线观看一区二区三区| 黄片无遮挡物在线观看| 91午夜精品亚洲一区二区三区| 婷婷色综合大香蕉| 中文精品一卡2卡3卡4更新| 99久久综合免费| 老熟女久久久| 大香蕉久久网| 中国三级夫妇交换| 免费观看的影片在线观看| 超碰97精品在线观看| 精品亚洲成国产av| 欧美精品人与动牲交sv欧美| 免费播放大片免费观看视频在线观看| 免费播放大片免费观看视频在线观看| 高清av免费在线| 观看美女的网站| 2021少妇久久久久久久久久久| 最近2019中文字幕mv第一页| 国产欧美日韩一区二区三区在线 | freevideosex欧美| 一区在线观看完整版| 国产一区亚洲一区在线观看| 免费大片18禁| 日韩免费高清中文字幕av| 青春草国产在线视频| 国产乱人偷精品视频| 精品人妻在线不人妻| 国精品久久久久久国模美| 亚洲欧美日韩卡通动漫| 女性被躁到高潮视频| 亚洲精品国产av蜜桃| 日韩精品免费视频一区二区三区 | 一本一本综合久久| 久久99热6这里只有精品| 精品亚洲乱码少妇综合久久| 少妇被粗大的猛进出69影院 | 丝袜脚勾引网站| 久久久久久人妻| 国产精品国产av在线观看| 色94色欧美一区二区| 汤姆久久久久久久影院中文字幕| 国产精品秋霞免费鲁丝片| 国国产精品蜜臀av免费| 美女福利国产在线| 久热久热在线精品观看| 一级毛片黄色毛片免费观看视频| 大片电影免费在线观看免费| 久久精品国产亚洲网站| videosex国产| √禁漫天堂资源中文www| 久久人人爽人人片av| 精品人妻熟女毛片av久久网站| 亚洲伊人久久精品综合| 国产成人a∨麻豆精品| 欧美成人午夜免费资源| 成年人免费黄色播放视频| 日本av免费视频播放| 久久影院123| 亚洲精品成人av观看孕妇| 精品久久久噜噜| 精品国产一区二区三区久久久樱花| 九九久久精品国产亚洲av麻豆| 高清不卡的av网站| 九九爱精品视频在线观看| 女人久久www免费人成看片| 国产伦精品一区二区三区视频9| av在线播放精品| 蜜桃国产av成人99| 99国产精品免费福利视频| 国产精品无大码| 日本爱情动作片www.在线观看| 蜜桃国产av成人99| 亚洲国产日韩一区二区| 国产精品一国产av| 日本爱情动作片www.在线观看| 国产精品久久久久久久久免| 99re6热这里在线精品视频| 制服诱惑二区| 久久青草综合色| 亚洲精品国产av蜜桃| 亚洲成色77777| 久久久久久伊人网av| 人妻人人澡人人爽人人| 免费黄频网站在线观看国产| 亚洲av二区三区四区| 少妇人妻久久综合中文| 日韩欧美精品免费久久| 亚洲丝袜综合中文字幕| 这个男人来自地球电影免费观看 | 最近中文字幕高清免费大全6| 91精品国产九色| 插逼视频在线观看| 亚洲av电影在线观看一区二区三区| 久久久久久久大尺度免费视频| 免费观看在线日韩| 久久综合国产亚洲精品| 午夜福利视频在线观看免费| 亚洲精品亚洲一区二区| 好男人视频免费观看在线| 亚洲三级黄色毛片| 亚洲精品aⅴ在线观看| 夫妻午夜视频| 51国产日韩欧美| 国产不卡av网站在线观看| 日本与韩国留学比较| 免费观看性生交大片5| 亚洲精品456在线播放app| 亚洲中文av在线| 一级毛片 在线播放| 美女内射精品一级片tv| 亚洲欧美色中文字幕在线| 人成视频在线观看免费观看| 欧美另类一区| 婷婷成人精品国产| 亚洲av二区三区四区| 久久久久视频综合| 日韩 亚洲 欧美在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 超色免费av| 亚洲欧美日韩卡通动漫| 97在线人人人人妻| 日本-黄色视频高清免费观看| 午夜激情久久久久久久| 成年人午夜在线观看视频| 在线精品无人区一区二区三| 日韩三级伦理在线观看| 韩国高清视频一区二区三区| 亚洲国产日韩一区二区| 黑人高潮一二区| 最近的中文字幕免费完整| 日本黄色日本黄色录像| 精品久久久久久电影网| 99久久精品国产国产毛片| 一区二区日韩欧美中文字幕 | 欧美97在线视频| 亚洲精品视频女| 91在线精品国自产拍蜜月| 国产视频内射| av视频免费观看在线观看| 黄片无遮挡物在线观看| 91精品伊人久久大香线蕉| 日本91视频免费播放| 男的添女的下面高潮视频| 国产日韩一区二区三区精品不卡 | 日韩成人伦理影院| 亚洲av成人精品一二三区| 99久国产av精品国产电影| 欧美变态另类bdsm刘玥| 一区在线观看完整版| 午夜老司机福利剧场| 青春草视频在线免费观看| 蜜臀久久99精品久久宅男| 老司机影院成人| 两个人免费观看高清视频| 哪个播放器可以免费观看大片| 一二三四中文在线观看免费高清| 在线观看免费视频网站a站| 国产 一区精品| 丰满少妇做爰视频| av视频免费观看在线观看| 精品国产乱码久久久久久小说| 欧美激情 高清一区二区三区| 欧美xxⅹ黑人| 国产日韩一区二区三区精品不卡 | 亚洲国产精品999| 亚洲av男天堂| 51国产日韩欧美| 国产白丝娇喘喷水9色精品| 色视频在线一区二区三区| 99久久人妻综合| 国产免费福利视频在线观看| 激情五月婷婷亚洲| 女的被弄到高潮叫床怎么办| 久久女婷五月综合色啪小说| 国产男人的电影天堂91| 51国产日韩欧美| 伊人久久精品亚洲午夜| 日韩制服骚丝袜av| 99久久精品国产国产毛片| 国产午夜精品一二区理论片| 亚洲美女黄色视频免费看| 久久午夜福利片| 在线观看免费视频网站a站| 国产午夜精品久久久久久一区二区三区| 日日摸夜夜添夜夜爱| 国产免费又黄又爽又色| 久久99精品国语久久久| 毛片一级片免费看久久久久| 大又大粗又爽又黄少妇毛片口| 男女啪啪激烈高潮av片| 国精品久久久久久国模美| 男人添女人高潮全过程视频| 国产黄色免费在线视频| 男男h啪啪无遮挡| 大又大粗又爽又黄少妇毛片口| 国产精品偷伦视频观看了| videossex国产| 秋霞伦理黄片| 少妇丰满av| 韩国av在线不卡| 夫妻性生交免费视频一级片| 观看美女的网站| 久久这里有精品视频免费| 久久久久久伊人网av| 国产精品蜜桃在线观看| 美女大奶头黄色视频| 在线观看免费高清a一片| 国产精品人妻久久久久久| 日韩大片免费观看网站| 国产一级毛片在线| 国产永久视频网站| 插阴视频在线观看视频| 国产探花极品一区二区| 国产精品人妻久久久影院| 亚洲av不卡在线观看| 国产精品国产三级专区第一集| av电影中文网址| 久久鲁丝午夜福利片| 国产男女内射视频| 晚上一个人看的免费电影| 久久免费观看电影| 日本爱情动作片www.在线观看| 中文字幕亚洲精品专区| 最新的欧美精品一区二区| 亚洲精品日本国产第一区| 国产精品一区二区在线观看99| 黄片播放在线免费| 男女无遮挡免费网站观看| 菩萨蛮人人尽说江南好唐韦庄| 国产69精品久久久久777片| 美女大奶头黄色视频| 国产在视频线精品| 女人久久www免费人成看片| 亚洲情色 制服丝袜| 久久婷婷青草| 女人精品久久久久毛片| 国产精品久久久久久久久免| 国产精品 国内视频| 大香蕉97超碰在线| 国产白丝娇喘喷水9色精品| 黄色怎么调成土黄色| 热99久久久久精品小说推荐| 在线观看三级黄色| av黄色大香蕉| 日韩精品有码人妻一区| 亚洲欧美日韩另类电影网站| 午夜福利在线观看免费完整高清在| 少妇猛男粗大的猛烈进出视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丰满迷人的少妇在线观看| 久热这里只有精品99| 777米奇影视久久| 青春草视频在线免费观看| 少妇的逼好多水| 成人午夜精彩视频在线观看| 国产成人精品婷婷| 成人综合一区亚洲| 成人国语在线视频| 日本91视频免费播放| 免费不卡的大黄色大毛片视频在线观看| 又黄又爽又刺激的免费视频.| 日日撸夜夜添| 丝袜喷水一区| 国产成人精品婷婷| 中文字幕久久专区| 国产男女内射视频| 亚洲精品久久成人aⅴ小说 | 亚洲,一卡二卡三卡| 99国产精品免费福利视频| 亚洲av电影在线观看一区二区三区| 两个人免费观看高清视频| 国产亚洲欧美精品永久| 亚洲第一av免费看| 黑人高潮一二区| 免费黄频网站在线观看国产| 久久久久久久久大av| 欧美成人午夜免费资源| 又大又黄又爽视频免费| 飞空精品影院首页| 哪个播放器可以免费观看大片| 热re99久久精品国产66热6| 久久99精品国语久久久| 少妇 在线观看| 亚洲av中文av极速乱| 国产黄色视频一区二区在线观看| 母亲3免费完整高清在线观看 | 成年女人在线观看亚洲视频| 国产男人的电影天堂91| 日本欧美视频一区| 精品午夜福利在线看| tube8黄色片| 久久精品人人爽人人爽视色| 欧美性感艳星| 亚洲国产精品999| 综合色丁香网| 国产在线视频一区二区| 国产av国产精品国产| 大香蕉久久网| 久久韩国三级中文字幕| 五月开心婷婷网| 热re99久久国产66热| 久久国产精品大桥未久av| 精品久久久精品久久久| 高清午夜精品一区二区三区| 欧美 亚洲 国产 日韩一| 国产免费福利视频在线观看| 五月天丁香电影| 成人毛片60女人毛片免费| 丝袜脚勾引网站| 亚洲精品日韩在线中文字幕| 国产精品国产三级国产专区5o| 久久99蜜桃精品久久| 美女视频免费永久观看网站| 日韩中字成人| 亚洲精品乱码久久久v下载方式| 最黄视频免费看| 嘟嘟电影网在线观看| 五月伊人婷婷丁香| 亚洲精品乱码久久久v下载方式| 亚洲国产精品专区欧美| 国产亚洲精品第一综合不卡 | 少妇人妻久久综合中文| 精品亚洲成国产av| 啦啦啦在线观看免费高清www| 日韩一区二区视频免费看| 久久久精品94久久精品|