• <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)用
    国产成人精品婷婷| 国产精品99久久99久久久不卡 | 亚洲精品乱久久久久久| 18禁裸乳无遮挡免费网站照片| 在线亚洲精品国产二区图片欧美 | 久久久a久久爽久久v久久| 在线精品无人区一区二区三 | 特大巨黑吊av在线直播| 欧美人与善性xxx| 1000部很黄的大片| 国产亚洲91精品色在线| 五月天丁香电影| 久久精品国产亚洲av天美| 高清视频免费观看一区二区| 哪个播放器可以免费观看大片| 五月玫瑰六月丁香| 一本一本综合久久| tube8黄色片| 插逼视频在线观看| 搡女人真爽免费视频火全软件| 久久99热这里只频精品6学生| 久久精品国产亚洲网站| 亚洲国产精品999| 亚洲精品自拍成人| 色哟哟·www| 久久久色成人| 久久精品熟女亚洲av麻豆精品| 成年美女黄网站色视频大全免费 | 久久这里有精品视频免费| 在线精品无人区一区二区三 | 亚洲精品国产av成人精品| av线在线观看网站| 久久女婷五月综合色啪小说| 日日撸夜夜添| 国产精品嫩草影院av在线观看| 中文在线观看免费www的网站| 少妇高潮的动态图| 精品一区二区免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产成人一精品久久久| 亚洲经典国产精华液单| 亚洲av成人精品一二三区| 女人十人毛片免费观看3o分钟| 晚上一个人看的免费电影| 内射极品少妇av片p| 街头女战士在线观看网站| 欧美国产精品一级二级三级 | 国产伦精品一区二区三区四那| 两个人的视频大全免费| 成人特级av手机在线观看| 国产成人aa在线观看| 国产欧美日韩精品一区二区| 日产精品乱码卡一卡2卡三| 亚洲精品国产色婷婷电影| 51国产日韩欧美| 秋霞伦理黄片| 国产精品偷伦视频观看了| 久久影院123| 日韩一区二区三区影片| 免费看不卡的av| 亚洲精品日韩在线中文字幕| 综合色丁香网| 亚洲国产成人一精品久久久| 国产熟女欧美一区二区| 国产亚洲91精品色在线| 久久久久性生活片| 人人妻人人爽人人添夜夜欢视频 | 亚洲av男天堂| 国产精品嫩草影院av在线观看| 成人综合一区亚洲| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va| 一级黄片播放器| 亚洲欧美成人精品一区二区| 亚洲精品日韩av片在线观看| 一级毛片我不卡| 精品亚洲乱码少妇综合久久| 成人毛片60女人毛片免费| 久久久久久伊人网av| 好男人视频免费观看在线| 国产乱人偷精品视频| 人妻制服诱惑在线中文字幕| 我的女老师完整版在线观看| 国产精品久久久久久精品电影小说 | 大陆偷拍与自拍| 国产伦精品一区二区三区视频9| av国产久精品久网站免费入址| av线在线观看网站| 欧美精品一区二区大全| 一级二级三级毛片免费看| 欧美3d第一页| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 国产人妻一区二区三区在| 精品一品国产午夜福利视频| 国产真实伦视频高清在线观看| 久久久久久久精品精品| 亚洲精品乱码久久久v下载方式| 亚洲成人av在线免费| 五月玫瑰六月丁香| 免费黄色在线免费观看| 久久av网站| 婷婷色综合大香蕉| 日本一二三区视频观看| 少妇的逼水好多| 91aial.com中文字幕在线观看| 啦啦啦中文免费视频观看日本| 国产欧美日韩一区二区三区在线 | 99国产精品免费福利视频| av国产免费在线观看| 少妇猛男粗大的猛烈进出视频| av黄色大香蕉| 永久免费av网站大全| 夫妻午夜视频| 超碰97精品在线观看| 岛国毛片在线播放| 亚洲国产欧美在线一区| 国产69精品久久久久777片| 99re6热这里在线精品视频| 少妇被粗大猛烈的视频| 国产av国产精品国产| 亚洲av在线观看美女高潮| 免费观看在线日韩| 国产精品精品国产色婷婷| 极品教师在线视频| 国产免费福利视频在线观看| 偷拍熟女少妇极品色| 国产91av在线免费观看| 亚洲精品aⅴ在线观看| 国产亚洲最大av| 国产精品伦人一区二区| 成年美女黄网站色视频大全免费 | av又黄又爽大尺度在线免费看| 乱系列少妇在线播放| 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 欧美xxxx黑人xx丫x性爽| 国产一区有黄有色的免费视频| 国产女主播在线喷水免费视频网站| 91久久精品国产一区二区三区| 国产精品久久久久久精品古装| 看十八女毛片水多多多| 欧美另类一区| 欧美丝袜亚洲另类| 色网站视频免费| 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 老司机影院毛片| 啦啦啦在线观看免费高清www| 最近的中文字幕免费完整| 日韩在线高清观看一区二区三区| 免费播放大片免费观看视频在线观看| 亚洲精品日本国产第一区| 亚洲精华国产精华液的使用体验| 久久99热这里只有精品18| 五月天丁香电影| 久久av网站| 久久99热这里只有精品18| 亚洲国产精品成人久久小说| 亚洲性久久影院| 草草在线视频免费看| 日本欧美国产在线视频| 日本黄色日本黄色录像| 日本猛色少妇xxxxx猛交久久| 色视频www国产| 久久久久久久久久久免费av| 国产男人的电影天堂91| 亚洲av不卡在线观看| 日本-黄色视频高清免费观看| 亚洲精品国产色婷婷电影| 久久久久久久国产电影| 日本-黄色视频高清免费观看| 激情 狠狠 欧美| 国产久久久一区二区三区| 午夜福利在线观看免费完整高清在| 成人特级av手机在线观看| 狠狠精品人妻久久久久久综合| 一级毛片黄色毛片免费观看视频| 黄色怎么调成土黄色| av在线老鸭窝| 亚洲欧美日韩无卡精品| a 毛片基地| 国产探花极品一区二区| 国产国拍精品亚洲av在线观看| 久久99热6这里只有精品| h日本视频在线播放| 精品久久久久久久久av| 欧美成人精品欧美一级黄| 日本-黄色视频高清免费观看| 久久影院123| 极品少妇高潮喷水抽搐| 国产真实伦视频高清在线观看| 婷婷色综合www| 亚洲欧美成人精品一区二区| 777米奇影视久久| 亚洲av综合色区一区| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| av视频免费观看在线观看| 最新中文字幕久久久久| 国产乱来视频区| 国产高清不卡午夜福利| 观看美女的网站| 国产大屁股一区二区在线视频| 亚洲精品久久久久久婷婷小说| av卡一久久| 男女无遮挡免费网站观看| 日本色播在线视频| 免费黄色在线免费观看| 大话2 男鬼变身卡| 亚洲aⅴ乱码一区二区在线播放| 男女无遮挡免费网站观看| 少妇被粗大猛烈的视频| 国产精品一二三区在线看| 亚洲成人av在线免费| 国产精品精品国产色婷婷| 国产大屁股一区二区在线视频| 一区二区三区精品91| 性色avwww在线观看| 欧美日本视频| 国产国拍精品亚洲av在线观看| 亚洲av不卡在线观看| 男的添女的下面高潮视频| 少妇猛男粗大的猛烈进出视频| 国产av一区二区精品久久 | 18禁裸乳无遮挡免费网站照片| 97超碰精品成人国产| av国产免费在线观看| 亚洲国产色片| 国产白丝娇喘喷水9色精品| 欧美国产精品一级二级三级 | 最近中文字幕2019免费版| 久久ye,这里只有精品| 九九爱精品视频在线观看| 寂寞人妻少妇视频99o| 一区二区三区四区激情视频| 精品视频人人做人人爽| 国产精品秋霞免费鲁丝片| 久久人人爽人人片av| 人体艺术视频欧美日本| 国产成人免费无遮挡视频| h视频一区二区三区| 亚洲欧美成人综合另类久久久| 51国产日韩欧美| 久久久久久久久久人人人人人人| 黄色怎么调成土黄色| 美女中出高潮动态图| 最新中文字幕久久久久| 嘟嘟电影网在线观看| 人妻夜夜爽99麻豆av| 国产精品久久久久久av不卡| 99re6热这里在线精品视频| 2018国产大陆天天弄谢| 亚洲av在线观看美女高潮| 日本与韩国留学比较| 99久国产av精品国产电影| 免费大片黄手机在线观看| 精品久久国产蜜桃| 国产成人a区在线观看| 少妇丰满av| 在线免费观看不下载黄p国产| 成人亚洲欧美一区二区av| 成人18禁高潮啪啪吃奶动态图 | 日本一二三区视频观看| av福利片在线观看| 卡戴珊不雅视频在线播放| 国产淫语在线视频| 国产一区亚洲一区在线观看| 老熟女久久久| 3wmmmm亚洲av在线观看| 嘟嘟电影网在线观看| 男女啪啪激烈高潮av片| 色婷婷久久久亚洲欧美| 成人亚洲精品一区在线观看 | 国产精品麻豆人妻色哟哟久久| 深夜a级毛片| 国产精品久久久久久精品电影小说 | 国产男人的电影天堂91| 97在线人人人人妻| 亚洲av成人精品一二三区| 热re99久久精品国产66热6| 日本黄大片高清| 欧美xxxx黑人xx丫x性爽| 内射极品少妇av片p| 欧美少妇被猛烈插入视频| 中文在线观看免费www的网站| 美女主播在线视频| 韩国av在线不卡| av在线蜜桃| 欧美xxxx黑人xx丫x性爽| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站| 亚洲人成网站高清观看| 最近最新中文字幕免费大全7| 国产v大片淫在线免费观看| 天天躁日日操中文字幕| 99国产精品免费福利视频| 精品一区在线观看国产| 小蜜桃在线观看免费完整版高清| 男女边吃奶边做爰视频| 99久久精品一区二区三区| 国产精品人妻久久久影院| 国产久久久一区二区三区| 亚洲熟女精品中文字幕| 美女视频免费永久观看网站| 亚洲欧美日韩无卡精品| 亚洲精品中文字幕在线视频 | 亚洲精品,欧美精品| 国产精品一区www在线观看| 一级毛片黄色毛片免费观看视频| 插阴视频在线观看视频| 大陆偷拍与自拍| 久久 成人 亚洲| 欧美人与善性xxx| 国产亚洲最大av| 成人一区二区视频在线观看| 国产精品免费大片| 国产欧美日韩精品一区二区| 一级a做视频免费观看| 精品亚洲成国产av| 日日撸夜夜添| 丝袜喷水一区| 亚洲精品色激情综合| 97在线人人人人妻| 国产精品一及| videos熟女内射| 视频区图区小说| 蜜桃亚洲精品一区二区三区| 亚洲精品日韩av片在线观看| 欧美成人精品欧美一级黄| 国产午夜精品久久久久久一区二区三区| 亚洲欧美日韩另类电影网站 | 国模一区二区三区四区视频| 亚洲av二区三区四区| 免费观看无遮挡的男女| 亚洲色图av天堂| 国产av国产精品国产| 99国产精品免费福利视频| 成人18禁高潮啪啪吃奶动态图 | 青春草国产在线视频| 亚洲精品亚洲一区二区| 日韩一区二区三区影片| 国产精品成人在线| 一级av片app| 久久99蜜桃精品久久| 一级黄片播放器| 午夜日本视频在线| 亚洲精品第二区| 91午夜精品亚洲一区二区三区| 久久久久精品性色| av网站免费在线观看视频| 日韩亚洲欧美综合| 国产女主播在线喷水免费视频网站| av不卡在线播放| 欧美日韩一区二区视频在线观看视频在线| 六月丁香七月| 最近中文字幕2019免费版| 欧美成人一区二区免费高清观看| 欧美日韩综合久久久久久| 一区二区三区免费毛片| 最近中文字幕2019免费版| 成人无遮挡网站| 黑丝袜美女国产一区| 最近中文字幕高清免费大全6| 欧美人与善性xxx| 青青草视频在线视频观看| 91精品一卡2卡3卡4卡| 亚洲精品国产成人久久av| 少妇的逼水好多| 精品一区在线观看国产| 亚洲性久久影院| 在线观看美女被高潮喷水网站| 亚洲精品一二三| 在线观看国产h片| 欧美老熟妇乱子伦牲交| 国产精品女同一区二区软件| 在线精品无人区一区二区三 | 你懂的网址亚洲精品在线观看| 日韩欧美精品免费久久| 欧美日韩综合久久久久久| 精品久久久精品久久久| 午夜福利在线观看免费完整高清在| 一二三四中文在线观看免费高清| 99久久精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 秋霞伦理黄片| 精品国产乱码久久久久久小说| 日韩亚洲欧美综合| 国产伦在线观看视频一区| 国产精品人妻久久久影院| 亚洲图色成人| 国产亚洲5aaaaa淫片| 美女内射精品一级片tv| 欧美变态另类bdsm刘玥| 久久久久久久精品精品| 国产精品一二三区在线看| 联通29元200g的流量卡| 国产精品不卡视频一区二区| 久久精品久久久久久噜噜老黄| 黄色视频在线播放观看不卡| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 超碰97精品在线观看| 22中文网久久字幕| 午夜福利视频精品| 97超碰精品成人国产| 七月丁香在线播放| 妹子高潮喷水视频| 精品国产露脸久久av麻豆| 成人18禁高潮啪啪吃奶动态图 | 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频| 午夜免费观看性视频| 色婷婷av一区二区三区视频| 国产精品国产三级国产专区5o| 九九在线视频观看精品| 亚洲色图综合在线观看| 国产在线男女| 久久99热这里只有精品18| 天堂俺去俺来也www色官网| 九九在线视频观看精品| 亚洲,一卡二卡三卡| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 建设人人有责人人尽责人人享有的 | 亚洲av二区三区四区| 欧美性感艳星| 99久久精品一区二区三区| 国产精品一及| 欧美成人精品欧美一级黄| 亚洲国产高清在线一区二区三| 久久久久久久精品精品| 99久久精品一区二区三区| 99热网站在线观看| 91精品国产九色| 日韩一本色道免费dvd| 国产在视频线精品| 国产在线免费精品| a级毛片免费高清观看在线播放| 国产乱人偷精品视频| 国产真实伦视频高清在线观看| 简卡轻食公司| 岛国毛片在线播放| 人人妻人人添人人爽欧美一区卜 | 午夜福利在线观看免费完整高清在| 校园人妻丝袜中文字幕| 联通29元200g的流量卡| videossex国产| 韩国av在线不卡| 成人毛片a级毛片在线播放| 久久 成人 亚洲| 联通29元200g的流量卡| 日韩三级伦理在线观看| 久久精品国产a三级三级三级| 国产精品成人在线| 亚洲人与动物交配视频| 亚洲电影在线观看av| 国产免费视频播放在线视频| 久久久久久久久久成人| 激情五月婷婷亚洲| 国产高清三级在线| 日韩精品有码人妻一区| 久久久久性生活片| 久久久久久九九精品二区国产| 午夜福利影视在线免费观看| 亚洲中文av在线| 国产高清国产精品国产三级 | 国产精品人妻久久久影院| 内地一区二区视频在线| 精品久久久久久久久亚洲| 在线免费观看不下载黄p国产| 熟女电影av网| 91精品国产国语对白视频| 91狼人影院| 舔av片在线| 国产av一区二区精品久久 | 人妻少妇偷人精品九色| 国产免费福利视频在线观看| 国精品久久久久久国模美| 91精品一卡2卡3卡4卡| 成人特级av手机在线观看| 成人一区二区视频在线观看| 亚洲欧美成人综合另类久久久| 欧美成人精品欧美一级黄| 91久久精品电影网| 99久国产av精品国产电影| 美女脱内裤让男人舔精品视频| 18禁在线无遮挡免费观看视频| 日本黄色日本黄色录像| 天堂8中文在线网| 久久毛片免费看一区二区三区| 国产精品三级大全| 久久久久久久亚洲中文字幕| 国产精品秋霞免费鲁丝片| 欧美 日韩 精品 国产| 免费观看a级毛片全部| 久久人人爽人人片av| 日韩人妻高清精品专区| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| 日韩成人av中文字幕在线观看| 日韩在线高清观看一区二区三区| 有码 亚洲区| 在线观看免费视频网站a站| 中文精品一卡2卡3卡4更新| 最近最新中文字幕大全电影3| 日日撸夜夜添| 在线观看人妻少妇| 日本vs欧美在线观看视频 | 一个人免费看片子| 99热网站在线观看| 狠狠精品人妻久久久久久综合| 一区在线观看完整版| 日本黄色片子视频| 精品人妻偷拍中文字幕| 国产淫片久久久久久久久| 国产高潮美女av| 麻豆成人午夜福利视频| 成人亚洲精品一区在线观看 | 久久精品国产鲁丝片午夜精品| 视频区图区小说| 免费黄网站久久成人精品| 亚洲国产av新网站| 精品亚洲乱码少妇综合久久| 菩萨蛮人人尽说江南好唐韦庄| 国产深夜福利视频在线观看| 99热国产这里只有精品6| 久久久久久久大尺度免费视频| 欧美bdsm另类| 特大巨黑吊av在线直播| 亚洲精品日韩av片在线观看| 欧美成人精品欧美一级黄| 欧美日韩精品成人综合77777| 狠狠精品人妻久久久久久综合| 国产一区二区在线观看日韩| 中文字幕制服av| 久久久久精品久久久久真实原创| 亚洲人成网站在线播| 亚洲怡红院男人天堂| 九九爱精品视频在线观看| 亚洲国产av新网站| 久久久a久久爽久久v久久| 一本—道久久a久久精品蜜桃钙片| 在线观看人妻少妇| 亚洲自偷自拍三级| 深爱激情五月婷婷| 蜜桃在线观看..| 中国国产av一级| 久久久久国产网址| 卡戴珊不雅视频在线播放| 久久久久久久久久人人人人人人| 丰满少妇做爰视频| 日本免费在线观看一区| 亚洲四区av| 人妻夜夜爽99麻豆av| 亚洲国产精品成人久久小说| 22中文网久久字幕| videossex国产| 久久人妻熟女aⅴ| 在线免费十八禁| 婷婷色麻豆天堂久久| 狂野欧美白嫩少妇大欣赏| 亚洲欧美精品专区久久| 香蕉精品网在线| 久久精品久久久久久噜噜老黄| 精品亚洲乱码少妇综合久久| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品50| 久久人人爽人人爽人人片va| 国产av一区二区精品久久 | 国产精品伦人一区二区| 婷婷色综合大香蕉| 久久久久久久久久久丰满| 中国三级夫妇交换| 日韩视频在线欧美| 伦理电影免费视频| 亚洲欧美成人精品一区二区| 国产黄片视频在线免费观看| 网址你懂的国产日韩在线| 亚洲国产最新在线播放| 人妻夜夜爽99麻豆av| av女优亚洲男人天堂| 亚洲av中文av极速乱| 国产69精品久久久久777片| 国内揄拍国产精品人妻在线| 精品人妻一区二区三区麻豆| 日韩中文字幕视频在线看片 | 国产成人精品婷婷| 免费看不卡的av| 天美传媒精品一区二区| 亚洲经典国产精华液单| 国内精品宾馆在线| 中文字幕av成人在线电影| 国产69精品久久久久777片| 欧美日韩精品成人综合77777| 国产乱人偷精品视频| 大片电影免费在线观看免费| av专区在线播放| 在线观看美女被高潮喷水网站| 久久鲁丝午夜福利片| 美女福利国产在线 | 国产精品国产三级国产av玫瑰| 成人一区二区视频在线观看| 在线观看国产h片| 国产欧美亚洲国产| 一区二区三区乱码不卡18| 国产高潮美女av| 日本黄色片子视频| 国产伦精品一区二区三区视频9| 免费少妇av软件| 欧美性感艳星| 女性生殖器流出的白浆| 99久久精品国产国产毛片| 男女啪啪激烈高潮av片| 亚洲av不卡在线观看|