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

    天宮2號POLAR探測器的低能X射線在軌定標(biāo)?

    2020-02-12 08:05:44HAJDASWojtek劉四明李友平
    天文學(xué)報 2020年1期
    關(guān)鍵詞:耀斑定標(biāo)高能

    張 平 HAJDAS Wojtek 劉四明 蘇 楊 李友平 陳 維

    (1 中國科學(xué)院紫金山天文臺南京210033)

    (2 中國科學(xué)院暗物質(zhì)和空間天文重點(diǎn)實(shí)驗(yàn)室南京210033)

    (3 中國科學(xué)院大學(xué)北京100049)

    (4 Paul Scherrer Institut Villigen 5232)

    1 背景介紹

    X/γ射線作為高能天體物理研究的一個窗口, 由于其輻射機(jī)制和局部等離子環(huán)境以及電磁場特性密切相關(guān), 因此是研究極端天體物理環(huán)境下粒子加速機(jī)制的重要手段.一般而言我們通過測量來自宇宙天體輻射的光子的能量、流強(qiáng)和方位以及偏振特性來全面了解輻射源的物理特性.理論研究表明, 大部分天體輻射產(chǎn)生的硬X射線具有一定的偏振特性, 其偏振特性依賴于輻射機(jī)制、局地的磁場結(jié)構(gòu)、高能粒子分布等輻射源的性質(zhì).在高能天體物理領(lǐng)域, 由于高能光子流量低, 背景復(fù)雜, 并且受限于探測器技術(shù), 大部分探測器主要探測輻射的能量、流強(qiáng)和方位, 偏振測量效果不理想, 誤差較大.

    太陽是離我們最近的高能輻射源, 其高能輻射主要來自于太陽耀斑和日冕物質(zhì)拋射過程中的非熱輻射過程.耀斑作為最劇烈的太陽活動, 其磁重聯(lián)過程可以在102–103s內(nèi)釋放~1032–1033erg的能量[1–2], 大量非熱電子在熱等離子體背景和電磁場中被加速, 通過電子的軔致輻射過程產(chǎn)生 10 keV的硬X射線輻射.一般認(rèn)為其非熱輻射的偏振度約為20%, 并且觀測上如果要區(qū)分不同的輻射加速模型, 偏振度測量精度至少要高于2%[3].近年來得益于多波段的成像成譜觀測, 對耀斑過程中的熱等離子體和非熱電子的加速有了較為精確的觀測, 但是由于觀測的投影效應(yīng), 我們依然無法準(zhǔn)確估測日冕硬X射線源的加速區(qū)熱等離子體的結(jié)構(gòu)和非熱電子的投擲角分布[4].

    太陽耀斑的第1次硬X射線偏振測量來自于前蘇聯(lián)的實(shí)驗(yàn), 其給出的偏振度大概在2%–3%到20%–40%.Zhitnik等[5]利用太陽活動的近地復(fù)雜軌道觀測衛(wèi)星(CORONASF)上的太陽分光偏振儀(SPR-N)探測到的25個樣本給出的偏振度上限在10–100 keV范圍內(nèi)為8%–40% (3σ), 不能對耀斑的粒子加速機(jī)制模型進(jìn)行有效的區(qū)分.RHESSI是主要用于觀測太陽硬X射線和γ射線輻射特性的衛(wèi)星[6], 雖然其攜帶了一個Be散射探測器可以用來重建硬X射線調(diào)制曲線, 實(shí)現(xiàn)測量來自伽馬暴(GRB)和太陽大耀斑的偏振, 但是對樣本的分析表明由于過高的本底和較小的有效散射面積, 沒有辦法有效區(qū)分偏振事件和背景事件[7–9].Suarez-Garcia等[10]對6個X級和1個M級太陽耀斑的RHESSI數(shù)據(jù)的偏振特性的研究給出大耀斑在100–350 keV能量區(qū)間1σ置信度內(nèi)的硬X射線偏振度約為2%–54%, 但是其誤差高達(dá)10%–26%, 不能作為有效的偏振探測器.

    POLAR是由中國、瑞士和波蘭合作研制, 搭載在天宮2號實(shí)驗(yàn)室上的一個γ射線偏振儀, 于2016年9月成功發(fā)射, 主要用來測量50–500 keV能段γ射線的線偏振度和線偏振方位角.其由1600個塑料閃爍體棒組成25個模塊(每個模塊為8×8的塑料閃爍體陣列),采用的是原子序數(shù)較低的塑料閃爍體材料.其測量原理是高能光子在低系數(shù)的塑料閃爍體上發(fā)生康普頓散射的效率較高, 在發(fā)生康普頓散射后, 反沖電子使閃爍棒發(fā)生退激發(fā), 產(chǎn)生熒光; 而散射光子則在其他棒上發(fā)生2次康普頓散射或者光電吸收, 直到散射光子的能量被完全沉積或者逃逸出探測器, 通過收集熒光得到探測器棒中所沉積的能量.基于重建最大光子沉積能量的棒的位置的統(tǒng)計(jì)分布可以給出源的γ光子的線偏振度和方向[11](圖1).POLAR探測的視場大概為1/3天區(qū), 并且具有~200 cm2的有效探測面積和~0.35的調(diào)制因子, 使其可以對GRB和太陽耀斑的硬X射線的偏振度實(shí)現(xiàn)精確測量.并且對于X級耀斑, 基于POLAR的蒙特卡洛模擬顯示, POLAR的偏振度測量精度可以高于5% (見文獻(xiàn)[12]).

    本文利用POLAR的在軌耀斑觀測數(shù)據(jù), 結(jié)合RHESSI的能譜觀測, 給出在耀斑期間峰值附近能譜的演化, 將這些能譜作為POLAR觀測的蒙特卡洛模擬的輸入, 通過擬合模擬結(jié)果和在軌觀測的結(jié)果對POLAR探測器的低能光子探測進(jìn)行相對定標(biāo), 并把有關(guān)結(jié)果與在軌和地面標(biāo)定的高能定標(biāo)結(jié)果做比較, 給出相應(yīng)的討論和分析.

    圖1 POLAR原理示意圖.星號代表在軌定標(biāo)源22Na, 淺藍(lán)色棒代表塑料閃爍體棒, 紅線代表光子在探測器中的徑跡.Fig.1 A cartoon for POLAR’s detection mechanism.Stars indicate the locations of the in-orbit calibration sources 22Na, light blue bars represent the plastic scintillator bars, and red lines show a photon trajectory in the detector.

    2 POLAR的低能X射線在軌相對標(biāo)定方法

    POLAR塑料閃爍體材料的能量分辨率參數(shù)已經(jīng)由實(shí)驗(yàn)室測定[13], 基于Dietze等[14]得到的經(jīng)驗(yàn)公式, 我們估測在小于50 keV的低能區(qū)塑料閃爍體的能量分辨率~20%, 在小于10 keV則低于50%.對于高能的光子數(shù)據(jù), 可以使用在軌的22Na源事例, 通過擬合康普頓邊的方法得到對應(yīng)的能量刻度因子[15].POLAR工作期間其高能在軌和地面定標(biāo)顯示POLAR的能量標(biāo)度在實(shí)驗(yàn)室、模擬和在軌差別很小, 儀器表現(xiàn)符合預(yù)期[16].其能量刻度因子和探測器模塊的高壓設(shè)置成線性關(guān)系, 而探測器溫度對探測器的能量刻度因子影響則較小.對于太陽耀斑來說, 因?yàn)橛^測到的耀斑都是小耀斑, 絕大部分耀斑的X射線光子都是源自熱等離子體的軔致輻射.低能光子在探測器里的康普頓散射效率較低, 閃爍體在低能部分的能量分辨率也較差, 考慮到探測器本身在低能端的非線性特性, 低能端無法直接采用在軌高能部分的標(biāo)定, 而且POLAR在地面也并未做有效的低能定標(biāo).我們要得到小耀斑的準(zhǔn)確能譜信息, 則需要找到一個合適的低能相對定標(biāo)源給出其相對的定標(biāo).在同時段的硬X射線觀測設(shè)備RHESSI能夠非常好地給出觀測到的太陽耀斑能譜信息, 并且在低能區(qū)RHESSI具有很高的分辨率.因此我們可以把RHESSI的能譜觀測作為POLAR在低能區(qū)相對定標(biāo)的參考, 結(jié)合對POLAR觀測的蒙特卡洛模擬, 來給POLAR進(jìn)行相對定標(biāo).

    為了得到較好的低能相對定標(biāo)示例, 我們需要選擇合適的耀斑事件, POLAR觀測到的耀斑列表如附錄表4所示.考慮到高能定標(biāo)的數(shù)據(jù)時間段為2016年11月19日, 我們主要考慮采用比較接近這個時間段, 并且探測器硬件基本參數(shù)設(shè)置一致的事件作為相對定標(biāo)源.綜合表4以及參考POLAR高能定標(biāo)的時間段[16], 我們選擇在軌時間比較接近, 儀器的高壓設(shè)定和溫度等參數(shù)一致并且同時有POLAR和RHESSI觀測的耀斑作為太陽耀斑的低能定標(biāo)源.空間本底的蒙特卡洛模擬[16]顯示彌漫X射線背景、正電子、中子、電子以及蟹狀星云的背景會產(chǎn)生~10%的背景信號, 那么我們選擇軌道遠(yuǎn)離大西洋異常區(qū)的事件, 有相對較低的軌道背景, 可以得到較高信噪比的耀斑數(shù)據(jù).綜合以上考慮, 我們選擇C7.5耀斑SOL2016112907作為低能相對標(biāo)定源.通過對比在軌高壓設(shè)置(如圖2), 可以看到POLAR耀斑期間和高能在軌定標(biāo)時間段內(nèi)的高壓設(shè)置一致.我們認(rèn)為這個耀斑期間儀器的工作狀態(tài)和高能在軌定標(biāo)的時間段的狀態(tài)是一致的, 并且可以利用RHESSI的觀測作為合理的相對定標(biāo)數(shù)據(jù).

    圖2 耀斑期間和在軌定標(biāo)期間POLAR高壓設(shè)置對比Fig.2 Comparison of POLAR high voltage settings during the flare and in-flight calibration

    圖3給出了SOL2016112907的扣除背景后的POLAR觸發(fā)光變和RHESSI光變, 這些觀測有非常好的對應(yīng)關(guān)系.為了更好地顯示高能光變, RHESSI 12–25 keV和25–50 keV的光變分別增大3倍和30倍.從光變上可以看出這個耀斑在25 keV以上的高能光子流量較小, POLAR觀測的峰值部分和RHESSI 12–25 keV的觀測對應(yīng)得較好.我們同時也給出在POLAR坐標(biāo)系中X射線輻射源的角度坐標(biāo)θ ~ 52?,? ~ 262?(圖3下圖).

    圖4給出在耀斑期間(07:02—07:16 UT) POLAR 25個模塊沉積事例的計(jì)數(shù)分布.POLAR下邊的模塊有比較高的計(jì)數(shù), 左上模塊的觸發(fā)是由07:12 UT左右的GRB事件導(dǎo)致的, 其方位坐標(biāo)為θ ~ 42.6?,? ~ 52.3?, ?值與耀斑相差~ 210?.RHESSI (總是指向太陽觀測)并沒有看到這個輻射源, 所以在我們的耀斑相對定標(biāo)分析中不再分析這個GRB事件.我們主要分析耀斑峰值附近響應(yīng)比較高的下邊的模塊, 特別是模塊14.

    為了得到可靠的定標(biāo)結(jié)果, 在25–50 keV峰值附近我們得到3個20 s的積分時間段07:09:00—07:09:20(T1),07:09:20—07:09:40(T2),07:09:40—07:10:00(T3)的RHESSI能譜.我們使用熱分量(vth)加截斷冪律譜(bpow)模型擬合峰值處硬X射線能譜,圖5為T2時刻的能譜擬合, 表1給出了3個時間段的擬合參數(shù)演化, 可以看到能譜擬合中熱分量的等離子體溫度(T)可達(dá)22 MK, 非熱分量主導(dǎo)了20 keV以上的流量, 其截斷能量(Ebreak)在15 keV左右.非熱冪律能譜非常軟, 譜指數(shù)α2~4.0.

    圖3 耀斑SOL2016112907背景扣除后的POLAR觸發(fā)光變(中間)和RHESSI (上)的觀測對應(yīng), 耀斑25–50 keV流量峰值附近陰影部分為選擇的能譜分析的時間片段; 下圖是耀斑在POLAR探測器坐標(biāo)系的位置變化.Fig.3 Background subtracted lightcurve of flare SOL2016112907 from POLAR trigger data (middle) and RHESSI observation (top), and the shaded areas near the flare 25–50 keV flux peak are selected time slices for spectral analysis; the bottom panel shows the evolution of flare position in the POLAR coordinate system.

    圖4 左圖為耀斑期間(07:02—07:16 UT)扣除背景后的POLAR 25個模塊的計(jì)數(shù)分布; 右圖為整個探測器2維計(jì)數(shù)投影圖, 顏色對應(yīng)于模塊計(jì)數(shù).Fig.4 Left: counts distribution of the 25 modules of POLAR after background subtraction during the flare (07:02—07:16 UT); right: 2D counts distribution of the whole detector.Colors correspond to module counts.

    圖5 耀斑流量峰值附近的RHESSI X射線的能譜擬合, 其中, Data-BK表示減去背景的能譜, Ebreak為截斷能量,EM為熱分量的發(fā)射度, T 為等離子體溫度, α1為低于截斷能量的冪律指數(shù), α2為高于截斷能量的冪律指數(shù).Fig.5 RHESSI X-ray spectral fit near the flux peak, here, Data-BK is the spectrum after background subtraction, Ebreak is the break energy, EM is the thermal emission measure, T is the plasma temperature,α1 is the power law index below break energy, and α2 is the power law index above break energy.

    表1 3個時間段(T1, T2, T3)的能譜擬合結(jié)果Table 1 The spectral fit results of 3 time slices (T1, T2, T3)

    相對于POLAR探測器, 太陽耀斑在這里可以看作是一個點(diǎn)源.表1還給出了耀斑位置在POLAR坐標(biāo)系的變化(θ ~ 52?,? ~ 267?), 可以看出源位置的變化對我們結(jié)果的影響可以忽略.而且閾值設(shè)定在5–40 keV范圍內(nèi), POLAR的探測效率的變化很小, 但是在低能部分塑料閃爍體的能量響應(yīng)表現(xiàn)出明顯的非線性特征[16], 為了盡可能減小非線性效應(yīng)的影響, 我們在使用POLAR GEANT4軟件包進(jìn)行模擬的時候, 入射的γ光子流的能譜范圍選擇為10–300 keV.這樣不僅可以提高模擬效率, 還可以在探測器里有效沉積高能事例.考慮到POLAR探測器在艙外的尺寸為450×450×250 mm3, 模擬中我們在距離探測器300 cm的球面上放置一個半徑40 cm的圓盤, 以表1給定的角度及其對應(yīng)的能譜入射了5×106個γ光子, 以得到足夠的事例和觀測進(jìn)行比較來實(shí)現(xiàn)相對定標(biāo).

    光子在探測器里被反射或吸收會造成不同程度的能量衰減, 而且并非所有的沉積能量都會轉(zhuǎn)變成熒光被收集.在耀斑的譜模擬過程中, GEANT4模擬軟件也考慮了這些相互作用過程, 并且閃爍體棒的能量分辨率參數(shù)已經(jīng)由實(shí)驗(yàn)室測定[15].那么對于同樣的γ光子譜輸入, 穩(wěn)定工作的POLAR探測器在蒙特卡洛模擬中所沉積的能量是可以和實(shí)際觀測到的能量沉積相比較的.但是由于光電倍增管的不均勻響應(yīng)以及探測器本身集成系統(tǒng)的串?dāng)_, 也會造成實(shí)際測量能量的一些非線性響應(yīng)[17–18].為了盡量減小系統(tǒng)以及背景誤差, 我們在具體做相對定標(biāo)的時候, 只選取在耀斑峰值附近大量被激發(fā)的閃爍體棒作為分析對象.在方法驗(yàn)證過程中, 主要是以模塊14閃爍體棒7 (后面簡稱為模塊14棒7)作為示例.在比較窄的能譜范圍內(nèi), 能量刻度參數(shù)是近似線性的.在篩選沉積的有效事例時, 我們參考在軌定標(biāo)對數(shù)據(jù)處理的流程[16,19], 對讀出數(shù)據(jù)的基礎(chǔ)噪音和共模噪音進(jìn)行矯正, 并且剔除科學(xué)數(shù)據(jù)中標(biāo)記的宇宙線事件(當(dāng)某個時刻產(chǎn)生沉積的棒較多或者事例數(shù)過高), 得到耀斑期間的有效數(shù)據(jù), 即我們所處理的事例.由于在分析GRB偏振的時候?qū)τ?jì)數(shù)的精度要求較高, 因此在得到需要分析的數(shù)據(jù)后需要矯正響應(yīng)的非線性影響,再計(jì)算閾值和進(jìn)行探測器整體之間串?dāng)_的修正.定標(biāo)源22Na產(chǎn)生的511 keV光子譜的康普頓邊為340 keV, 擬合其模擬數(shù)字轉(zhuǎn)換器(ADC)得到的道數(shù)(channel)能譜的康普頓邊, 就可以確定探測器的能量標(biāo)度因子.在耀斑事件過程中, 雖然耀斑持續(xù)時間較長,但是大部分硬X射線的光子在50 keV以下, 其能譜沒有一個特征能量的譜峰或康普頓邊,并且考慮到閃爍體棒在低能的分辨率較差, 所以在低能端無法通過簡單的特征康普頓邊做標(biāo)定.在相對定標(biāo)時, 我們通過RHESSI能譜分析已經(jīng)確定了入射到POLAR的光子譜以及方向, 如果探測器的響應(yīng)是穩(wěn)定的, 那么其所沉積的能量(EkeV)譜(EkeV)和經(jīng)過矯正基礎(chǔ)噪音和共模噪音的閃爍體棒上的沉積能量(EADC)譜是可以對應(yīng)的, 即在EADC=fadc2kevEkeV的能量刻度下并且其能量通道刻度因子fadc2kev即是我們所需要得到的.注意這里我們已經(jīng)忽略了探測器串?dāng)_對結(jié)果的影響.

    3 定標(biāo)結(jié)果

    圖6左給出了在這個耀斑期間, POLAR模塊14棒7的基準(zhǔn)噪音.可以看出在耀斑期間POLAR模塊14閃爍棒7的基準(zhǔn)噪音分布可以用高斯分布較好擬合, 其平均值(Mean)~177.7 ADC道數(shù)(ADC channel), 均方差(σ)~34.5 ADC channel.由于基準(zhǔn)噪音的讀出頻率為1 Hz, 每一秒都會有一個基礎(chǔ)噪音的數(shù)據(jù), 那么在相對較短的時間范圍內(nèi)可以認(rèn)為閃爍體棒所在的模塊的系統(tǒng)基準(zhǔn)噪音是相對穩(wěn)定的.在后續(xù)的分析中, 我們將篩選出積分時間區(qū)間內(nèi)的事例, 對每個事例通過減去前一秒的基準(zhǔn)噪音來消除讀出通道的基準(zhǔn)噪音.對扣除基礎(chǔ)噪音后的事例, 我們再去除宇宙線事例以及在宇宙線事例后死時間< 10?4s內(nèi)的事例, 它們主要是宇宙線事件引起的次級響應(yīng), 并非正常的有效事例.并且在單個模塊中一些事例同時會有多個觸發(fā)響應(yīng), 圖6中表明這類事例的平均最大沉積能譜在< 70 ADC channel有一個峰.所以為了盡可能地減少非正常事例, 我們將所有平均最大沉積能量(AVMADC) < 70 ADC channel和宇宙線事例后死時間里的事例都去除, 從而得到沉積在POLAR中的正常事例.篩選出我們所需要的正常事例后, 我們可以給出探測器模塊的整體噪音水平即共模噪音.這里我們篩選出沒有觸發(fā)并且沒有和觸發(fā)事例相鄰的探測器閃爍體棒, 通過讀出平均值的分布來給出單個模塊的共模噪音分布[16,19], 圖6右給出耀斑前背景(07:02:00—07:03:00 UT)、峰值(07:09:00—07:10:00 UT)和耀斑后背景(07:15:00—07:16:00 UT)分別為一分鐘時間段的共模噪音分布.可以看出耀斑前后的背景共模噪音基本一致, 在峰值處由于受耀斑事件的影響有一定偏移.

    圖6 左圖為耀斑期間POLAR模塊14棒7的基準(zhǔn)噪音分布(灰線)及其高斯擬合(黑線); 中間為扣除基準(zhǔn)噪音后平均最大沉積能量的能譜圖(實(shí)線)和扣除宇宙射線之后的次級響應(yīng)的平均最大沉積能譜(虛線), 點(diǎn)線為70 ADC channel分割線; 右圖為模塊14在耀斑前背景(黑線), 峰值(紅線)以及耀斑后背景(藍(lán)線)的共模噪音分布.Fig.6 Left: distribution of the pedestal noise of module 14 bar 7 during the flare (gray line) and its Gaussian fit (black line).Middle: distribution of the averaged maximum energy deposition of all events(solid line) and post cosmic ray events (dashed line), the dotted line stands for 70 ADC channel line.Right: Common noise of module 14 at background before the flare (black line), peak (red line), and background after the flare (blue line).

    圖7左為我們矯正基礎(chǔ)噪音和共模噪音后, 篩選出的正常事例隨時間演化, 可以看出在耀斑期間有明顯的讀出通道增強(qiáng)和事例數(shù)增多, 但是依然有一些和耀斑無關(guān)的背景事例.這些事例無法通過事例的篩選來剔除, 需要從耀斑能譜中直接扣除.圖7右圖給出了對這個閃爍體棒背景的線性擬合, 其斜率(slope) ~ ?0.0001 counts·bar?1·s?1,背景流量較低~ 10 counts·bar?1·s?1.根據(jù)空間背景的蒙特卡洛模擬分析[20]在南大西洋異常區(qū)(SAA)之外的背景主要由彌散的宇宙X射線背景主導(dǎo), 其產(chǎn)生的背景流量~5.4 counts·bar?1·s?1, 另外還有正電子、中子、電子以及蟹狀星云所產(chǎn)生的背景, 總的背景大約會產(chǎn)生~16 counts·bar?1·s?1, 和我們觀測到的背景相當(dāng).圖8左圖顯示耀斑前(07:02:00—07:03:00 UT), 耀斑后(07:15:00—07:16:00 UT)背景能譜沒有明顯變化,我們因此可以結(jié)合擬合曲線扣除我們分析時間段里的非耀斑背景事例譜.

    圖8右圖為我們得到的扣除背景譜后在峰值處(07:09:20—07:09:40 UT)的積分沉積能譜, 即模塊14棒7的未定標(biāo)ADC channel譜.探測器閃爍體棒有一個觸發(fā)閾值, 在能譜上則可以利用誤差函數(shù)擬合低能端ADC channel譜來確定沉積譜的閾值范圍[19].在圖8右圖中可以看到誤差函數(shù)的擬合, 誤差函數(shù)的中心(center) ~192 ADC channel,寬度(width) ~61 ADC channel, 那么可以計(jì)算得到模塊14棒7的通道閾值范圍~278 ADC channel.

    圖7 左圖為耀斑期間模塊14棒7扣除基準(zhǔn)噪音和共模噪音并篩選出正常事件后所有探測到的事件隨時間的演化; 右圖為經(jīng)過事件篩選后模塊14棒7的光變曲線(黑色線)以及其背景的線性擬合(灰色線).Fig.7 Left: temporal evolution of all normal events recorded by module 14 bar 7 after removing pedestal noise and common noise during the flare; right: the flare lightcurve (black line) of module 14 bar 7 after events selection and a linear fit to the background (gray line).

    圖8 左圖為耀斑前(07:02:00—07:03:00 UT)(黑線)和耀斑后(07:15:00—07:16:00 UT)(藍(lán)線)的背景譜; 右圖為模塊14棒7在07:09:20 UT到07:09:40 UT扣除背景譜后的20 s積分沉積譜, 紅色曲線為閾值附近的誤差函數(shù)的擬合.Fig.8 Left: background spectra before the flare (07:02:00—07:03:00 UT) (black line) and after the flare(07:15:00—07:16:00 UT) (blue line).Right: 20 s integration deposition spectrum of module 14 bar 7 from 07:09:20 UT to 07:09:40 UT after background subtracted, the red line is a error function fitting near the threshold.

    通過線性變換EADC= fadc2kevEkeV, 由觀測得到的ADC譜可以和基于RHESSI觀測能譜的蒙特卡洛模擬沉積能譜做比較.這里我們還需要考慮閃爍體棒激發(fā)閾值對沉積能譜的影響, 即利用前面擬合ADC譜得到的誤差函數(shù)對模擬沉積能譜加權(quán).然后通過最小化能量閾值之上的觀測和模擬能譜的χ2給出相應(yīng)的相對定標(biāo)參數(shù)fadc2kev和能量閾值Ethreshold.圖9給出了在峰值處(07:09:20—07:09:40)觀測和模擬的能譜擬合, 其誤差棒為歸一化的統(tǒng)計(jì)誤差, 擬合的χ2為0.25.由Xiao等[16]在軌定標(biāo)的結(jié)果知道模塊14棒7高能端的閾值~ 18.57 keV、fadc2kev~ 14.4 ADC channel · keV?1, 而由RHESSI譜重建的相對定標(biāo)得到的3個時間片段(T1, T2, T3)的標(biāo)定參數(shù), 如表2所示, 擬合得到的閾值~14 keV, 小于高能端的閾值.擬合得到fadc2kev的均值和方差分別為20.5 ADC channel· keV?1和0.07 ADC channel · keV?1.在峰值時間段內(nèi)低能刻度因子較為穩(wěn)定.但是擬合給出的相對定標(biāo)參數(shù)和高能端使用22Na源事例的康普頓邊得到的結(jié)果差別顯著.

    圖9 T2觀測能譜和模擬能譜的比較, 虛線為閾值能量Ethreshold, 藍(lán)色線代表模擬能譜, 紅色線代表觀測標(biāo)定能譜.Fig.9 Comparison of observed spectrum and simulated spectrum for T2, dashed line is the threshold energy Ethreshold, blue line is simulation spectrum, and red line is calibration observation spectrum.

    表2 3個時間段(T1, T2, T3)的定標(biāo)參數(shù)Table 2 Calibration parameters of 3 time intervals (T1, T2, T3)

    應(yīng)用同樣的定標(biāo)方法我們對耀斑流量峰值附近20 s積分時間間隔里能譜峰值計(jì)數(shù)率> 1 counts·s?1(激活)的閃爍體棒進(jìn)行了定標(biāo).表3可以看到不同的棒定標(biāo)給出的能量閾值和通道能量轉(zhuǎn)換因子也不同, 這個結(jié)果和高能端的在軌定標(biāo)結(jié)果類似, 其能量和ADC通道的線性關(guān)系如圖10所示.所有激活的閃爍體棒的定標(biāo)及其和高能端的對比結(jié)果如圖11, 激活的閃爍體棒約為41個, 分布于5個模塊(M4、M9、M14、M19、M24),在POLAR最外面的3排, 約占POLAR整個探測器的塑料閃爍體總數(shù)的~3%, 它們探測到的有效事例和相對定標(biāo)的結(jié)果在不同閃爍體棒上沒有明顯的規(guī)律性, 對比能量閾值低能端(Low threshold)和高能端(High threshold)沒有表現(xiàn)出明顯的差異, 符合較好; 能量通道刻度因子fadc2kev低能端的定標(biāo)(Low factor)對比高能端(High factor)差異同樣沒有表現(xiàn)出明顯的規(guī)律性.這一結(jié)果表明POLAR矩陣式的設(shè)計(jì)導(dǎo)致探測效率依賴源的方向.在這個耀斑期間只有一側(cè)的探測器模塊有較強(qiáng)的響應(yīng), 并且其中被激活的閃爍體棒也主要集中于源方向一側(cè)的棒.

    表3 耀斑期間有顯著信號的閃爍體棒的相對定標(biāo)和在軌定標(biāo)參數(shù)的對比Table 3 Comparison of cross calibration and in-flight calibration parameters in plastic bars with a strong signal during the flare

    圖10 圖中為能量標(biāo)定后, 能量和ADC通道的線性關(guān)系, 紅色為低能端, 灰色為誤差棒, 藍(lán)色為高能端.左圖為模塊14棒7的線性關(guān)系圖; 右圖為所有標(biāo)定的閃爍體棒的線性關(guān)系圖.Fig.10 The linear relation between energy and ADC channel, red is low energy band, gray is error bar,and blue is high energy band.Left: the linear relation of module 14 bar 7; right: the linear relation of all calibration bars.

    圖11 POLAR下3排棒的定標(biāo)參數(shù).上圖為每個棒耀斑期間記錄的有效事例數(shù); 中圖和下圖分別為高低能端定標(biāo)閾值的能量以及對應(yīng)的定標(biāo)的能量通道刻度因子的對比.其中M代表POLAR模塊, 灰色虛線代表不同排分割線.Fig.11 Calibration parameters for the bars in the bottom three rows of POLAR.The top panel gives the total number of events recorded by each bar; the middle panel and the bottom panel are the comparison of threshold energy and conversion factors obtained via calibration in low and high energy ranges,respectively.M stands for POLAR module, and the dashed line is the line between different rows.

    4 總結(jié)與展望

    經(jīng)過半年多的觀測, POLAR成功探測到了多個太陽耀斑, 但是由于探測到的耀斑在高能段的輻射較弱, 無法使用地面或在軌定標(biāo)的參數(shù)對所探測到的耀斑進(jìn)行比較準(zhǔn)確的能譜和輻射偏振測量.我們通過結(jié)合RHESSI的硬X射線能譜觀測, 給出了通過模擬對觀測能譜進(jìn)行相對定標(biāo)的方法.分析表明在非康普頓模式下POLAR塑料閃爍體可以探測到低能的γ光子, 并且由于具有較大的有效面積相對于RHESSI可以探測到更多的高能光子事例.為了和更高能段的地面實(shí)驗(yàn)室、模擬和在軌標(biāo)定結(jié)果做比較, 我們選擇和POLAR在軌標(biāo)定相近的時間段且儀器設(shè)置不變情形下的耀斑事件對POLAR的低能段進(jìn)行相對定標(biāo).在確定耀斑的位置和能譜后, 我們利用RHESSI觀測得到的能譜作為POLAR蒙特卡洛模擬的輸入, 得到相應(yīng)的探測器沉積能譜.通過合理的事例篩選和背景扣除得到的POLAR ADC channel能譜可以和模擬得到的沉積能譜比較得到相應(yīng)探測器棒的能量觸發(fā)閾值和能量通道刻度因子fadc2kev.雖然在耀斑期間這些定標(biāo)參數(shù)比較穩(wěn)定, 但是和高能端的定標(biāo)結(jié)果相比, 兩者有較為顯著的差別.我們的分析還表明陣列式的模塊排列導(dǎo)致探測器之間相互遮擋, 顯著壓低了中間9個模塊的計(jì)數(shù)率.未來探測器設(shè)計(jì)上要盡量避免各模塊相互遮擋.

    附錄 POLAR對太陽耀斑的觀測

    雖然是在太陽活動谷年發(fā)射運(yùn)行, POLAR依然探測到了十幾個太陽耀斑事件,結(jié)合RHESSI硬X射線觀測以及靜止環(huán)境觀測衛(wèi)星(GOES)軟X射線觀測的耀斑列表,我們得到POLAR所探測到的耀斑列表4.從表4中, 可以看到POLAR耀斑列表中探測到的16個耀斑中有11個是C級, 5個是B級耀斑, 都是較小的太陽耀斑.有3個耀斑沒有對應(yīng)的RHESSI觀測, RHESSI有對應(yīng)觀測的耀斑50–300 keV能段的計(jì)數(shù)率都低于102counts/s, 但是POLAR更大的視場和有效探測面積使其能探測到大量的光子.

    表4 POLAR觀測到的耀斑列表Table 4 Solar flares list observed by POLAR

    猜你喜歡
    耀斑定標(biāo)高能
    前方高能!戰(zhàn)機(jī)懟瞼
    我國為世界大豆精準(zhǔn)選種“定標(biāo)”
    基于恒星的電離層成像儀在軌幾何定標(biāo)
    搞笑秀
    意林彩版(2020年3期)2020-11-22 20:54:12
    太陽表面的微型耀斑
    “最強(qiáng)”太陽耀斑爆發(fā)會帶來什么?
    《高能少年團(tuán)》少年 未來可期
    基于角反射器的機(jī)載毫米波云雷達(dá)外定標(biāo)實(shí)驗(yàn)
    4m直徑均勻擴(kuò)展定標(biāo)光源
    Duang!6·18巾幗館前方高能
    海峽姐妹(2015年7期)2015-02-27 15:11:35
    亚洲中文字幕日韩| 亚洲avbb在线观看| 欧美精品人与动牲交sv欧美| 亚洲一区二区三区欧美精品| 欧美精品一区二区免费开放| 校园春色视频在线观看| 咕卡用的链子| 建设人人有责人人尽责人人享有的| 亚洲av成人不卡在线观看播放网| 国产又色又爽无遮挡免费看| 欧美黑人欧美精品刺激| 亚洲av成人av| 黑人巨大精品欧美一区二区蜜桃| 两个人免费观看高清视频| 久久人人爽av亚洲精品天堂| 欧美精品人与动牲交sv欧美| 亚洲五月天丁香| 中文字幕人妻丝袜制服| 叶爱在线成人免费视频播放| 交换朋友夫妻互换小说| 亚洲人成电影观看| 一级,二级,三级黄色视频| 亚洲精华国产精华精| 男人的好看免费观看在线视频 | 纯流量卡能插随身wifi吗| 免费在线观看影片大全网站| 亚洲国产欧美日韩在线播放| 久久久久久免费高清国产稀缺| av不卡在线播放| av天堂在线播放| 亚洲熟妇熟女久久| 99精品久久久久人妻精品| 亚洲在线自拍视频| 九色亚洲精品在线播放| 热99久久久久精品小说推荐| 亚洲精华国产精华精| 亚洲avbb在线观看| 午夜两性在线视频| 精品国产美女av久久久久小说| 国产成人精品无人区| 亚洲精品乱久久久久久| 欧美激情久久久久久爽电影 | 婷婷成人精品国产| 国产精品九九99| 麻豆国产av国片精品| 亚洲情色 制服丝袜| 国产精品亚洲一级av第二区| 国产欧美日韩一区二区三| 母亲3免费完整高清在线观看| 色综合欧美亚洲国产小说| 亚洲专区字幕在线| av天堂在线播放| 亚洲性夜色夜夜综合| 黄色片一级片一级黄色片| 脱女人内裤的视频| 亚洲伊人色综图| 老司机影院毛片| 国产精品秋霞免费鲁丝片| 久久国产亚洲av麻豆专区| 91字幕亚洲| 欧美激情极品国产一区二区三区| 国产97色在线日韩免费| 激情在线观看视频在线高清 | 久久人人爽av亚洲精品天堂| 欧美最黄视频在线播放免费 | 久久精品亚洲av国产电影网| 成年人午夜在线观看视频| 国产成人一区二区三区免费视频网站| 中国美女看黄片| 1024视频免费在线观看| 久久这里只有精品19| 99热网站在线观看| 婷婷丁香在线五月| 欧美国产精品一级二级三级| 人人澡人人妻人| 成人永久免费在线观看视频| 又黄又爽又免费观看的视频| 黑人巨大精品欧美一区二区mp4| 日本黄色视频三级网站网址 | 深夜精品福利| 免费在线观看黄色视频的| 欧美最黄视频在线播放免费 | 亚洲国产看品久久| 男人的好看免费观看在线视频 | 色综合婷婷激情| 精品人妻熟女毛片av久久网站| 亚洲成人免费电影在线观看| 欧美日韩亚洲高清精品| 少妇粗大呻吟视频| 日本五十路高清| 亚洲欧美激情在线| 新久久久久国产一级毛片| 成年动漫av网址| 老汉色av国产亚洲站长工具| 国产精品av久久久久免费| 国产精品成人在线| 十分钟在线观看高清视频www| 成年女人毛片免费观看观看9 | 国产精华一区二区三区| 欧美黄色淫秽网站| 日韩欧美在线二视频 | 亚洲成av片中文字幕在线观看| 欧美日韩中文字幕国产精品一区二区三区 | av视频免费观看在线观看| 日韩制服丝袜自拍偷拍| 精品国产乱码久久久久久男人| e午夜精品久久久久久久| 久久精品亚洲精品国产色婷小说| 无遮挡黄片免费观看| 国产精品av久久久久免费| 99久久综合精品五月天人人| 麻豆成人av在线观看| 深夜精品福利| 亚洲一区高清亚洲精品| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看www视频免费| 国产精品免费一区二区三区在线 | 黑人巨大精品欧美一区二区mp4| 男女床上黄色一级片免费看| 成人国语在线视频| 中文亚洲av片在线观看爽 | 国产亚洲精品一区二区www | 午夜福利免费观看在线| 亚洲伊人色综图| 91麻豆精品激情在线观看国产 | 国产男女超爽视频在线观看| 国产免费现黄频在线看| 天天操日日干夜夜撸| 69精品国产乱码久久久| 成人三级做爰电影| 黄色毛片三级朝国网站| 美女国产高潮福利片在线看| 日本五十路高清| 午夜福利,免费看| 成人亚洲精品一区在线观看| 在线观看免费日韩欧美大片| 午夜福利在线观看吧| 亚洲国产欧美一区二区综合| 成年人黄色毛片网站| 高清av免费在线| 我的亚洲天堂| 国产精品久久久久成人av| 建设人人有责人人尽责人人享有的| 身体一侧抽搐| 欧美av亚洲av综合av国产av| 岛国毛片在线播放| 人成视频在线观看免费观看| 热99re8久久精品国产| 在线播放国产精品三级| 中文字幕制服av| 久久精品人人爽人人爽视色| 久久性视频一级片| 亚洲熟女精品中文字幕| 国产亚洲欧美98| 久久ye,这里只有精品| 国产免费男女视频| 女人高潮潮喷娇喘18禁视频| 久久久久国内视频| 好看av亚洲va欧美ⅴa在| 亚洲片人在线观看| 99热国产这里只有精品6| 国产在视频线精品| 国产欧美日韩精品亚洲av| 久99久视频精品免费| 好看av亚洲va欧美ⅴa在| 欧美黑人精品巨大| 日本撒尿小便嘘嘘汇集6| 97人妻天天添夜夜摸| 99精品欧美一区二区三区四区| 在线观看免费午夜福利视频| 无限看片的www在线观看| 两个人免费观看高清视频| 别揉我奶头~嗯~啊~动态视频| 在线观看舔阴道视频| 一区二区日韩欧美中文字幕| 99热网站在线观看| 丰满饥渴人妻一区二区三| 在线永久观看黄色视频| 中文欧美无线码| 黄片播放在线免费| 免费少妇av软件| 一级,二级,三级黄色视频| 国产成人免费无遮挡视频| 波多野结衣av一区二区av| 亚洲黑人精品在线| 色播在线永久视频| 在线免费观看的www视频| 三级毛片av免费| 国产精品99久久99久久久不卡| 久久香蕉国产精品| 午夜成年电影在线免费观看| 午夜影院日韩av| a在线观看视频网站| 色播在线永久视频| 国产91精品成人一区二区三区| 日本一区二区免费在线视频| 麻豆乱淫一区二区| 黄色视频,在线免费观看| 国产高清videossex| 亚洲一区二区三区不卡视频| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 777久久人妻少妇嫩草av网站| 国产一区二区激情短视频| www.精华液| 欧美+亚洲+日韩+国产| 啪啪无遮挡十八禁网站| 亚洲午夜理论影院| 99热只有精品国产| 国产精品 欧美亚洲| 悠悠久久av| 国产有黄有色有爽视频| 午夜免费鲁丝| 国产又色又爽无遮挡免费看| 99久久精品国产亚洲精品| 黑人巨大精品欧美一区二区mp4| 91成年电影在线观看| 法律面前人人平等表现在哪些方面| 日韩免费av在线播放| 啪啪无遮挡十八禁网站| 亚洲自偷自拍图片 自拍| 一本一本久久a久久精品综合妖精| 高清毛片免费观看视频网站 | 99国产精品一区二区三区| 亚洲av成人av| 校园春色视频在线观看| 久久国产乱子伦精品免费另类| 又大又爽又粗| 精品久久久久久久久久免费视频 | 亚洲一区二区三区欧美精品| 老鸭窝网址在线观看| 色综合欧美亚洲国产小说| www日本在线高清视频| 亚洲国产精品sss在线观看 | 亚洲视频免费观看视频| 国产精品久久久久久人妻精品电影| 最新在线观看一区二区三区| 在线观看免费高清a一片| ponron亚洲| 黄片小视频在线播放| 欧美丝袜亚洲另类 | 免费日韩欧美在线观看| 久久精品人人爽人人爽视色| 精品国内亚洲2022精品成人 | 国产麻豆69| 亚洲成人免费电影在线观看| 日韩人妻精品一区2区三区| 看片在线看免费视频| av视频免费观看在线观看| 色婷婷av一区二区三区视频| 999久久久国产精品视频| 一级,二级,三级黄色视频| 搡老熟女国产l中国老女人| 少妇被粗大的猛进出69影院| 十分钟在线观看高清视频www| 老汉色∧v一级毛片| 夜夜爽天天搞| 99国产精品一区二区三区| 大香蕉久久网| 制服人妻中文乱码| 日韩精品免费视频一区二区三区| 日本vs欧美在线观看视频| 超碰97精品在线观看| 久久久久久久精品吃奶| 两性午夜刺激爽爽歪歪视频在线观看 | 中出人妻视频一区二区| 丝袜美足系列| 一级a爱片免费观看的视频| 啪啪无遮挡十八禁网站| 人人澡人人妻人| 欧美精品啪啪一区二区三区| 国产欧美日韩精品亚洲av| 亚洲成a人片在线一区二区| 国产视频一区二区在线看| 捣出白浆h1v1| 亚洲免费av在线视频| 老司机亚洲免费影院| 午夜福利乱码中文字幕| 夜夜爽天天搞| 丰满人妻熟妇乱又伦精品不卡| 变态另类成人亚洲欧美熟女 | 色播在线永久视频| 天天影视国产精品| 亚洲熟妇熟女久久| 日本vs欧美在线观看视频| 韩国av一区二区三区四区| 国产一区有黄有色的免费视频| 中文字幕另类日韩欧美亚洲嫩草| 大片电影免费在线观看免费| 国产亚洲av高清不卡| av有码第一页| 性少妇av在线| 国产精品一区二区精品视频观看| 99在线人妻在线中文字幕 | 国产免费现黄频在线看| 真人做人爱边吃奶动态| 国产亚洲欧美在线一区二区| 精品少妇久久久久久888优播| 日韩欧美一区二区三区在线观看 | 黑人巨大精品欧美一区二区蜜桃| 欧美日韩视频精品一区| 999精品在线视频| 久久精品91无色码中文字幕| 曰老女人黄片| 久久人妻福利社区极品人妻图片| 国产成人免费无遮挡视频| 欧美黑人欧美精品刺激| 国产在线精品亚洲第一网站| 人妻一区二区av| 日韩欧美在线二视频 | 精品久久久久久久久久免费视频 | 亚洲人成77777在线视频| 我的亚洲天堂| 精品乱码久久久久久99久播| 制服人妻中文乱码| 欧美黄色片欧美黄色片| 最新美女视频免费是黄的| 两个人免费观看高清视频| 人妻久久中文字幕网| 人人妻,人人澡人人爽秒播| 亚洲一卡2卡3卡4卡5卡精品中文| 不卡一级毛片| 男人操女人黄网站| 亚洲精品国产区一区二| 1024视频免费在线观看| 一边摸一边做爽爽视频免费| 亚洲五月婷婷丁香| 午夜两性在线视频| 中文欧美无线码| 最新美女视频免费是黄的| 免费在线观看影片大全网站| 别揉我奶头~嗯~啊~动态视频| av不卡在线播放| 欧美日韩亚洲高清精品| 中文字幕最新亚洲高清| 一本大道久久a久久精品| 久久久水蜜桃国产精品网| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久精品古装| 国产xxxxx性猛交| 天天操日日干夜夜撸| 欧美精品av麻豆av| 欧美日韩亚洲高清精品| 俄罗斯特黄特色一大片| 国产成人啪精品午夜网站| 国产99白浆流出| 亚洲精品美女久久久久99蜜臀| 国产极品粉嫩免费观看在线| 久久青草综合色| 久久国产乱子伦精品免费另类| 国产成人精品久久二区二区91| 成熟少妇高潮喷水视频| 免费在线观看影片大全网站| 国产av精品麻豆| av一本久久久久| 黄色 视频免费看| 51午夜福利影视在线观看| 精品第一国产精品| 伊人久久大香线蕉亚洲五| 老汉色∧v一级毛片| a级毛片在线看网站| 在线十欧美十亚洲十日本专区| 久久久久精品人妻al黑| 免费久久久久久久精品成人欧美视频| 国产成人欧美| 亚洲av熟女| 美女福利国产在线| 国产免费av片在线观看野外av| 欧美人与性动交α欧美精品济南到| 成年版毛片免费区| 日韩欧美免费精品| 少妇猛男粗大的猛烈进出视频| 午夜影院日韩av| 别揉我奶头~嗯~啊~动态视频| 亚洲三区欧美一区| 51午夜福利影视在线观看| 9191精品国产免费久久| 国产淫语在线视频| 久久久久久亚洲精品国产蜜桃av| 亚洲三区欧美一区| 交换朋友夫妻互换小说| 色综合婷婷激情| 成年动漫av网址| 高清毛片免费观看视频网站 | 国产淫语在线视频| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品合色在线| 欧美日本中文国产一区发布| 久久久久国产一级毛片高清牌| 欧美亚洲日本最大视频资源| 色94色欧美一区二区| 精品午夜福利视频在线观看一区| 热re99久久精品国产66热6| 国产亚洲精品久久久久5区| 丰满迷人的少妇在线观看| 老熟女久久久| 国产黄色免费在线视频| 咕卡用的链子| 在线免费观看的www视频| 免费在线观看亚洲国产| 建设人人有责人人尽责人人享有的| 欧美亚洲 丝袜 人妻 在线| 一边摸一边做爽爽视频免费| 久9热在线精品视频| 妹子高潮喷水视频| 99国产精品一区二区蜜桃av | 国产亚洲欧美98| 亚洲国产精品合色在线| 捣出白浆h1v1| 中文字幕最新亚洲高清| 精品一品国产午夜福利视频| 极品人妻少妇av视频| 韩国精品一区二区三区| 国产精华一区二区三区| 一级a爱视频在线免费观看| 亚洲精品中文字幕一二三四区| 香蕉久久夜色| 成年人免费黄色播放视频| 精品第一国产精品| 80岁老熟妇乱子伦牲交| 视频区图区小说| 欧美激情久久久久久爽电影 | av福利片在线| 亚洲在线自拍视频| 亚洲一区二区三区不卡视频| avwww免费| 国产熟女午夜一区二区三区| 欧美不卡视频在线免费观看 | 久久国产精品影院| 国产精品av久久久久免费| www日本在线高清视频| 啪啪无遮挡十八禁网站| 欧美日韩黄片免| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| 自拍欧美九色日韩亚洲蝌蚪91| av在线播放免费不卡| 亚洲综合色网址| 国产又色又爽无遮挡免费看| 黑人欧美特级aaaaaa片| 一个人免费在线观看的高清视频| 久久久久国产精品人妻aⅴ院 | 老熟女久久久| 欧美另类亚洲清纯唯美| 日日夜夜操网爽| 亚洲人成电影免费在线| 亚洲欧美精品综合一区二区三区| 亚洲成人国产一区在线观看| 亚洲成国产人片在线观看| 精品一区二区三区视频在线观看免费 | 亚洲成人免费电影在线观看| 亚洲欧美一区二区三区黑人| xxx96com| avwww免费| 老司机福利观看| 国产日韩欧美亚洲二区| 亚洲色图 男人天堂 中文字幕| 国产视频一区二区在线看| 亚洲黑人精品在线| 免费在线观看视频国产中文字幕亚洲| 亚洲av成人av| 欧美久久黑人一区二区| 国产蜜桃级精品一区二区三区 | 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 亚洲精品一卡2卡三卡4卡5卡| 少妇裸体淫交视频免费看高清 | 热re99久久精品国产66热6| 亚洲人成电影观看| 天堂俺去俺来也www色官网| 12—13女人毛片做爰片一| 国产精品自产拍在线观看55亚洲 | 亚洲精品美女久久av网站| 别揉我奶头~嗯~啊~动态视频| 三上悠亚av全集在线观看| 两性夫妻黄色片| 午夜视频精品福利| 国产精品美女特级片免费视频播放器 | 王馨瑶露胸无遮挡在线观看| 在线永久观看黄色视频| 精品久久蜜臀av无| 国产一区有黄有色的免费视频| 一进一出抽搐动态| 王馨瑶露胸无遮挡在线观看| 亚洲中文av在线| 性色av乱码一区二区三区2| 欧美中文综合在线视频| 午夜免费成人在线视频| 久久这里只有精品19| 成人国产一区最新在线观看| 免费在线观看黄色视频的| 国产精品久久久久久精品古装| 18禁裸乳无遮挡免费网站照片 | 一边摸一边做爽爽视频免费| 1024视频免费在线观看| 97人妻天天添夜夜摸| 国产精品 国内视频| 亚洲三区欧美一区| 欧美日韩视频精品一区| 999久久久精品免费观看国产| 久久性视频一级片| 国产深夜福利视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 精品乱码久久久久久99久播| 国产亚洲精品一区二区www | 乱人伦中国视频| 一级片'在线观看视频| 亚洲成av片中文字幕在线观看| 好看av亚洲va欧美ⅴa在| 视频在线观看一区二区三区| 国产深夜福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 美女高潮到喷水免费观看| 不卡av一区二区三区| 久久九九热精品免费| 亚洲美女黄片视频| 黑人欧美特级aaaaaa片| 91av网站免费观看| 激情视频va一区二区三区| 女人精品久久久久毛片| 91麻豆精品激情在线观看国产 | 色在线成人网| 人人澡人人妻人| 91av网站免费观看| 精品国产亚洲在线| 国产成人欧美在线观看 | 好看av亚洲va欧美ⅴa在| 丝袜人妻中文字幕| 新久久久久国产一级毛片| 国产精品欧美亚洲77777| 成人免费观看视频高清| 97人妻天天添夜夜摸| 18禁黄网站禁片午夜丰满| 日本欧美视频一区| 国产精品电影一区二区三区 | 亚洲国产欧美日韩在线播放| 国产精品亚洲av一区麻豆| 亚洲成人国产一区在线观看| 精品久久久精品久久久| 嫁个100分男人电影在线观看| 男女午夜视频在线观看| 正在播放国产对白刺激| 国产精品一区二区在线观看99| 一本大道久久a久久精品| 成年人黄色毛片网站| 一级a爱片免费观看的视频| 免费观看精品视频网站| 最近最新中文字幕大全免费视频| 欧美日韩亚洲高清精品| 国产成人啪精品午夜网站| 色94色欧美一区二区| 又紧又爽又黄一区二区| 女人精品久久久久毛片| 亚洲成av片中文字幕在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧美色视频一区免费| 欧美日韩亚洲高清精品| www.精华液| 午夜福利在线免费观看网站| 亚洲熟女毛片儿| 少妇的丰满在线观看| 国产一区二区三区综合在线观看| 欧美国产精品va在线观看不卡| 青草久久国产| 国产成人系列免费观看| 亚洲色图综合在线观看| 亚洲精华国产精华精| 伊人久久大香线蕉亚洲五| 亚洲欧洲精品一区二区精品久久久| 露出奶头的视频| 在线观看66精品国产| 国产97色在线日韩免费| 十八禁网站免费在线| 手机成人av网站| 很黄的视频免费| 亚洲男人天堂网一区| 久久精品亚洲精品国产色婷小说| 欧美色视频一区免费| 久久精品aⅴ一区二区三区四区| 亚洲av片天天在线观看| 国产色视频综合| 国产在线观看jvid| 久久 成人 亚洲| 夜夜爽天天搞| 国产日韩欧美亚洲二区| 怎么达到女性高潮| 国产精品 国内视频| 免费看a级黄色片| 国产精品一区二区免费欧美| 男女下面插进去视频免费观看| 丝瓜视频免费看黄片| 亚洲情色 制服丝袜| 男女下面插进去视频免费观看| 国产深夜福利视频在线观看| 两个人看的免费小视频| 一夜夜www| 91字幕亚洲| 久久国产精品人妻蜜桃| 国产亚洲欧美在线一区二区| 亚洲三区欧美一区| 高清av免费在线| 亚洲性夜色夜夜综合| 欧美亚洲日本最大视频资源| 多毛熟女@视频| 免费观看人在逋| 欧美亚洲日本最大视频资源| 色尼玛亚洲综合影院| 香蕉丝袜av| 亚洲专区中文字幕在线| 成人影院久久| 午夜亚洲福利在线播放| 一进一出抽搐gif免费好疼 | 激情视频va一区二区三区| 国产男女内射视频|