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

    利用轉移熵研究引起磁暴擾動的太陽風參數(shù)重要性排序*

    2022-06-20 09:03:54于佳斌佟繼周方少峰胡曉彥
    空間科學學報 2022年3期
    關鍵詞:磁層磁暴太陽活動

    于佳斌 佟繼周 方少峰 胡曉彥

    1(中國科學院國家空間科學中心 北京 100190)

    2(中國科學院大學 北京 100049)

    0 引言

    太陽風–磁層相互作用形成了一個多尺度耦合的復雜系統(tǒng),太陽風是引起地磁暴的最主要驅動源,例如當行星際磁場南向時,行星際磁場可與地球磁場發(fā)生磁場重聯(lián),較易引發(fā)磁暴。因而,理解太陽風變化對地磁場的影響機制,是預報地磁暴的關鍵之一。

    有關太陽風與地磁擾動之間的相關性研究,新方法不斷出現(xiàn),這些方法主要基于太陽風和地磁參數(shù)的解析關系、相關系數(shù)或預報模型,線性回歸、統(tǒng)計相關等方法,已被證實是認識磁暴時地磁變化的有效方法,例如Liu 等[1]研究了太陽活動上升年的地磁Kp指數(shù)與行星際磁場南向分量Bz、太陽風速度vsw、太陽風溫度Tsw與太陽風數(shù)密度Dsw的相關性,相關系數(shù)分別為0.66、0.03、0.58 和–0.45,發(fā)現(xiàn)有一類磁暴與Bz地磁場重聯(lián)無關。Cane 等[2]統(tǒng)計結果表明行星際磁場南向分量Bz與 磁暴擾動指數(shù)最小值Dstmin的相關系數(shù)為0.74;Wu[3]統(tǒng)計分析了1995-1998 年Wind 數(shù)據(jù),證實vswBzmin與Dstmin的相關系數(shù)為0.87,略高于行星際磁場南向分量最小值Bzmin與磁暴擾動指數(shù)最小值Dstmin的相關系數(shù)0.81。Zhang 等[4]發(fā)現(xiàn)行星際電場E與Dst指數(shù)有很好的相關性,并且在中、強等地磁活動期間,存在顯著的突變特征曲線,相對于vsw,,vswB2和 能量耦合函數(shù)ε=,行星際電場E的突變特征曲線更易識別。Zhao 等[5]的研究表明在諸多參數(shù)中,就單一因素來說,行星際電場E對磁暴強度影響最大,行星際磁場南向分量Bz對磁暴強度影響次之;Khabarova等[6]和Ahmed 等[7]也得到了相似的結果。除此之外,Iyemori 等[8]首次成功地應用線性預測濾波方法,基于太陽風參數(shù)預測了地磁活動,Ji 等[9]和Rast?tter等[10]用經(jīng)驗和物理模型、Ahmed 等[7]用神經(jīng)網(wǎng)絡對磁暴進行預報研究。然而,地磁系統(tǒng)對太陽風變化的響應是非線性的,線性統(tǒng)計不足以表達其完備的相關性,同時在使用各種預報模型時,不同的輸入?yún)?shù)對磁暴的影響和重要性仍不明確,關于哪個太陽風參數(shù)對磁暴的發(fā)生更重要,還未有統(tǒng)一的認識。

    近年來,基于信息論的新方法[11]被提出并成功地應用于揭示地球磁層–電離層系統(tǒng)響應太陽風變化的復雜動力學特征。Michelis 等[12]以地磁活動指數(shù)AL和Sym-H分別定量化亞暴和磁暴過程,進行了雙變量轉移熵分析,結果表明亞暴和磁暴之間的信息流動方向取決于全球地磁活動水平。Wing 等[13]利用信息論進行太陽風–地球輻射帶系統(tǒng)研究,利用互信息、條件互信息和轉移熵對能量范圍為1.8~3.5 MeV的地球同步軌道高能電子通量的太陽風驅動因素進行了研究,得出了基于信息論分析的重要性排序。轉移熵可以提供強大的無模型統(tǒng)計量,在太陽風–地球輻射帶系統(tǒng)研究已經(jīng)得到了驗證,但目前還未在太陽風–磁層整體耦合系統(tǒng)中得到應用。本文利用轉移熵研究太陽風參數(shù)與地磁指數(shù)之間的非線性關系,應用轉移熵來衡量磁暴時的太陽風驅動因素并進行量化,對不同太陽風參數(shù)的重要性進行研究。

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)來源

    用于研究的太陽風及地磁活動指數(shù)的數(shù)據(jù)均來源于美國國家航空航天局(NASA)的OMNI*http://omniweb.gsfc.nasa.gov/公開數(shù)據(jù),具體包括分鐘分辨率數(shù)據(jù)(HRO)和小時分辨率數(shù)據(jù)(LRO),時間跨度從1996 年到2018 年。OMNI 整合了來自ACE、WIND、IMP8 等航天器的太陽風磁場和等離子體數(shù)據(jù)集,并根據(jù)實際的太陽風速度時移到地球的弓激波鼻尖處。此外,本研究還使用了世界地磁數(shù)據(jù)中心(WDC)提供的地磁Sym-H指數(shù)數(shù)據(jù)、比利時SILSO 中心提供的太陽黑子數(shù)據(jù)、和加拿大自然資源部的太陽10.7 cm 射電流量F10.7數(shù)據(jù)。

    本研究利用美國空間天氣研究中心(SWRC)提供的空間天氣數(shù)據(jù)庫(DONKI)**https://kauai.ccmc.gsfc.nasa.gov/DONKI/中提取的磁暴事件列表,共計93 個事例,主要分布在第24 太陽活動周期間的2010 年到2018 年。

    1.2 相關系數(shù)與轉移熵

    1.2.1 相關系數(shù)算法

    相關系數(shù)又稱皮爾遜積矩相關系數(shù),適用于度量兩個變量x與y之間的線性相關性,其值介于1 與–1之間,即

    直接推導可知

    時移相關系數(shù)[13]是皮爾遜積矩相關系數(shù)度量向兩個時間序列x(t)和y(t+τ)的拓展,即

    其中τ為時間延遲。當相關函數(shù)有多個峰或沒有明顯的不對稱時,這種分析的結果可能不是特別清楚。此外,相關系數(shù)是兩個變量之間線性關系的一個度量,不能用于描述非線性關系。

    1.2.2 轉移熵算法

    度量兩個變量x和y之間依賴關系,更通用的方法是通過考慮

    其中,P(x,y)為x和y的聯(lián)合概率。如果式(2)等號成立,則表示變量x與y相互獨立;相比式(1),除線性關系外,式(2)可以表達更多的非線性關系。

    信息論中互信息[14]表示一個隨機變量中包含的關于另一個隨機變量的信息量,代表兩個隨機變量的依賴程度。

    對于離散的概率分布,互信息的定義如下:

    從式(3)和(4)可以看出,互信息M(x,y)對式(2)進行了定量描述,相比于相關系數(shù),更好地表達了兩個變量x與y之間的依賴關系。進一步,如果x與y服從聯(lián)合高斯概率分布,互信息具有如下表達式[14]:

    互信息的取值范圍為[ 0,∞),經(jīng)過轉化Λ(x,y)的取值范圍位于 [ 0,1]。特別地,當Λ(x,y)=0 時,互信息M(x,y)=0,從其定義可知p(x,y)=p(x)p(y),即x與y相互獨立。

    Λ(x,y)是一個同時包含線性和非線性關系的度量,而互相關系數(shù)r(x,y)和Λ(x,y)之間的差異

    可用于揭示變量之間除線性關系之外更復雜的關聯(lián)關系。

    當已知z∈{z1,z2,...,zn}≡N3時,衡量x與y的互信息即條件互信息[13]有

    轉移熵是一個量化兩個系統(tǒng)相關程度或混沌程度的參數(shù),可用于處理有限長度的信號數(shù)據(jù)。轉移熵定義由Schreiber[15]于2000 年提出,同時考慮源序列和目標序列的數(shù)據(jù)長度,通用的轉移熵定義為(默認單位為nats)

    其 中,yp(t)=[y(t),y(t ??),...,y(t ?k?)],k+1為歷史時間序列的長度,?為最小的取樣間隔,τ為時間延遲。轉移熵可以被視為一種特殊形式的條件互信息[13],即

    轉移熵可以衡量系統(tǒng)中已知目標歷史序列yp的情況下,有多少信息從輸入序列X傳遞到目標序列Y的下一狀態(tài)。轉移熵與相關系數(shù)的不同在于方向性,即TX→Y TY →X。

    轉移熵算法容易受到有限大小樣本效應的影響,也容易受到數(shù)值效應的影響,因此建立一個閾值來界定轉移熵計算結果是否顯著是必要的。為了檢驗計算結果的統(tǒng)計學意義,將源序列x隨機打亂重新抽樣,這一步驟保證了置換時間序列具有與原始序列相同的均值、方差、自相關函數(shù),但破壞了非線性關系。這個過程的目的是消除兩個序列數(shù)據(jù)x與y之間的所有潛在關聯(lián)關系,因此計算得到的Tsur(X)→Y理論值為零。但是在有限時間序列中,由于有限樣本效應Tsur(X)→Y很少為零,故而可以得到轉移熵TX→Y高于Tsur(X)→Y即可認為是顯著的。通過100 次隨機的零假設抽樣計算Tsur(X)→Y,若根據(jù)置換序列計算出的新Tsur(X)→Y值的95%小于原來TX→Y,則可以認為原始序列計算得來的TX→Y顯著。

    1.2.3 算法比較

    運用時移相關系數(shù)、轉移熵分析行星際磁場南向分量Bz與地磁指數(shù)Sym-H之間的相關性。使用1996-2018 年小時精度數(shù)據(jù)對時移相關系數(shù)和轉移熵進行逐年計算。針對每一年的數(shù)據(jù),將時移相關系數(shù)rx→y(τ)=r[x(t),y(t+τ)]的最大時間延遲設置為120 h,以步長1 h 逐步計算,即每年共120 個數(shù)據(jù)點,然后取相關系數(shù)絕對值最大的數(shù)據(jù)點作為衡量該年Bz和地磁指數(shù)Sym-H之間的相關系數(shù)計算結果,有

    同理,計算轉移熵時采用同樣的時間延遲設置和最大值選取方法,由此可得圖1。圖1 中1996 年行星際磁場南向分量與地磁指數(shù)的最大時移相關系數(shù)r(Bz,Sym-H)=0.39;在太陽活動高年2000 年r(Bz,Sym-H)=0.35,同極小年的差異率為11%。另外,1996 年和2000 年行星際磁場南向分量與地磁Sym-H指數(shù)之間的轉移熵T(Bz →Sym-H)分別為0.18 nats、0.26 nats,差異率為31%。相關系數(shù)在太陽活動高低年的差異性相對較小,轉移熵在太陽活動高低年表現(xiàn)出明顯的差異。與傳統(tǒng)相關系數(shù)的對比表明,轉移熵可以有效表征太陽風–磁層關系的相關性。

    如圖1 所示,轉移熵的變化趨勢表現(xiàn)出類似太陽活動周變化的雙峰結構。太陽活動水平在一定程度上影響著空間天氣,太陽黑子數(shù)(Nss)是描述太陽活動水平的重要指標,計算得到太陽黑子和行星際磁場南向分量與地磁Sym-H指數(shù)之間的轉移熵TBz→Sym-H=T(Bz →Sym-H)的相關系數(shù)

    圖1 行星際磁場南向分量Bz 與地磁Sym-H 指數(shù)的轉移熵和相關系數(shù)Fig.1 Transfer entropy and correlation coefficient between I MF Bz and Sym-H

    太陽黑子和行星際磁場南向分量與地磁Sym-H指數(shù)之間的轉移熵TBz→Sym-H的轉移熵

    太陽黑子和行星際磁場南向分量與地磁Sym-H指數(shù)之間的相關系數(shù)rBz,Sym-H=r(Bz,Sym-H)的相關系數(shù):

    對比結果得知,

    以太陽活動周期的角度看,各太陽風參數(shù)向地磁Sym-H指數(shù)信息傳遞的強弱與太陽活動水平的周期性變化一致。另外,太陽10.7 cm 射電流量(F10.7)是綜合衡量太陽色球、過渡區(qū)和日冕極紫外輻射強度的一種常用指數(shù),以F10.7作為衡量太陽活動水平的指標可得

    以F10.7刻畫太陽活動強度,得到與太陽黑子數(shù)相同的結果。因此,轉移熵可以更好地定量化太陽活動與太陽風–磁層系統(tǒng)的相關性。

    2 太陽活動水平對太陽風–磁層相互作用的影響

    以往太陽風–磁層耦合關系多參數(shù)影響因素的重要性研究,往往僅集中在個別太陽風參數(shù),例如行星際磁場南向分量(Bz)或行星際電場(E),而與其他參數(shù),例如行星際磁場大?。˙)、太陽風速度(vsw)、太陽風數(shù)密度(Dsw)、太陽風溫度(Tsw)、太陽風動壓(Psw)等的綜合研究較少涉及。Schwenn 等[16]發(fā)現(xiàn),地磁擾動在大小和方向上都與行星際磁場波動密切相關。長期南向的行星際磁場與地球磁場之間的相互連接允許太陽風能量傳輸?shù)降厍虼艑覽17]。包括Kissinger等[18]在內(nèi)的多項研究表明了太陽風速度在磁暴產(chǎn)生中的作用,事實上持續(xù)和增強的太陽風速度以及南向或北向的行星際磁場分量通常與行星際物質拋射有關,其被認為是磁暴的重要原因。另外,增強的太陽風密度也是一個經(jīng)常影響磁暴強度的重要參數(shù)。太陽風密度的增大會導致日側磁層頂?shù)膲嚎s,從而導致磁層頂流、場向流和橫向流的增加。Ahmed 等[7]的研究表明太陽風密度與其他參數(shù)協(xié)同具有更好的相關性,而太陽風溫度則表現(xiàn)較差。Xie 等[19]多篇論文描述了高速太陽風動壓與磁暴的關系。Mcpherron等[20]的研究表明磁暴時的環(huán)電流積聚原則上可以直接由太陽風電場驅動。Kane[21]表示磁暴期間太陽風速度與行星際磁場南向分量vswBz的共同作用效果明顯優(yōu)于vsw;Wang 等[22]研究表明,對于磁暴的形成,的重要性強于行星際磁場南向分量增強的持續(xù)時間 ?t。

    標識地磁活動水平的典型地磁指數(shù),例如Dst指數(shù),表征了全球對稱赤道環(huán)電流的強度,是四個低緯天文臺的地球磁場水平分量的小時平均值[7],而Sym-H為1 min 分辨率的Dst指數(shù)。

    2.1 太陽風參數(shù)與地磁指數(shù)的轉移熵

    實驗使用1996-2018 年小時精度的太陽風數(shù)據(jù),主要分析參數(shù)包括行星際磁場(B)、行星際磁場南向分量(Bz)、太陽風速度(vsw)、太陽風等離子體數(shù)密度(Dsw)、太陽風溫度(Tsw)、太陽風動壓(Psw)、行星際電場(E),以及小時精度的地磁Sym-H指數(shù)數(shù)據(jù)。

    圖2 展示了利用轉移熵方法分析各太陽風參數(shù)與地磁Sym-H指數(shù)的相關性結果,各太陽風參數(shù)向地磁指數(shù)的轉移熵的相對值表現(xiàn)出一定的穩(wěn)定性,即:1996-2018 年中,行星際電場與地磁指數(shù)之間的信息傳遞TE→Sym-H均最高,隨后是和,最低是和位于中部且相互差異較小,排名在各年中略有不同。各太陽風參數(shù)與地磁Sym-H指數(shù)轉移熵的變化趨勢表現(xiàn)類似于圖1 中行星際磁場南向分量與地磁Sym-H指數(shù)的轉移熵相似的特性:(1)在太陽活動高年,Tsw→Sym-H更為顯著;(2)對于任一太陽風參數(shù),Tsw→Sym-H在第24 太陽活動周的峰值比23 太陽活動周要低。

    圖2 各太陽風參數(shù)對地磁Sym-H 指數(shù)的轉移熵Fig.2 Transfer entropy of multiple solar wind parameters to Sym-H

    2.2 轉移熵與太陽活動水平相關性分析

    如圖2 所示,除行星際磁場南向分量外,行星際磁場、太陽風速度、太陽風等離子體數(shù)密度、行星際電場及地磁Sym-H指數(shù)轉移熵表現(xiàn)出與太陽活動周類似的雙峰結構。同時將Tsw→Sym-H與太陽活動指數(shù)太陽黑子數(shù)、F10.7進行相關性分析,結果列于表1。

    由 表1 可 知,r(Nss,Tsw→Sym-H)位 于[0.61,0.83]之間、r(F10.7,Tsw→Sym-H)位于[0.58,0.82]之間,均表現(xiàn)出較強的相關性。同時獲得T(Nss,Tsw→Sym-H)位于[0.15,0.36]之間、T(F10.7,Tsw→Sym-H)位于[0.12,0.33]之間,即太陽活動指數(shù)與Tsw→Sym-H存在明顯的信息傳遞。由此可知,以小時精度的數(shù)據(jù)分析,行星際磁場南向分量與太陽活動之間的強相關并非偶然,其他各太陽風參數(shù)與太陽活動水平也存在較強的關聯(lián),佐證了Liu 等[23]對于地球近地空間環(huán)境變化對太陽活動的依賴性。

    表1 T sw→Sym-H同太陽活動指數(shù)的相關性Table 1 Correlation betweenTsw→Sym-H and solar radiation

    3 引起磁暴擾動的太陽風參數(shù)重要性排序

    3.1 數(shù)據(jù)選擇

    磁暴是指整個地球磁層發(fā)生的持續(xù)十幾個小時到幾十個小時的一種劇烈地磁擾動,中低緯度地磁臺站水平分量的顯著減小為磁暴的主要特征。Wu 等[24]的研究表明,太陽擾動日變化在主相期間最強,在恢復相期間逐漸減弱。在磁暴初相和主相期間,一般會出現(xiàn)太陽風–磁層發(fā)電機效應,將導致磁層對流電場快速滲透到中低緯度地區(qū);高緯焦耳加熱產(chǎn)生全球熱層擾動場,主要出現(xiàn)在恢復相期間[25]。磁暴主相期間,低緯水平磁場擾動(?H)在各磁地方時(MLT)均明顯下降并且出現(xiàn)晨昏不對稱分布;磁暴恢復相期間,行星際電場(E)為負值時,?H較主相期間弱。已有的研究表明:主相階段主要源于環(huán)電流的增強,而影響恢復相的主要為電荷損失機制,包括電荷交換、庫侖散射和波粒相互作用等[26]。因此,由于不同的物理成因,應分別研究磁暴主相和恢復相,目前本文將研究范圍限定于磁暴主相階段。

    一般根據(jù)Dst峰值的不同,將地磁暴分為小磁暴(?50 nT

    3.2 重要性排序方法

    以往研究結果表明,地磁指數(shù)對太陽風變化的響應時間普遍小于3 h[11,28],本研究將響應時間范圍設定在0 至500 min,與Stumpo 等[11]研究磁層電離層對太陽風的響應時間范圍一致。以往的轉移熵往往集中在單個時間序列,而磁暴期間的重要性系數(shù)為磁暴集合(見表2)上的太陽風參數(shù)對地磁Sym-H指數(shù)的轉移熵平均值,故而定義

    表2 2010-2018 年磁暴事件Table 2 Geomagnetic storm events from 2010 to 2018

    其中〈·〉代表在表2 磁暴集合上的平均值。

    相應地由轉移熵定義的重要性系數(shù)為

    以及定義信噪比(S/N),顯著性檢驗指標

    轉移熵算法通過置換時間序列[15]來創(chuàng)建零假設下的概率分布來計算噪聲,噪聲的均值和標準偏差是由100 次隨機的零假設抽樣計算Tsur(X)→Sym-H平均而來。

    3.3 重要性排序

    對行星際磁場南向分量(Bz)和地磁Sym-H指數(shù)之間的轉移熵進行計算,分析兩者之間的信息流動。首先對每個磁暴事件進行獨立的轉移熵分析,而后按照式(12)以相同的時間延遲τ計算多個磁暴事件轉移熵(τ)的平均值(τ)。由圖3 可以看出,轉移熵在τ=60 min時達到最大值Tpeak= 0.244 nats、Γ=0.196 nats,信噪比δS/N=4.15和顯著性αsig=36.83σ。轉移熵明顯大于噪聲,即信息傳遞顯著。而基本都處于背景噪聲的內(nèi)部或邊緣,不顯示任何明顯的時間延遲特征,信息傳遞不顯著。

    圖3 行星際磁層南向分量 Bz和 地磁Sym-H指數(shù)的轉移熵Fig.3 Transfer entropy of B z and Sym-H

    經(jīng)比較,上述結果與Runge 等[28]的結論基本一致,即行星際磁場南向分量(Bz)擾動是導致行星際條件變化的重要部分;Stumpo 等[11]的研究也表明,在時間延遲τ=60 min 時行星際磁場南向分量(Bz)與Sym-H指數(shù)之間有較強的信息傳遞。但是已有研究并未涉及太陽風溫度(Tsw)、行星際電場(E)等太陽風參數(shù)的綜合比較。進一步利用轉移熵分析磁暴期間其他太陽風參數(shù),包括行星際磁場(B)、太陽風速度(vsw)、太陽風等離子體數(shù)密度(Dsw)、太陽風溫度(Tsw)、太陽風動壓(Psw)和行星際電場(E),對Sym-H指數(shù)的影響,結果如圖4 所示。

    由圖4 可知,各太陽風參數(shù)對地磁Sym-H指數(shù)的轉移熵均高于噪聲;而則處于噪聲邊緣,缺乏顯著性。最大值出現(xiàn)的時間普遍小于3 h,反應出磁層對太陽風變異性的響應是突變式的激增,表現(xiàn)為短時間內(nèi)一個強烈的信息傳遞。表3 根據(jù)各太陽風參數(shù)到地磁Sym-H指數(shù)的轉移熵,得出了各太陽風參數(shù)的重要性系數(shù)序列。

    由表3 可知,以轉移熵構建的重要性系數(shù)為指標,排在第一位的是行星際電場(E),ΓE=0.200 nats。第二位是行星際磁場南向分量(Bz),ΓBz=0.196 nats,略低于ΓE。然后排名依次是太陽風速度(vsw)、太陽風溫度(Tsw)、太陽風密度(Dsw)、行星際磁場(B)。最低的是太陽風動壓(Psw),ΓPsw=0.135 nats。

    值得注意的是圖4(f),行星際電場(E)對地磁Sym-H指數(shù)的轉移熵在τ=60 min時達到最大值Tpeak= 0.252 nats、Γ=0.200 nats,信噪比δS/N= 3.92 和顯著性αsig=37.43σ,也有明顯的信息傳遞,甚至比更強烈。

    表3 中排在第三位的是太陽風速度(vsw),結合圖4(d)可知在τ=140 min 時達到最大值Tpeak= 0.210 nats、Γ=0.159 nats,信噪比δS/N=3.15和顯著性αsig=29.92σ。行星際磁場南向分量和磁層的重聯(lián)是引發(fā)磁暴等空間天氣事件的重要因素。Zhang 等[4]進一步得出,行星際磁場南向分量與太陽風速度的聯(lián)合作用又體現(xiàn)于行星際電場的時變曲線,統(tǒng)一表現(xiàn)為行星際電場,太陽風速度的變化影響著磁層對流電場,從而影響環(huán)電流的形成快慢與強度。因此,從信息傳遞的角度分析,行星際電場擾動是磁暴發(fā)生的主要驅動因素。以往太陽風溫度(Tsw)、太陽風密度(Dsw)、行星際磁場(B)、太陽風動壓(Psw)等綜合研究較少,從轉移熵的角度出發(fā),本文給出了這些太陽風參數(shù)在引發(fā)磁暴時的相對重要性。

    圖4 行星際磁場(B )、太陽風速度(vsw)、太陽風等離子體數(shù)密度(D sw )、太陽風溫度(Tsw)、太陽風動壓(Psw)、行星際電場(E)對地磁Sym-H 指數(shù)的轉移熵Fig.4 Transfer entropy of interplanetary magnetic field (B),solar wind speed (vsw),solar wind plasma number density (D sw) ,solar wind temperature (Tsw) solar wind dynamic pressure (Psw),the interplanetary electric field (E) and Sym-H

    表3 太陽風參數(shù)對地磁Sym-H 指數(shù)重要性系數(shù)排序Table 3 Importance ranking of solar wind parameters to geomagnetic index Sym-H

    以相關系數(shù)為基礎進行排序,由表3 中的|r|列可知,前三位仍是行星際電場(E)、行星際磁場南向分量(Bz)和太陽風速度(vsw),但是行星際磁場(B)的相對重要性提升至第四位,太陽風溫度(Tsw)下降至第六位。根據(jù)1.2 節(jié)的討論可知,相關系數(shù)僅能挖掘變量之間的線性關系,而從概率角度出發(fā)的轉移熵可以挖掘變量之間更一般的非線性依賴關系,轉移熵和相關系數(shù)排序的差異,表明了太陽風參數(shù)與地磁Sym-H指數(shù)之間關系的復雜性。進一步通過式(7)計算DM來對差異進行定量刻畫,DM >0 揭示太陽風–磁層之間線性關系之外更復雜的關聯(lián)性。由此可知,相對于相關系數(shù),轉移熵方法可以更多地挖掘參數(shù)間的非線性關系,故而導致了重要性排序的差異。

    綜上所述,行星際電場(E)的變化是引發(fā)地磁暴的主要驅動因素,這與Zhao 等[5]的研究結果基本一致。行星際磁場南向分量(Bz)與地磁Sym-H指數(shù)之間的重要性系數(shù)略低于ΓE,但是明顯高于其他太陽風參數(shù)的Γsw。太陽風動壓(Psw)與地磁Sym-H指數(shù)之間的重要性系數(shù)最低;其他各太陽風參數(shù)在引發(fā)磁暴時的相對關聯(lián)程度由表3 給出。

    4 結論

    應用轉移熵研究了地磁擾動的太陽風驅動因素,可以得出如下結論。

    (1)從太陽活動周期的角度分析,各太陽風參數(shù)向地磁Sym-H指數(shù)信息傳遞的強弱與太陽活動水平的周期性變化一致。

    (2)磁暴期間各太陽風參數(shù)向地磁Sym-H指數(shù)信息傳遞,最相關的變量是行星際電場(E),行星際磁場南向分量(Bz)次之,然后排名依次是太陽風速度(vsw)、太陽風溫度(Tsw)、太陽風等離子體數(shù)密度(Dsw)、行星際磁場(B),最弱的是太陽風動壓(Psw),各參數(shù)信息傳遞強弱排名見表3。

    (3)行星際電場與地磁指數(shù)之間的信息傳遞TE→Sym-H在τmax=60 min 處達到峰值,比TBz→Sym-H略高,佐證了Zhao 等[5]的研究。

    在研究太陽風–磁層系統(tǒng)響應太陽風條件變化時,運用轉移熵可以推斷太陽風–磁層的非線性關系。本文將轉移熵方法推廣至磁暴事件集合上的平均轉移熵,定義了磁暴期間太陽風參數(shù)與地磁Sym-H指數(shù)的重要性系數(shù),并綜合考慮了7 種太陽風參數(shù),以重要性系數(shù)為依據(jù)獲得了引發(fā)磁暴擾動的太陽風參數(shù)重要性排序。轉移熵定量描述太陽風與磁層參數(shù)之間的非線性信息傳遞,為探索磁層對行星際太陽風擾動的響應過程提供一條數(shù)據(jù)驅動的研究途徑,為模型研究提供參數(shù)選擇的依據(jù),由此可加速建模過程。未來轉移熵在空間天氣領域可以通過更多事件統(tǒng)計研究開展廣泛應用。

    致謝數(shù)據(jù)分析環(huán)境由國家科技資源共享服務平臺之一的國家空間科學數(shù)據(jù)中心(http://www.nssdc.ac.cn)提供。

    猜你喜歡
    磁層磁暴太陽活動
    小諸葛漫“話”天氣之極光與地磁暴
    科學24小時(2024年5期)2024-05-17 12:59:56
    航天器的“隱形殺手”:地磁暴
    軍事文摘(2023年14期)2023-07-28 08:39:58
    《從太陽風暴到行星磁暴》專刊(Ⅲ)寄語
    第24太陽活動周中國地區(qū)電離層閃爍統(tǒng)計特性研究
    基于對抗神經(jīng)網(wǎng)絡有限角度的磁層邊界CT重構技術
    第23和24太陽活動周高緯地磁感應電流分布特性
    磁鞘快速流引起的磁層頂凹陷事件
    等離子體鞘層效應對磁層探測電場儀設計的影響
    地球磁層亞暴統(tǒng)計分析
    磁暴恢復相后期超低頻波觀測分析
    一区二区三区高清视频在线| 久久人人爽人人爽人人片va | 免费看日本二区| 国产黄色小视频在线观看| 一区福利在线观看| 亚洲最大成人中文| 欧美黄色淫秽网站| 国产精品自产拍在线观看55亚洲| 国产不卡一卡二| 三级国产精品欧美在线观看| 久久久久久久亚洲中文字幕 | 中文字幕高清在线视频| x7x7x7水蜜桃| 两个人视频免费观看高清| 一区福利在线观看| 在线观看66精品国产| 日韩欧美国产在线观看| 能在线免费观看的黄片| 欧美日韩乱码在线| 久久婷婷人人爽人人干人人爱| 人妻夜夜爽99麻豆av| 亚洲欧美日韩高清在线视频| 久久性视频一级片| 88av欧美| 国产一级毛片七仙女欲春2| av专区在线播放| 如何舔出高潮| 麻豆国产av国片精品| 一级毛片久久久久久久久女| 亚洲精品一区av在线观看| 国产探花在线观看一区二区| 非洲黑人性xxxx精品又粗又长| 18禁在线播放成人免费| 午夜a级毛片| 91午夜精品亚洲一区二区三区 | 一本精品99久久精品77| 免费在线观看亚洲国产| 啪啪无遮挡十八禁网站| 99热只有精品国产| 成人国产一区最新在线观看| 日韩精品中文字幕看吧| 99热这里只有是精品50| 久99久视频精品免费| 日本a在线网址| 婷婷六月久久综合丁香| 波野结衣二区三区在线| 欧美中文日本在线观看视频| 十八禁网站免费在线| 免费无遮挡裸体视频| 尤物成人国产欧美一区二区三区| 99国产精品一区二区蜜桃av| 国产精品一及| 欧美日韩黄片免| 午夜日韩欧美国产| 可以在线观看的亚洲视频| 黄片小视频在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产色片| 99久久精品热视频| 免费看光身美女| 国产成人啪精品午夜网站| 久久亚洲真实| 一区福利在线观看| 国产伦在线观看视频一区| 国产日本99.免费观看| 五月伊人婷婷丁香| 久久久久久久久久成人| 日韩亚洲欧美综合| 欧美精品国产亚洲| 在线观看免费视频日本深夜| 午夜精品一区二区三区免费看| 成年人黄色毛片网站| 淫秽高清视频在线观看| 少妇人妻一区二区三区视频| 国产精品亚洲一级av第二区| 少妇的逼水好多| 99热这里只有是精品在线观看 | www.熟女人妻精品国产| 亚洲人与动物交配视频| 亚洲av第一区精品v没综合| 国产高清视频在线观看网站| 久久亚洲精品不卡| av天堂在线播放| 一进一出抽搐gif免费好疼| av天堂中文字幕网| 9191精品国产免费久久| 国内精品久久久久精免费| 精品日产1卡2卡| 日韩有码中文字幕| 毛片女人毛片| 桃红色精品国产亚洲av| 一区二区三区四区激情视频 | 美女免费视频网站| 婷婷精品国产亚洲av在线| 女同久久另类99精品国产91| netflix在线观看网站| 亚洲在线自拍视频| 免费观看的影片在线观看| 国产伦精品一区二区三区视频9| 又紧又爽又黄一区二区| 天堂av国产一区二区熟女人妻| 亚洲男人的天堂狠狠| 免费在线观看成人毛片| 午夜免费激情av| 亚洲欧美日韩卡通动漫| 亚洲最大成人av| 九色成人免费人妻av| 日韩欧美免费精品| 精品一区二区免费观看| 好男人电影高清在线观看| 日韩中文字幕欧美一区二区| 伦理电影大哥的女人| 国产精品人妻久久久久久| 99精品在免费线老司机午夜| 国产精品永久免费网站| 色av中文字幕| 日本在线视频免费播放| 免费在线观看成人毛片| 热99re8久久精品国产| 久久午夜亚洲精品久久| 国产精华一区二区三区| 毛片一级片免费看久久久久 | 国产真实伦视频高清在线观看 | 日韩免费av在线播放| 在线观看免费视频日本深夜| 国产精品98久久久久久宅男小说| 亚洲av五月六月丁香网| 欧美精品啪啪一区二区三区| 国产探花极品一区二区| 国产精品人妻久久久久久| 欧美高清成人免费视频www| 免费av观看视频| 国产伦在线观看视频一区| 午夜免费男女啪啪视频观看 | 全区人妻精品视频| 两性午夜刺激爽爽歪歪视频在线观看| 真人做人爱边吃奶动态| 黄色丝袜av网址大全| 老司机午夜十八禁免费视频| 亚洲精品久久国产高清桃花| 我要看日韩黄色一级片| 国产一区二区亚洲精品在线观看| 日韩欧美三级三区| 国产精品一及| 国产成年人精品一区二区| 国产不卡一卡二| 午夜福利高清视频| 国产午夜福利久久久久久| 伦理电影大哥的女人| 国内精品一区二区在线观看| 国产成+人综合+亚洲专区| 亚洲成av人片在线播放无| 国产伦精品一区二区三区四那| 夜夜躁狠狠躁天天躁| 日本一本二区三区精品| 国产激情偷乱视频一区二区| 欧美日韩综合久久久久久 | 久久中文看片网| 精品熟女少妇八av免费久了| 白带黄色成豆腐渣| 色综合站精品国产| 色av中文字幕| 国产三级中文精品| 97热精品久久久久久| 国产精品野战在线观看| 美女免费视频网站| 一进一出抽搐动态| 日本熟妇午夜| 88av欧美| 99久久精品国产亚洲精品| 看十八女毛片水多多多| 中文字幕熟女人妻在线| 又爽又黄a免费视频| 久久久久久九九精品二区国产| 免费看美女性在线毛片视频| 欧美一级a爱片免费观看看| 2021天堂中文幕一二区在线观| 欧美性猛交╳xxx乱大交人| 色尼玛亚洲综合影院| 国产免费一级a男人的天堂| 听说在线观看完整版免费高清| 日韩欧美在线二视频| 人妻夜夜爽99麻豆av| 天美传媒精品一区二区| 日日夜夜操网爽| 18禁黄网站禁片免费观看直播| 69人妻影院| 欧美丝袜亚洲另类 | 亚洲真实伦在线观看| 亚洲av日韩精品久久久久久密| 高清毛片免费观看视频网站| 国产精品久久电影中文字幕| 看免费av毛片| 国产一区二区三区视频了| 88av欧美| 久久天躁狠狠躁夜夜2o2o| 黄色丝袜av网址大全| ponron亚洲| 亚洲七黄色美女视频| 国产色婷婷99| 国产伦精品一区二区三区视频9| 男人和女人高潮做爰伦理| 久久久久久久久久成人| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| 亚洲 国产 在线| 嫩草影院入口| 看黄色毛片网站| 热99在线观看视频| 桃色一区二区三区在线观看| 日韩欧美 国产精品| 51午夜福利影视在线观看| 日本一本二区三区精品| 无遮挡黄片免费观看| 一本一本综合久久| 一级黄色大片毛片| av视频在线观看入口| 窝窝影院91人妻| 精品国产三级普通话版| 国产亚洲欧美98| 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| 首页视频小说图片口味搜索| 黄色日韩在线| a级毛片a级免费在线| 国产在视频线在精品| 99国产极品粉嫩在线观看| 中国美女看黄片| 9191精品国产免费久久| 国产成人福利小说| 国内久久婷婷六月综合欲色啪| 亚洲精品在线观看二区| 精品熟女少妇八av免费久了| 亚洲电影在线观看av| 天美传媒精品一区二区| 成人一区二区视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 别揉我奶头~嗯~啊~动态视频| 日韩精品中文字幕看吧| 亚洲av一区综合| 怎么达到女性高潮| 国产伦精品一区二区三区视频9| 免费在线观看成人毛片| 午夜福利高清视频| 亚洲精品久久国产高清桃花| 伦理电影大哥的女人| 亚洲av电影在线进入| 亚洲人成网站高清观看| 亚洲av.av天堂| 国产单亲对白刺激| 日本五十路高清| 在线观看午夜福利视频| 欧美绝顶高潮抽搐喷水| 欧美性猛交╳xxx乱大交人| www.熟女人妻精品国产| 中文字幕人成人乱码亚洲影| 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸| 51午夜福利影视在线观看| 精品久久久久久久人妻蜜臀av| 亚洲综合色惰| 别揉我奶头 嗯啊视频| 亚洲真实伦在线观看| 真人一进一出gif抽搐免费| 国产精品一区二区三区四区免费观看 | 国产在线精品亚洲第一网站| 免费av毛片视频| 亚洲一区高清亚洲精品| 亚洲精品亚洲一区二区| 亚洲国产色片| 99视频精品全部免费 在线| а√天堂www在线а√下载| 成人精品一区二区免费| 国产国拍精品亚洲av在线观看| 国产av一区在线观看免费| 色在线成人网| 色视频www国产| 91久久精品电影网| 1024手机看黄色片| 国产精品久久久久久久电影| 久久午夜福利片| 久久精品国产亚洲av涩爱 | 精品福利观看| 日本a在线网址| 亚洲av中文字字幕乱码综合| 精品乱码久久久久久99久播| 中文在线观看免费www的网站| 人人妻,人人澡人人爽秒播| 欧美色视频一区免费| 欧美不卡视频在线免费观看| 露出奶头的视频| 99久久成人亚洲精品观看| 亚洲国产精品sss在线观看| 亚洲人与动物交配视频| 久久久久久九九精品二区国产| 一本久久中文字幕| 婷婷亚洲欧美| 亚洲国产精品成人综合色| 国产午夜精品论理片| 久久精品综合一区二区三区| 日日摸夜夜添夜夜添小说| 美女高潮喷水抽搐中文字幕| 直男gayav资源| 日韩中字成人| 男人狂女人下面高潮的视频| 人人妻,人人澡人人爽秒播| 中亚洲国语对白在线视频| 欧美日韩瑟瑟在线播放| av视频在线观看入口| 搡老熟女国产l中国老女人| 国产一区二区亚洲精品在线观看| 国产综合懂色| 色吧在线观看| 男女那种视频在线观看| 国产91精品成人一区二区三区| 人妻久久中文字幕网| 国产男靠女视频免费网站| 一进一出抽搐动态| 欧美+亚洲+日韩+国产| 日韩亚洲欧美综合| 999久久久精品免费观看国产| 亚洲精品久久国产高清桃花| 久久香蕉精品热| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| a级毛片a级免费在线| 国产欧美日韩精品一区二区| 琪琪午夜伦伦电影理论片6080| 亚洲三级黄色毛片| 国产成人福利小说| 午夜视频国产福利| 欧美性猛交黑人性爽| 免费观看精品视频网站| 麻豆国产97在线/欧美| 亚洲在线观看片| 欧美黄色片欧美黄色片| 十八禁国产超污无遮挡网站| 波多野结衣高清作品| 亚洲天堂国产精品一区在线| 免费av毛片视频| 直男gayav资源| 精品国产三级普通话版| 色尼玛亚洲综合影院| 久久精品夜夜夜夜夜久久蜜豆| 特大巨黑吊av在线直播| 一区福利在线观看| 丰满人妻熟妇乱又伦精品不卡| 搡女人真爽免费视频火全软件 | 老鸭窝网址在线观看| 国产毛片a区久久久久| 日日干狠狠操夜夜爽| 天天一区二区日本电影三级| 一夜夜www| 国产亚洲精品综合一区在线观看| 舔av片在线| 午夜视频国产福利| 乱人视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲aⅴ乱码一区二区在线播放| 757午夜福利合集在线观看| 69人妻影院| 国产精品乱码一区二三区的特点| 丝袜美腿在线中文| 三级男女做爰猛烈吃奶摸视频| 啪啪无遮挡十八禁网站| 国产老妇女一区| 最后的刺客免费高清国语| 欧美区成人在线视频| 成人无遮挡网站| 免费av不卡在线播放| 99久久无色码亚洲精品果冻| 国产精品亚洲一级av第二区| 男人和女人高潮做爰伦理| 身体一侧抽搐| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av在线| 免费在线观看成人毛片| 亚洲精品日韩av片在线观看| 亚洲av电影不卡..在线观看| 国产伦精品一区二区三区视频9| 久久精品国产亚洲av天美| 国产伦精品一区二区三区视频9| 国产精品永久免费网站| 深夜精品福利| 欧美潮喷喷水| 亚洲性夜色夜夜综合| 亚洲成av人片在线播放无| 欧美性猛交╳xxx乱大交人| 欧美日韩国产亚洲二区| 亚洲欧美日韩卡通动漫| 成人国产综合亚洲| 欧美xxxx黑人xx丫x性爽| 国产精品亚洲一级av第二区| 久9热在线精品视频| 999久久久精品免费观看国产| 国产精华一区二区三区| 日韩免费av在线播放| 精品熟女少妇八av免费久了| 国产三级中文精品| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| 中出人妻视频一区二区| 日本 av在线| 日本三级黄在线观看| 亚洲av.av天堂| 淫妇啪啪啪对白视频| 精品久久久久久久久av| 高潮久久久久久久久久久不卡| 长腿黑丝高跟| 国产精品不卡视频一区二区 | 少妇高潮的动态图| 久久久久久大精品| 亚洲va日本ⅴa欧美va伊人久久| 国产精品女同一区二区软件 | 中文字幕人成人乱码亚洲影| 校园春色视频在线观看| 国产91精品成人一区二区三区| 天堂网av新在线| 国产亚洲欧美在线一区二区| 搞女人的毛片| 深爱激情五月婷婷| 激情在线观看视频在线高清| 天天躁日日操中文字幕| 国内精品美女久久久久久| 一进一出抽搐动态| 国产一级毛片七仙女欲春2| 一区二区三区免费毛片| 国产精品99久久久久久久久| 我要看日韩黄色一级片| 欧美不卡视频在线免费观看| 青草久久国产| 91字幕亚洲| 两人在一起打扑克的视频| 18禁黄网站禁片免费观看直播| 久久精品国产亚洲av天美| 日韩欧美三级三区| 久久久久久久久中文| 女生性感内裤真人,穿戴方法视频| 国产精品人妻久久久久久| 嫩草影院新地址| 婷婷亚洲欧美| 一本久久中文字幕| 欧美+亚洲+日韩+国产| 久久99热6这里只有精品| 亚洲精品日韩av片在线观看| 国产精品久久电影中文字幕| 香蕉av资源在线| 亚洲美女搞黄在线观看 | 久久婷婷人人爽人人干人人爱| 观看免费一级毛片| 国产一区二区亚洲精品在线观看| 99国产精品一区二区三区| 色精品久久人妻99蜜桃| 亚洲中文字幕日韩| 在线免费观看不下载黄p国产 | 免费人成视频x8x8入口观看| 全区人妻精品视频| 国产精品美女特级片免费视频播放器| 欧美又色又爽又黄视频| 成人亚洲精品av一区二区| 国产精品,欧美在线| 淫妇啪啪啪对白视频| 久久久久久大精品| 免费无遮挡裸体视频| 亚洲最大成人手机在线| 国产三级中文精品| 亚洲久久久久久中文字幕| 美女高潮喷水抽搐中文字幕| 桃色一区二区三区在线观看| 中文字幕av成人在线电影| 国产69精品久久久久777片| 欧美又色又爽又黄视频| 亚洲不卡免费看| 国产欧美日韩精品亚洲av| 国产熟女xx| 性插视频无遮挡在线免费观看| 精品一区二区三区人妻视频| 国产精品亚洲av一区麻豆| www.色视频.com| 欧美最黄视频在线播放免费| 天美传媒精品一区二区| 亚洲中文字幕日韩| 99热这里只有是精品在线观看 | 国产av麻豆久久久久久久| 我要搜黄色片| 成人特级黄色片久久久久久久| 美女 人体艺术 gogo| 美女免费视频网站| 天堂网av新在线| 欧美一区二区亚洲| 国产亚洲精品久久久久久毛片| 三级毛片av免费| 一区福利在线观看| 国产av麻豆久久久久久久| 一夜夜www| 日韩欧美精品v在线| 亚洲不卡免费看| 国产成人福利小说| 国产精品国产高清国产av| 久久午夜亚洲精品久久| 亚洲欧美日韩无卡精品| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 九九在线视频观看精品| xxxwww97欧美| 国产精品久久久久久精品电影| 午夜视频国产福利| 老鸭窝网址在线观看| 国产爱豆传媒在线观看| 99热这里只有精品一区| 国产精品三级大全| 免费无遮挡裸体视频| 简卡轻食公司| 午夜a级毛片| 欧美潮喷喷水| 亚洲黑人精品在线| 亚洲专区中文字幕在线| 亚洲美女搞黄在线观看 | 欧美黑人欧美精品刺激| 国产精品99久久久久久久久| 亚洲最大成人手机在线| 美女大奶头视频| 国产 一区 欧美 日韩| 最近视频中文字幕2019在线8| 欧美最新免费一区二区三区 | 精品一区二区三区人妻视频| 99久国产av精品| 我的女老师完整版在线观看| 在线播放无遮挡| 国产成人影院久久av| 人人妻人人澡欧美一区二区| 天堂影院成人在线观看| 欧美一区二区亚洲| 国产极品精品免费视频能看的| 性欧美人与动物交配| 国产精品免费一区二区三区在线| 国产午夜精品论理片| 禁无遮挡网站| 18禁黄网站禁片免费观看直播| 国内精品久久久久精免费| 69人妻影院| 国产高潮美女av| 亚洲男人的天堂狠狠| 亚洲精品粉嫩美女一区| 亚洲va日本ⅴa欧美va伊人久久| 日本 欧美在线| 亚洲成人免费电影在线观看| 最近在线观看免费完整版| 精品无人区乱码1区二区| 一进一出抽搐gif免费好疼| 亚洲人成伊人成综合网2020| 男女下面进入的视频免费午夜| 九色国产91popny在线| 日韩欧美精品v在线| 内地一区二区视频在线| ponron亚洲| 亚洲人与动物交配视频| 色在线成人网| 亚洲第一电影网av| 日本精品一区二区三区蜜桃| 免费在线观看日本一区| 国产高清激情床上av| 亚洲人成网站高清观看| 露出奶头的视频| 日韩国内少妇激情av| 成人无遮挡网站| 亚洲欧美日韩东京热| 一区二区三区免费毛片| 国产高清视频在线播放一区| 成人性生交大片免费视频hd| 亚洲精品亚洲一区二区| 夜夜夜夜夜久久久久| 午夜a级毛片| 日韩大尺度精品在线看网址| 一个人免费在线观看的高清视频| 国产探花极品一区二区| 丝袜美腿在线中文| 国内精品一区二区在线观看| www.www免费av| 亚洲一区二区三区色噜噜| 能在线免费观看的黄片| 国产野战对白在线观看| 国产v大片淫在线免费观看| www.色视频.com| 免费观看精品视频网站| 久久久精品欧美日韩精品| 757午夜福利合集在线观看| 夜夜爽天天搞| 国产午夜福利久久久久久| 国产麻豆成人av免费视频| 成人三级黄色视频| 黄色丝袜av网址大全| 国产欧美日韩精品亚洲av| 国产亚洲精品av在线| 男女视频在线观看网站免费| 永久网站在线| 在线免费观看不下载黄p国产 | 日韩中文字幕欧美一区二区| 国产美女午夜福利| 亚洲成人中文字幕在线播放| 日韩欧美国产在线观看| 国产一区二区激情短视频| 欧美绝顶高潮抽搐喷水| 午夜老司机福利剧场| 1000部很黄的大片| 三级毛片av免费| 欧美午夜高清在线| 成人永久免费在线观看视频| 简卡轻食公司| 免费电影在线观看免费观看| 成人高潮视频无遮挡免费网站| 国内精品久久久久久久电影| 欧美精品国产亚洲| 在线观看66精品国产| 如何舔出高潮|