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

    EAST 上中性束注入和離子回旋共振加熱下快離子分布函數(shù)層析反演*

    2023-11-24 05:05:46孫延旭黃娟高偉常加峰張偉史唱李云鶴
    物理學(xué)報 2023年21期
    關(guān)鍵詞:層析正則反演

    孫延旭 黃娟 高偉 常加峰 張偉 史唱 李云鶴

    1) (中國科學(xué)院合肥物質(zhì)科學(xué)研究院,等離子體物理研究所,合肥 230031)

    2) (中國科學(xué)技術(shù)大學(xué),合肥 230026)

    3) (北京大學(xué),核物理與技術(shù)國家重點實驗室,北京 100871)

    聚變等離子體中快離子分布函數(shù)的速度空間層析反演(tomography)是研究磁約束核聚變裝置中快離子分布和輸運的重要手段.在東方超環(huán)(experimental advanced superconducting tokamak,EAST)中性束注入(neutral beam injection,NBI)與離子回旋共振加熱(ion cyclotron range of frequencies heating,ICRF)協(xié)同加熱實驗中,快離子診斷測量信號以及TRANSP 模擬獲得的快離子分布函數(shù)中觀測到協(xié)同效應(yīng)產(chǎn)生的高于中性束注入能量的高能粒子.為了研究快離子分布行為,獲得不同加熱條件下診斷測量的快離子分布函數(shù),采用不同的正則化方法,增加先驗信息以及將快離子Dα 光譜診斷(fast-ion Dα spectroscopy,FIDA)與中子發(fā)射譜儀(neutron emission spectrometers,NES)相結(jié)合等方式,有效提高快離子診斷信噪比以及在速度空間權(quán)重覆蓋率,實現(xiàn)在單獨NBI 加熱以及NBI 和ICRF 協(xié)同加熱條件下基于診斷測量的層析反演,獲得真實可靠的快離子分布函數(shù).這為下一步提高NBI 與ICRF 協(xié)同加熱效率,研究協(xié)同加熱機制以及相關(guān)的快離子分布和輸運行為奠定基礎(chǔ).

    1 引言

    在磁約束核聚變裝置中,由中性束注入(NBI)和離子回旋共振加熱(ICRF)等輔助加熱產(chǎn)生的快離子在提高電流驅(qū)動和等離子體加熱效率中起著重要作用[1].協(xié)同效應(yīng)是通過ICRF 的快波電場梯度加速位于高次諧波共振層附近的束離子,而協(xié)同效應(yīng)產(chǎn)生高能離子尾巴,提高等離子體性能的機制在DIII-D[2],ASDEX Upgrade[3],JET[4,5]等各大裝置上得到實驗驗證.各種高能粒子診斷,如快離子Dα(FIDA)診斷[6,7]、相干湯姆遜散射(CTS)[8,9]、伽馬射線光譜(GRS)[10,11]和中子發(fā)射光譜(NES)[12]等被應(yīng)用于協(xié)同加熱下的快離子分布信息的測量.

    高能粒子診斷獲得的是快離子分布函數(shù)局部的一維投影或其他衍生量,無法直接得到分布函數(shù).通過速度空間層析反演(tomography)方法[13,14]對快離子診斷測量進行反演獲得速度空間二維快離子分布函數(shù)已經(jīng)成功應(yīng)用于各大托卡馬克裝置,包括ASDEX Upgrade[15?17],DIII-D[18]和MAST[19].特別是在ASDEX Upgrade,擁有五視場FIDA 診斷,為反演提供豐富的數(shù)據(jù)集,在速度空間有足夠的覆蓋率,提高了反演的精度.以前的研究表明,在速度空間探測范圍不完全覆蓋的情況下,從診斷系統(tǒng)中重建速度分布函數(shù)是十分困難的.而通過在快離子速度空間層析反演方法中結(jié)合不同診斷方法的測量,可以提高速度空間的分辨率和覆蓋率,在ASDEX Upgrade 上結(jié)合FIDA 和CTS 診斷[20]、JET 上已經(jīng)成功地結(jié)合NES 和GRS 診斷實現(xiàn)速度空間分布函數(shù)的反演[21].在EAST 托卡馬克上,盡管雙視場FIDA 診斷的速度空間覆蓋率有限,但通過在零階Tikhonov 正則化層析反演方法中添加先驗信息以及將FIDA 與NES 診斷結(jié)合,實現(xiàn)了NBI 加熱下快離子分布函數(shù)重建[22,23].同時,隨著EAST 上NBI 以及ICRF 加熱效率的提高[24],最近在EAST 進行的實驗中觀察到通過NBI-ICRF協(xié)同加熱,NBI 的能量低于80 keV 的束離子被加速至300 keV 以上,同時提高了儲能、中子產(chǎn)額等等離子體參數(shù),進一步擴展了EAST 裝置運行區(qū)間[25],相應(yīng)的協(xié)同加熱下的快離子分布的測量與反演亟待解決.

    為了提高EAST 上基于快離子診斷測量的速度空間層析反演方法的可靠性,獲得NBI-ICRF協(xié)同加熱下的快離子分布函數(shù).本文基于快離子診斷實驗測量數(shù)據(jù)采用一階Tikhonov 正則化方法,通過加入額外的先驗信息,例如空值測量(null measurement)、速度空間NBI 峰值位置(known peaks)等,以及將FIDA 與NES 診斷相結(jié)合的方法提高EAST 上反演快離子分布函數(shù)的精度,獲得NBI 單獨加熱以及NBI-ICRF 協(xié)同加熱下的快離子分布函數(shù),驗證EAST 上NBI-ICRF 協(xié)同加熱效應(yīng),為開展協(xié)同加熱下快離子行為物理研究,探索提高協(xié)同效應(yīng)以及加熱效率的方法奠定基礎(chǔ).

    2 快離子速度空間層析反演

    2.1 FIDA 和NES 診斷

    FIDA 診斷系統(tǒng)通過測量快離子與注入的中性束發(fā)生電荷交換反應(yīng)后發(fā)出的Balmer-alpha 輻射,從測量光譜中的多普勒頻移獲得關(guān)于快離子速度分布的信息.而安裝在H 窗口的NES 診斷可以測量從視線中逃出的聚變產(chǎn)生的中子[26].圖1(a)俯視圖和圖1(b)極向截面圖描述了EAST 上FIDA 和NES 診斷系統(tǒng)的總體布局,以及本文中使用放電.在EAST 上包括安裝在O 窗口的切向視圖和B 窗口的垂直視圖兩套FIDA 診斷系統(tǒng)[27],FIDA 以A 窗口NBI1L 為診斷束,診斷視場與托卡馬克內(nèi)局部磁場的夾角范圍是10°—20° (O 窗口),75°—84° (B 窗口),中平面上的測量光斑尺寸為3 cm.分別覆蓋徑向位置177—237 cm 和177—205 cm 的區(qū)域,兩個窗口觀測位置重疊以實現(xiàn)快離子分布函數(shù)的反演.由于FIDA 測量的信噪比通常在靠近等離子體芯部區(qū)域較好,同時那里的快離子密度通常是最高的,所以本文采用O 窗口和B 窗口FIDA 視圖的芯部測量,其大半徑R分別約為195 cm 和194 cm(歸一化小半徑ρ~ 0.2),對應(yīng)圖1 中紅色視線觀測位置.EAST 上FIDA 診斷使用NBI1L 束調(diào)制的時間片減法(束開-束關(guān))去除復(fù)雜的背景輻射,同時在光譜儀中使用2 nm 寬的陷波濾波器過濾掉來自邊界的波長為656.1 nm的冷Dα光.

    圖1 EAST 托卡馬克上四條中性束束線位置、O 和B 窗口FIDA 視圖以及NES 的布局 (a)俯視圖;(b) 極向截面圖.視場中的紅線是本文所分析的信號的觀測位置Fig.1.Layout of four neutral beams,the O and B-port view of FIDA,the location of the NES on EAST: (a) Top view;(b) poloidal view.The red lines-of-sight are analyzed in this paper for shot #113648.

    為了從診斷測量光譜中獲得快離子分布函數(shù),需要引入一個權(quán)重函數(shù)W,權(quán)重函數(shù)代表診斷測量的速度空間靈敏度,在特定速度空間的快離子與注入的中性束發(fā)生電荷交換反應(yīng)被中和,并在特定波長范圍內(nèi)發(fā)射 Dα光子的概率.通過構(gòu)建的權(quán)重函數(shù)W將快離子速度空間分布函數(shù)與特定光譜范圍內(nèi)的信號聯(lián)系起來:

    式中S(λ) 為測量信號;f(E,p) 是任意的快離子速度分布函數(shù);E代表快離子能量;p是俯仰角(v∥/v),其中v為快離子速度,v∥是與磁場平行的速度分量,λ 代表波長.圖2(a)、圖2(b)和圖2(c)分別為NES 診斷以及FIDA 切向視圖、垂直視圖在速度空間權(quán)重覆蓋率,由該視圖的所有測量相關(guān)的權(quán)重函數(shù)相加得到.以前的研究表明,速度空間權(quán)重函數(shù)的充分覆蓋對于利用層析反演提供可靠的重建至關(guān)重要[9].

    圖2 歸一化的權(quán)重函數(shù)總和 (a) NES H 窗口;(b) FIDA O 窗口;(c) FIDA B 窗口Fig.2.Normalized weight function coverage: (a) NES H-port;(b) FIDA O-port;(c) FIDA B-port.

    2.2 速度空間快離子分布函數(shù)層析反演

    速度空間層析反演需要解決一個不適定問題,依賴于確定尋求方程F*,S=WF*的數(shù)學(xué)上不穩(wěn)定的解.在這個方程中S代表測量信號和快離子分布函數(shù)F*之間的關(guān)系,傳遞矩陣W為已知的權(quán)重函數(shù).為了獲得這個不適定問題的穩(wěn)定解決方案,需要進行正則化.以前的研究表明,使用Tikhonov 正則化方法[28]可以在托卡馬克等離子體中實現(xiàn)快離子分布函數(shù)的可靠反演.這種方法的關(guān)鍵在于解決最小化問題:

    這里,F是任意的快離子速度分布函數(shù),根據(jù)懲罰矩陣L的不同分為零階(L為單位矩陣)和一階(L是對梯度進行有限差分近似的矩陣算子)正則化[28],而λ1表示自由正則化參數(shù).正則化參數(shù)平衡解對含有噪聲的測量數(shù)據(jù)的擬合程度和對正則化約束的服從程度.但是由于測量光譜經(jīng)常受到軔致輻射和雜質(zhì)線輻射的背景干擾,以及FIDA 診斷視場、實驗測量數(shù)據(jù)集不足的限制,基于標(biāo)準(zhǔn)的Tikhonov正則化反演方法一直受到偽影的困擾,出現(xiàn)負(fù)的快離子相空間密度等情況.因此通過添加各種已知的先驗信息,彌補測量數(shù)據(jù)的缺乏,改善偽影問題[13].為了增強反演的效果,尤其是在EAST上雙視圖FIDA,權(quán)重函數(shù)在速度空間覆蓋率不足的情況下,在使用正則化方法時添加額外的先驗信息.聚變等離子體中快離子的速度空間層析反演已經(jīng)考慮了非負(fù)值、空測量和已知峰值位置等先驗信息[19].首先,為了解決使用Tikhonov 正則化方法帶來的負(fù)相空間密度,需要施加非負(fù)性的約束條件,即不允許負(fù)相空間密度的區(qū)域的出現(xiàn),即F*≥ 0.此外實驗上FIDA 診斷測量到的多普勒頻移是有上限的.在光譜中觀察到的低于噪聲底線的FIDA光的部分被稱為空值測量,其波長范圍可以通過權(quán)重函數(shù)與速度空間區(qū)域相關(guān),與空值測量區(qū)域?qū)?yīng)的速度空間內(nèi)不可能出現(xiàn)較高的快離子密度.而NBI注入粒子的全能量、半能量和三分之一能量成分比值以及能量大小已知,以之可獲得中性束已知峰值位置的先驗信息.

    空值測量曾被作為硬約束納入反演.確定速度空間的空值測量區(qū)域,忽略空測量權(quán)重函數(shù)所覆蓋區(qū)域的快離子相空間密度,使用F(E0,p0)=0,其中E0和p0是空測權(quán)重函數(shù)覆蓋的能量和間距[19].雖然這種約束可以避免空值測量區(qū)域偽影的出現(xiàn),并減少反演中的未知數(shù),但不能排除空值測量區(qū)域的快離子的存在,因此采用另一種方法,通過引入一個旨在抑制速度空間空值測量區(qū)域中規(guī)范解的懲罰函數(shù),允許在速度空間的空測區(qū)域存在快離子[19]:

    單位矩陣用L0表示,而懲罰強度由正則化參數(shù)λ0決定,ξ=ξ(E,p)代表空值測量區(qū)域的懲罰函數(shù).懲罰函數(shù)ξ 被選擇為隨著空值測量區(qū)域能量的增加而呈指數(shù)增長,有效地抑制了空值測量區(qū)域的速度空間中的明顯偽影.這種方法也實現(xiàn)了密集和稀疏速度空間區(qū)域之間的平滑過渡,消除了空值測量速度空間內(nèi)的尖銳邊界[13,19].一階Tikhonov正則化,對速度空間的陡峭梯度進行懲罰,可能會減少NBI 的峰值,因此反演在NBI 條件下產(chǎn)生的快離子的分布函數(shù)的峰值,特別是在NBI 能量的全量、半量和三分之一處是十分困難的.由于NBI 峰值位置的E,p是已知的.通過采用二維函數(shù)κ(E,p)來編碼NBI 的已知幾何形狀,減少了峰值位置的懲罰強度,允許分布函數(shù)中更大梯度的存在[13].結(jié)合非負(fù)性和空值測量的約束,分布函數(shù)的解表示為

    3 NBI 以及NBI-ICRF 協(xié)同加熱下的快離子分布函數(shù)

    在EAST 中性束注入系統(tǒng)以及ICRF 天線升級后,針對NBI 單獨加熱以及NBI-ICRF 協(xié)同加熱開展了實驗研究,基于快離子診斷測量數(shù)據(jù)開展協(xié)同加熱以及NBI 單獨加熱下的快離子分布函數(shù)反演工作,研究表明ICRF 對注入的束離子進行加速,產(chǎn)生了高于NBI 能量的快離子,EAST 上協(xié)同加熱效果得到實驗驗證.圖3 為EAST 上#113648放電的實驗參數(shù)設(shè)置以及中子產(chǎn)額和極向比壓βp的測量結(jié)果.此炮放電的等離子體電流Ip是400 kA,環(huán)向磁場Bt是2.43 T,電流方向和磁場方向標(biāo)注于圖1(a),分別是逆時針方向和順時針方向.

    圖3 #113648 放電參數(shù)時序演化圖 (a) 中子產(chǎn)額;(b) 極向比壓;(c) 中性束注入功率;(d) ICRF 加熱功率Fig.3.Time traces of shot #113648: (a) Neutron yield;(b) poloidal beta;(c) neutral beam injection power;(d) ICRF heating power.

    診斷束NBI1L 以調(diào)制注入的方式實現(xiàn)FIDA診斷的有效測量,NBI2L 和NBI2R 以40 kV 的束壓注入進行加熱,ICRF 功率為1 MW,通過對加熱時序的調(diào)整,在保持其他等離子體參數(shù)不變的情況下實現(xiàn)一炮內(nèi)NBI 單獨加熱(圖3(c)中黑色箭頭,2.72 s)和NBI-ICRF 協(xié)同加熱(圖3(c)中藍色箭頭,3.72 s).實驗發(fā)現(xiàn)在NBI-ICRF 協(xié)同加熱時等離子體中子產(chǎn)額和βp分別提高約57%和32%,等離子體性能得到提高.為了驗證ICRF 加速束離子的協(xié)同加熱機制,研究不同加熱下的快離子分布函數(shù)的區(qū)別,接下來將對NBI 單獨以及協(xié)同加熱下的診斷測量數(shù)據(jù)進行反演.

    NBI-ICRF 協(xié)同加熱時,利用FIDA 以及NES測量協(xié)同效應(yīng)產(chǎn)生的高能粒子的分布信息,如圖4 雙視圖FIDA 測量結(jié)果所示,相較于NBI 單獨加熱,協(xié)同加熱時FIDA 光譜強度的提升反映快離子密度的增加,圖中誤差棒為診斷束打開區(qū)間內(nèi)測量光譜的整體標(biāo)準(zhǔn)差.同時NBI 能量為40 keV的快離子在忽略斯塔克分裂的情況下最大多普勒頻移為4.1 nm,對應(yīng)于Dα線左側(cè)的652 nm 和右側(cè)的660.2 nm.在NBI 單獨注入時,僅在這些波長范圍內(nèi)測量到FIDA 光.在NBI-ICRH 加熱階段,可以觀測到這些邊界(圖4 中藍色虛線)以上的FIDA 光譜尾巴,尤其在對捕獲離子敏感的垂直觀測方向,這表明存在高于注入能量的快離子.同時,如圖5 所示,在NBI-ICRF 加熱階段,NES 測量到快離子與背景熱離子反應(yīng)產(chǎn)生的中子,在強度和能量上限上都遠超于NBI 單獨加熱階段,EAST上協(xié)同加熱效果得到驗證.由于#113648 放電過程中沒有磁流體不穩(wěn)定性(MHD)行為的出現(xiàn),不用考慮快離子的反常輸運,因此,可以認(rèn)為此次放電滿足新經(jīng)典輸運模型.在完成動理學(xué)平衡處理后,將湯姆孫散射診斷測量的電子溫度、彎晶譜儀測量的離子溫度剖面以及偏振干涉儀、微波反射計系統(tǒng)測量的電子密度剖面等基本的剖面分布以及電流、環(huán)電壓等放電基本信息輸入至TRANSP,利用NUBEAM[29],TORIC 模塊對不同加熱階段的快離子分布函數(shù)進行模擬,獲得的速度空間快離子分布函數(shù)如圖6 所示.

    圖5 單獨NBI 以及NBI-ICRF 協(xié)同加熱時NES 能譜Fig.5.NES spectrum of NBI only and NBI-ICRF synergetic heating.

    圖6 TRANSP 模擬的快離子分布函數(shù) (a)單獨NBI 加熱;(b) NBI-ICRF 協(xié)同加 熱;黑色虛線為NBI 能 量,黑 色實線為單獨NBI 加熱時快離子能量上限Fig.6.TRANSP simulation of fast ion distribution function: (a) NBI only heating;(b) NBI-ICRF synergetic heating.Black dashed line is the NBI injection energy,black solid line is the upper limit of fast ion energy for NBI only.

    由于EAST 上的中性束皆為同電流方向注入,單獨NBI 加熱時,快離子主要分布在負(fù)的俯仰角區(qū)域(以Bt方向為正),圖中的黑色實線為NBI 單獨加熱階段快離子分布函數(shù)的能量上限,從圖6(b)模擬的分布函數(shù)以及圖7 沿負(fù)俯仰角積分的快離子分布能譜圖也可以看出,NBI-ICRF 協(xié)同加熱產(chǎn)生了高于NBI 能量的快離子,提高了高能量快離子的份額.

    圖7 沿負(fù)俯仰角積分的快離子分布能譜圖.Fig.7.Fast ion distribution energy spectrum integrated along the negative pitch angle.

    利用2.2 節(jié)的速度空間層析反演方法,對NBI單獨加熱階段的雙視角的FIDA 測量數(shù)據(jù),在不添加先驗條件的情況下,分別利用零階和一階Tikhonov 正則化方法進行反演,由于層析反演方法獲得的分布函數(shù)的能量、俯仰角分辨率即描述分布函數(shù)的速度空間(E,p)網(wǎng)格大小,取決于用于反演的實驗數(shù)據(jù)測量點的數(shù)量[14],如果網(wǎng)格點多于測量值,權(quán)重函數(shù)的特征也會出現(xiàn)在層析成像反演中[30],影響反演結(jié)果的真實性.而EAST 上FIDA診斷僅有兩個視場,用于反演的實驗數(shù)據(jù)點為446,因此設(shè)置速度空間(E,p)反演分布函數(shù)的網(wǎng)格大小為20×20 (小于測量數(shù)據(jù)點446),此時能量分辨率為5 keV.如圖8(a)和圖8(b)所示,無論是零階或一階Tikhonov 正則化方法,其反演結(jié)果都存在嚴(yán)重的偽影問題,即在高能量區(qū)域存在較高的遠大于NBI 能量的快離子密度.但是一階Tikhonov正則化反演結(jié)果在較低能量區(qū)域(E小于50 keV)反演精度要高于零階正則化方法.在一階Tikhonov正則化方法的基礎(chǔ)上,通過添加非負(fù)性條件以及基于FIDA 診斷測量光譜獲得的空值測量區(qū)域的先驗信息,反演帶來的偽影問題得到了解決,如圖8(c),反演與TRANSP 模擬的新經(jīng)典的快離子分布函數(shù)基本一致(圖6(a)).但是一階正則化導(dǎo)致在速度空間中NBI 峰值位置的陡峭梯度被削弱,峰值信息未能反演出來.通過引入已知的NBI 的中性粒子的全能量、半能量、三分之一能量信息作為先驗條件,能夠有效地重建出NBI 峰值信息,如圖8(d)所示.采用這種方法,實現(xiàn)了在EAST 上NBI單獨注入條件下實現(xiàn)快離子分布函數(shù)的可靠重建,提高了反演的精度,但是由于分布函數(shù)的能量、俯仰角分辨率由實驗測量數(shù)據(jù)點的數(shù)量決定,因此能量分辨率并未提高.

    圖8 單獨NBI 時基于FIDA 診斷測量的快離子分布函數(shù)層析反演 (a)零階Tikhonov 正則化;(b) 一階Tikhonov 正則化;(c) 一階Tikhonov 正則化,空值測量區(qū)域;(d) 一階Tikhonov 正則化,已知峰值位置Fig.8.Reconstructions of the fast ion distribution function based on FIDA diagnostic measurements during NBI only heating:(a) Zero-order Tikhonov regularization;(b) first-order Tikhonov regularization;(c) first-order Tikhonov regularization,null measurement;(d) first-order Tikhonov regularization,known peak.

    FIDA 診斷數(shù)據(jù)空值測量區(qū)域這個先驗信息適用于快離子速度分布函數(shù)中存在一個較為急劇的躍變到零的情況.這個條件在NBI 單獨加熱時得到了滿足.而在NBI-ICRF 協(xié)同加熱的實驗中,由于產(chǎn)生了高于NBI 能量的快離子,在FIDA 測量的能量范圍內(nèi)快離子分布函數(shù)不會有躍變,此時空值測量區(qū)域不再適用于基于FIDA 診斷測量的層析反演.由于#113648 放電,加熱束的束壓較低僅為40 keV,同時基于TRANSP 模擬和FIDA 測量發(fā)現(xiàn)了在百keV 量級內(nèi)快離子分布的能量上限,因此,本文在對NBI-ICRF 協(xié)同加熱階段快離子分布函數(shù)反演中,使用TRANSP 模擬的快離子分布函數(shù)能量上限作于空值測量的約束條件,驗證EAST 上速度空間層析反演算法在協(xié)同加熱時的可行性.

    首先采用合成信號進行反演,即根據(jù)(1)式將TRANSP 模擬的快離子分布函數(shù)與權(quán)重函數(shù)結(jié)合,獲得合成信號,并添加10%的噪聲進行反演.在NBI-ICRF 協(xié)同加熱階段,在不添加先驗信息的情況下,如圖9(b)所示,基于一階Tikhonov 正則化方法的層析反演效果優(yōu)于零階(圖9(a)),反演結(jié)果與模擬結(jié)果基本一致,反演的誤差在20%以內(nèi),但是在俯仰角為?0.5的區(qū)域,仍然存在較大的偽影,這個問題在添加基于TRANSP 模擬的空值測量區(qū)域后得到改善(圖9(c)).進一步基于EAST上雙視圖FIDA 測量數(shù)據(jù)展開一階Tikhonov 正則化層析反演,不添加先驗信息的反演結(jié)果如圖10(a)所示,圖中虛線代表TRANSP 模擬的快離子能量上限.由于EAST 上僅有兩個視場的FIDA,如圖2(a)和圖2(b)所示,在速度空間俯仰角為0 附近,僅有垂直B 窗口一個視場覆蓋,這給反演帶來了挑戰(zhàn).圖10(a)中,在權(quán)重函數(shù)覆蓋不足的俯仰角0 至1 范圍內(nèi),反演存在嚴(yán)重的偽影問題.在權(quán)重覆蓋率相對較大的俯仰角為?0.7 附近,偽影相對較小且反演與TRANSP 模擬吻合.為了增加速度空間的權(quán)重覆蓋率,采用將FIDA 與NES 診斷測量相結(jié)合進行反演,結(jié)果如圖10(b)所示,在速度空間內(nèi)反演的偽影問題得到了明顯的改善,與模擬的快離子分布函數(shù)誤差進一步減小,但是由于NES 在百keV 內(nèi)的探測效率較低(從圖2(a) NES的總的權(quán)重函數(shù)可以看出),因此在這種快離子能量上限較低的情況下,并不能從NES 測量中獲得空值測量的先驗信息,偽影并不能完全消除,因此,使用TRANSP 模擬的能量上限作為空值測量的先驗條件,極大地提升了反演的精度,如圖10(c)和圖10(d)所示,同時結(jié)合FIDA 與NES 診斷提升了反演在較高能量范圍(50 keV)的反演效果.因此,基于EAST 上現(xiàn)有的FIDA 與NES 診斷,利用一階Tikhonov 正則化,FIDA 與NES 相結(jié)合的層析反演方法,實現(xiàn)在NBI-ICRF 協(xié)同加熱條件下,獲得真實可靠的快離子分布函數(shù).

    圖9 NBI-ICRF 協(xié)同加熱階段基于TRANSP 模擬合成信號的層析反演 (a) 零階Tikhonov 正則化;(b) 一階Tikhonov 正則化;(c) 一階Tikhonov 正則化,空值測量區(qū)域Fig.9.Reconstructions of synthetic signal based on the TRANSP distributions: (a) Zero-order Tikhonov regularization;(b) first-order Tikhonov regularization;(c) first-order Tikhonov regularization,null measurement.

    圖10 基于快離子診斷數(shù)據(jù)的分布函數(shù)層析反演 (a),(b) 不添加先驗信息;(c),(d) TRANSP 模擬的快離子能量上限作為先驗信息;其中(a),(c)僅使用FIDA 測量;(b),(d) FIDA 與NES 測量結(jié)合Fig.10.Reconstructions based on fast ion diagnostic data: (a),(b) Without adding a priori information;(c),(d) fast ion energy upper bounds for TRANSP simulations as prior information.Only FIDA measurement is used in panel (a) and (c);FIDA in combination with NES measurements are used in panel (b) and (d).

    4 結(jié)論

    聚變等離子體中通過速度空間層析反演方法獲得快離子分布函數(shù)通常受到快離子Dα(FIDA)診斷視場和信噪比的限制.本文在單獨中性束加熱的情況下,通過對比不同正則化方法,采用反演精度更高的一階Tikhonov 正則化層析反演方法,進一步通過添加空值測量區(qū)域,已知的NBI 峰值位置等先驗信息,彌補EAST 上雙視圖FIDA 在速度空間測量數(shù)據(jù)的不足,解決偽影問題,在速度空間預(yù)期的位置上重建了NBI 峰,反演結(jié)果得到了明顯的改善.同時,隨著EAST 上中性束注入與離子回旋共振加熱系統(tǒng)的性能,加熱效率的提升,快離子診斷FIDA 和NES 等測量發(fā)現(xiàn)了明顯的NBIICRF 協(xié)同效應(yīng),產(chǎn)生了遠高于中性束注入能量的快離子.為了驗證NBI-ICRF 協(xié)同效應(yīng),研究協(xié)同加熱下快離子分布機制,本文在一階Tikhonov 正則化層析反演方法的基礎(chǔ)上,通過將FIDA 與NES診斷結(jié)合,偽影問題得到改善,獲得了真實可靠的基于實驗測量的快離子速度空間分布函數(shù),并將TRANSP 模擬得到的快離子能量上限作為先驗信息,彌補NES 診斷在百keV 內(nèi)探測效率的不足,提升了反演獲得的快離子分布函數(shù)精度,對協(xié)同加熱效果進行了驗證.證明層析反演方法在EAST協(xié)同加熱高參數(shù)運行條件下獲得快離子分布函數(shù)的可行性,為下一步開展氫氘比、中性束注入能量對協(xié)同加熱效率的影響,協(xié)同加熱下快離子分布和輸運機制的實驗研究奠定基礎(chǔ).

    感謝丹麥科技大學(xué)的Mirko Salewski 教授對速度空間層析反演方法研究的支持,感謝EAST團隊對實驗運行和數(shù)據(jù)獲取的支持.

    猜你喜歡
    層析正則反演
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    犬細(xì)小病毒量子點免疫層析試紙條的研制
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    A族鏈球菌膠體金免疫層析試紙條的制備及應(yīng)用
    新型B族鏈球菌膠體金免疫層析試紙條的臨床應(yīng)用評價
    有限秩的可解群的正則自同構(gòu)
    疊前同步反演在港中油田的應(yīng)用
    亚洲丝袜综合中文字幕| 亚洲国产成人一精品久久久| 精品久久久精品久久久| 国产免费福利视频在线观看| 国产成人精品福利久久| 亚洲国产日韩欧美精品在线观看| 久久久久久久亚洲中文字幕| 中文字幕久久专区| 免费av观看视频| 欧美激情久久久久久爽电影| 国产精品1区2区在线观看.| 永久网站在线| 日韩av在线免费看完整版不卡| 亚洲熟女精品中文字幕| 精品久久久久久电影网| 777米奇影视久久| 乱人视频在线观看| 熟女人妻精品中文字幕| 国产免费又黄又爽又色| 免费观看av网站的网址| 久久人人爽人人爽人人片va| 亚洲人与动物交配视频| 91精品一卡2卡3卡4卡| 久久久久久九九精品二区国产| av又黄又爽大尺度在线免费看| 午夜福利在线观看吧| 高清毛片免费看| 日韩,欧美,国产一区二区三区| 婷婷色麻豆天堂久久| 国产一级毛片在线| 欧美日韩亚洲高清精品| 亚洲色图av天堂| 日本色播在线视频| 久久99热这里只频精品6学生| 三级国产精品片| 中文字幕av在线有码专区| 看免费成人av毛片| 精品久久久久久久久亚洲| 久久综合国产亚洲精品| 18+在线观看网站| 国产高清不卡午夜福利| 夫妻性生交免费视频一级片| 中文字幕亚洲精品专区| 十八禁国产超污无遮挡网站| 又黄又爽又刺激的免费视频.| a级一级毛片免费在线观看| 亚洲图色成人| 亚洲图色成人| 日韩电影二区| 十八禁国产超污无遮挡网站| 亚洲丝袜综合中文字幕| 国产精品99久久久久久久久| 夫妻性生交免费视频一级片| 午夜日本视频在线| 最新中文字幕久久久久| 777米奇影视久久| 嫩草影院新地址| 亚洲精品日本国产第一区| 欧美3d第一页| 在线观看人妻少妇| 亚洲,欧美,日韩| 国产午夜福利久久久久久| 亚洲精品色激情综合| 亚洲av免费高清在线观看| 成人美女网站在线观看视频| 亚洲av中文av极速乱| 男女国产视频网站| 国产精品国产三级国产专区5o| 91精品一卡2卡3卡4卡| 国产亚洲5aaaaa淫片| 乱人视频在线观看| 日本三级黄在线观看| 最近中文字幕2019免费版| 一级毛片 在线播放| 亚洲va在线va天堂va国产| 丝袜喷水一区| 日韩av不卡免费在线播放| 国产视频首页在线观看| 亚洲人成网站在线播| 夫妻性生交免费视频一级片| 久久人人爽人人爽人人片va| 啦啦啦啦在线视频资源| 少妇高潮的动态图| 国产精品人妻久久久影院| 精品欧美国产一区二区三| 黄色欧美视频在线观看| 只有这里有精品99| 毛片一级片免费看久久久久| 国产日韩欧美在线精品| av播播在线观看一区| 国产黄色视频一区二区在线观看| 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 日韩av免费高清视频| 水蜜桃什么品种好| 内射极品少妇av片p| 日本午夜av视频| 免费观看无遮挡的男女| 亚洲av成人精品一区久久| 国产高潮美女av| 丰满乱子伦码专区| 欧美丝袜亚洲另类| 亚洲欧美成人精品一区二区| 日韩精品有码人妻一区| 国产人妻一区二区三区在| 日韩在线高清观看一区二区三区| 亚洲欧洲国产日韩| 成年人午夜在线观看视频 | 国产极品天堂在线| 毛片一级片免费看久久久久| 日韩一区二区视频免费看| 亚洲精品乱码久久久久久按摩| 亚洲av成人av| 国产国拍精品亚洲av在线观看| 精品不卡国产一区二区三区| 肉色欧美久久久久久久蜜桃 | 国产不卡一卡二| 精品久久久久久久末码| 精品国内亚洲2022精品成人| av国产免费在线观看| 成人午夜高清在线视频| 日韩电影二区| 日本与韩国留学比较| 日本爱情动作片www.在线观看| 成年av动漫网址| 高清午夜精品一区二区三区| 乱人视频在线观看| 观看免费一级毛片| av播播在线观看一区| 亚洲国产精品成人综合色| 老女人水多毛片| 日韩欧美 国产精品| 国产高清国产精品国产三级 | 乱系列少妇在线播放| 国产亚洲精品av在线| av免费观看日本| 最后的刺客免费高清国语| 精品国内亚洲2022精品成人| 男女下面进入的视频免费午夜| 亚洲va在线va天堂va国产| av专区在线播放| 国产一区二区三区av在线| 六月丁香七月| kizo精华| 一级毛片 在线播放| 水蜜桃什么品种好| 国产乱人偷精品视频| 国产精品1区2区在线观看.| 亚洲18禁久久av| 三级国产精品片| 一个人看的www免费观看视频| 亚洲18禁久久av| 久久久久性生活片| 高清午夜精品一区二区三区| 欧美性感艳星| 国产成人a区在线观看| 亚洲精品日韩在线中文字幕| 熟女人妻精品中文字幕| 天天一区二区日本电影三级| 又粗又硬又长又爽又黄的视频| 国产乱来视频区| 国产精品一区二区三区四区免费观看| 日韩三级伦理在线观看| 国产精品久久久久久久久免| 夫妻性生交免费视频一级片| 国产精品三级大全| 直男gayav资源| 我的老师免费观看完整版| 舔av片在线| 国产一级毛片在线| 寂寞人妻少妇视频99o| 亚洲av电影在线观看一区二区三区 | 老司机影院成人| 久久亚洲国产成人精品v| 哪个播放器可以免费观看大片| 97热精品久久久久久| 久久久精品免费免费高清| 美女大奶头视频| 午夜福利在线观看吧| 天堂影院成人在线观看| 国产在线男女| 最近2019中文字幕mv第一页| 青春草国产在线视频| 国产大屁股一区二区在线视频| 国产精品一区二区在线观看99 | 99久久精品热视频| 日韩av在线大香蕉| 如何舔出高潮| 联通29元200g的流量卡| 久久99精品国语久久久| 亚洲国产精品成人综合色| 欧美xxxx性猛交bbbb| 国产一区有黄有色的免费视频 | 只有这里有精品99| 91精品伊人久久大香线蕉| 人妻少妇偷人精品九色| 国产色婷婷99| 伊人久久精品亚洲午夜| 日韩大片免费观看网站| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 国产白丝娇喘喷水9色精品| 亚洲精品久久久久久婷婷小说| 超碰av人人做人人爽久久| 人妻系列 视频| 久久久久免费精品人妻一区二区| 日韩欧美 国产精品| 最近视频中文字幕2019在线8| 18禁动态无遮挡网站| 日韩国内少妇激情av| 蜜臀久久99精品久久宅男| av在线天堂中文字幕| 国产亚洲av嫩草精品影院| 在线播放无遮挡| freevideosex欧美| 丰满乱子伦码专区| 国产中年淑女户外野战色| 久久久精品免费免费高清| 日本免费在线观看一区| or卡值多少钱| 久久99热6这里只有精品| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 嫩草影院新地址| 成人二区视频| 欧美激情国产日韩精品一区| 街头女战士在线观看网站| 欧美成人a在线观看| av线在线观看网站| 夜夜看夜夜爽夜夜摸| 日本wwww免费看| 日本一二三区视频观看| 久久久欧美国产精品| .国产精品久久| 久久久久久久久大av| 国语对白做爰xxxⅹ性视频网站| 亚洲丝袜综合中文字幕| 国产白丝娇喘喷水9色精品| 日产精品乱码卡一卡2卡三| 欧美一区二区亚洲| 国产成人aa在线观看| 国产在视频线在精品| 欧美xxⅹ黑人| 久久这里有精品视频免费| 国产黄片美女视频| 99久久精品热视频| av国产久精品久网站免费入址| 特大巨黑吊av在线直播| 99久久九九国产精品国产免费| 日韩中字成人| 亚洲高清免费不卡视频| videos熟女内射| 高清午夜精品一区二区三区| 日韩av免费高清视频| 日韩一区二区三区影片| 久久久久免费精品人妻一区二区| 亚洲国产精品sss在线观看| 91久久精品国产一区二区成人| 国产在线一区二区三区精| 国产有黄有色有爽视频| 亚洲18禁久久av| 国产中年淑女户外野战色| 亚洲在久久综合| 有码 亚洲区| 伊人久久国产一区二区| 在线观看一区二区三区| 嘟嘟电影网在线观看| 久久久亚洲精品成人影院| 性色avwww在线观看| 人妻制服诱惑在线中文字幕| 免费大片18禁| 尾随美女入室| or卡值多少钱| 人人妻人人看人人澡| 亚州av有码| 亚洲人成网站在线播| 2021天堂中文幕一二区在线观| 国产淫片久久久久久久久| 国产精品一区二区在线观看99 | 久久99热6这里只有精品| 两个人的视频大全免费| 男人和女人高潮做爰伦理| 成年版毛片免费区| 啦啦啦中文免费视频观看日本| 波多野结衣巨乳人妻| 亚洲18禁久久av| 国产在线一区二区三区精| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类| 日韩精品青青久久久久久| 大香蕉97超碰在线| 亚洲第一区二区三区不卡| 国产成人a∨麻豆精品| 亚洲性久久影院| 久久精品综合一区二区三区| a级一级毛片免费在线观看| 麻豆久久精品国产亚洲av| 国产黄片视频在线免费观看| 中国美白少妇内射xxxbb| 日本猛色少妇xxxxx猛交久久| av黄色大香蕉| 日日撸夜夜添| 中文字幕人妻熟人妻熟丝袜美| 最新中文字幕久久久久| 中文字幕免费在线视频6| 久久99精品国语久久久| 五月天丁香电影| 久久精品国产亚洲av涩爱| 爱豆传媒免费全集在线观看| 一级黄片播放器| 国产黄片美女视频| 国产高清国产精品国产三级 | 国产一区二区亚洲精品在线观看| 婷婷色av中文字幕| 夫妻性生交免费视频一级片| 亚洲国产欧美在线一区| 国产精品国产三级国产专区5o| 少妇高潮的动态图| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲精品av在线| 亚洲国产欧美人成| 赤兔流量卡办理| 九色成人免费人妻av| 国产日韩欧美在线精品| 亚洲av成人av| 老女人水多毛片| 国产精品国产三级专区第一集| 久久久久久久久久人人人人人人| 2021天堂中文幕一二区在线观| 在线观看免费高清a一片| av线在线观看网站| 肉色欧美久久久久久久蜜桃 | av免费观看日本| 69av精品久久久久久| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看| 少妇人妻一区二区三区视频| av一本久久久久| 乱系列少妇在线播放| 欧美日本视频| 热99在线观看视频| 国产精品久久久久久精品电影小说 | 成人毛片60女人毛片免费| 在线观看美女被高潮喷水网站| 日韩欧美国产在线观看| freevideosex欧美| 嫩草影院入口| 国产伦精品一区二区三区四那| 日本色播在线视频| 大香蕉97超碰在线| 亚洲性久久影院| 爱豆传媒免费全集在线观看| a级毛色黄片| 三级男女做爰猛烈吃奶摸视频| 色综合站精品国产| 国产v大片淫在线免费观看| 亚洲自拍偷在线| 免费看av在线观看网站| 精华霜和精华液先用哪个| 久久草成人影院| 亚洲精品日韩av片在线观看| 亚洲人成网站在线播| 成人国产麻豆网| 九九爱精品视频在线观看| 免费播放大片免费观看视频在线观看| 亚洲精品一区蜜桃| 亚州av有码| 国产成人精品一,二区| 国产精品人妻久久久影院| 国产日韩欧美在线精品| 国产成人福利小说| 97超碰精品成人国产| 狂野欧美激情性xxxx在线观看| 18+在线观看网站| 中国国产av一级| 人妻一区二区av| 青青草视频在线视频观看| 亚洲性久久影院| 偷拍熟女少妇极品色| 伦精品一区二区三区| 久久99热这里只频精品6学生| 亚洲在久久综合| 免费观看性生交大片5| 少妇熟女aⅴ在线视频| 搡老妇女老女人老熟妇| 精品酒店卫生间| 免费观看无遮挡的男女| 国产老妇女一区| 午夜免费激情av| 永久网站在线| 三级毛片av免费| 国产精品国产三级专区第一集| 韩国高清视频一区二区三区| 五月伊人婷婷丁香| 尤物成人国产欧美一区二区三区| 国产成年人精品一区二区| 日韩电影二区| 国产成人福利小说| 高清在线视频一区二区三区| 日本免费在线观看一区| 精品不卡国产一区二区三区| 淫秽高清视频在线观看| av线在线观看网站| 国产三级在线视频| 亚洲精品国产av蜜桃| eeuss影院久久| 精品99又大又爽又粗少妇毛片| 婷婷六月久久综合丁香| 国产高清三级在线| 国产黄a三级三级三级人| 人人妻人人看人人澡| 欧美一区二区亚洲| 看十八女毛片水多多多| 男插女下体视频免费在线播放| 精品一区二区三区人妻视频| 在线免费观看不下载黄p国产| 国产黄频视频在线观看| 熟妇人妻不卡中文字幕| 亚洲精品色激情综合| 可以在线观看毛片的网站| 狂野欧美激情性xxxx在线观看| 国产精品99久久久久久久久| 亚洲婷婷狠狠爱综合网| 精品一区二区三区人妻视频| 午夜激情福利司机影院| 国产真实伦视频高清在线观看| 18禁动态无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 老师上课跳d突然被开到最大视频| av又黄又爽大尺度在线免费看| 美女被艹到高潮喷水动态| 最近视频中文字幕2019在线8| 亚洲,欧美,日韩| 欧美zozozo另类| 色网站视频免费| 99久久精品国产国产毛片| 成人毛片a级毛片在线播放| 午夜激情欧美在线| 五月伊人婷婷丁香| 尤物成人国产欧美一区二区三区| 国产成人福利小说| 国产探花极品一区二区| 国语对白做爰xxxⅹ性视频网站| 国产av码专区亚洲av| 成人av在线播放网站| 日本与韩国留学比较| 成年av动漫网址| 免费不卡的大黄色大毛片视频在线观看 | 国产高清不卡午夜福利| 午夜精品在线福利| 国产黄色免费在线视频| 狠狠精品人妻久久久久久综合| 婷婷六月久久综合丁香| 免费人成在线观看视频色| av在线天堂中文字幕| 亚洲色图av天堂| 最近手机中文字幕大全| 神马国产精品三级电影在线观看| 国产精品蜜桃在线观看| 亚洲av免费在线观看| 久久久欧美国产精品| 啦啦啦中文免费视频观看日本| 国产老妇女一区| 久99久视频精品免费| 亚洲熟女精品中文字幕| 国产综合懂色| 国产高清国产精品国产三级 | 亚洲欧美成人精品一区二区| 午夜福利视频精品| 纵有疾风起免费观看全集完整版 | 美女国产视频在线观看| 免费黄网站久久成人精品| 亚洲精品色激情综合| 国产成人午夜福利电影在线观看| 日韩av在线大香蕉| 你懂的网址亚洲精品在线观看| 六月丁香七月| 国内少妇人妻偷人精品xxx网站| 啦啦啦啦在线视频资源| 国产成人免费观看mmmm| 赤兔流量卡办理| 一区二区三区四区激情视频| 美女脱内裤让男人舔精品视频| 国产人妻一区二区三区在| 少妇人妻一区二区三区视频| 免费黄频网站在线观看国产| 我要看日韩黄色一级片| 亚洲成色77777| 久久久久精品性色| av国产免费在线观看| 十八禁国产超污无遮挡网站| 99九九线精品视频在线观看视频| 色网站视频免费| 熟女人妻精品中文字幕| 在线观看一区二区三区| 亚洲丝袜综合中文字幕| av在线播放精品| 免费观看av网站的网址| 啦啦啦啦在线视频资源| 九九久久精品国产亚洲av麻豆| 亚洲丝袜综合中文字幕| 免费观看无遮挡的男女| 欧美精品一区二区大全| 成人亚洲精品av一区二区| 777米奇影视久久| 伦精品一区二区三区| 久久久久九九精品影院| 国产av在哪里看| 高清在线视频一区二区三区| 淫秽高清视频在线观看| 国产高清不卡午夜福利| 免费在线观看成人毛片| 亚洲av男天堂| 亚洲精品一二三| 欧美成人一区二区免费高清观看| 国产黄色免费在线视频| 精品久久久久久久久久久久久| av一本久久久久| 亚洲成色77777| 欧美xxxx性猛交bbbb| 搡女人真爽免费视频火全软件| 亚洲三级黄色毛片| 嫩草影院新地址| 大话2 男鬼变身卡| 欧美日韩亚洲高清精品| 国产国拍精品亚洲av在线观看| 少妇熟女aⅴ在线视频| 五月天丁香电影| 精华霜和精华液先用哪个| 国产色婷婷99| 日韩强制内射视频| 欧美激情国产日韩精品一区| 国产高清有码在线观看视频| 黄色一级大片看看| 视频中文字幕在线观看| 少妇高潮的动态图| 一区二区三区乱码不卡18| 日韩国内少妇激情av| 午夜老司机福利剧场| 永久免费av网站大全| 麻豆av噜噜一区二区三区| 男女下面进入的视频免费午夜| 亚洲不卡免费看| 国产精品一及| 成人午夜精彩视频在线观看| 成人漫画全彩无遮挡| 九草在线视频观看| 一区二区三区乱码不卡18| 亚洲精品aⅴ在线观看| 99热全是精品| 3wmmmm亚洲av在线观看| 老女人水多毛片| 国产成人福利小说| 欧美日本视频| 综合色丁香网| 少妇高潮的动态图| 毛片女人毛片| freevideosex欧美| 亚洲精品456在线播放app| 26uuu在线亚洲综合色| 中文字幕人妻熟人妻熟丝袜美| 在线免费观看的www视频| 国产探花在线观看一区二区| 国产精品精品国产色婷婷| 国产精品.久久久| 成人无遮挡网站| 日本免费在线观看一区| 简卡轻食公司| 2021天堂中文幕一二区在线观| 欧美高清性xxxxhd video| 亚洲av成人av| 欧美潮喷喷水| 最近视频中文字幕2019在线8| 少妇高潮的动态图| 午夜老司机福利剧场| 男的添女的下面高潮视频| 亚洲精品成人av观看孕妇| 黄片wwwwww| 国产 一区 欧美 日韩| 国产成人福利小说| 久久久久免费精品人妻一区二区| 最近中文字幕高清免费大全6| 网址你懂的国产日韩在线| 日韩欧美精品免费久久| 免费观看在线日韩| 久久精品夜色国产| 欧美日韩在线观看h| 高清毛片免费看| 免费不卡的大黄色大毛片视频在线观看 | 国产午夜精品论理片| 精品久久久久久久末码| 成人毛片a级毛片在线播放| 建设人人有责人人尽责人人享有的 | 少妇的逼好多水| 日日摸夜夜添夜夜添av毛片| 男女视频在线观看网站免费| 欧美日韩一区二区视频在线观看视频在线 | 亚洲自拍偷在线| 亚洲精品乱码久久久v下载方式| 国产成人精品婷婷| 亚洲性久久影院| 美女内射精品一级片tv| 大话2 男鬼变身卡| 亚洲一区高清亚洲精品| 久久久成人免费电影| 午夜免费观看性视频| 一区二区三区四区激情视频| 免费观看av网站的网址| 亚洲精品成人av观看孕妇| 日韩精品有码人妻一区| 国产一区有黄有色的免费视频 | 51国产日韩欧美| 国产精品福利在线免费观看| 日韩欧美国产在线观看| 综合色av麻豆|