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

    連續(xù)波腔衰蕩光譜技術(shù)中模式篩選的數(shù)值方法*

    2019-12-24 08:21:42王金舵余錦貊澤強3何建國3代守軍孟晶晶王曉東劉洋3
    物理學(xué)報 2019年24期
    關(guān)鍵詞:基模高階閾值

    王金舵 余錦 ? 貊澤強3) 何建國3) 代守軍 孟晶晶 王曉東 劉洋3)

    1) (中國科學(xué)院光電研究院, 北京 100094)

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

    3) (中國科學(xué)院計算光學(xué)成像技術(shù)重點實驗室, 北京 100094)

    連續(xù)波腔衰蕩光譜(CW-CRDS)測量中多模衰蕩的產(chǎn)生是嚴(yán)重影響痕量氣體測量靈敏度的重要因素.本文針對衰蕩腔內(nèi)無光闌或光闌濾模不徹底的CRDS裝置, 通過分析腔誤調(diào)時的能量耦合規(guī)律以及受關(guān)斷時間影響的衰蕩過程, 提出閾值選擇和擬合度判定兩種非光闌模式篩選方法, 利用數(shù)值方法達到抑制多模衰蕩及篩選基模衰蕩(優(yōu)衰蕩)的目的.首先對CW-CRDS實驗中平均采樣和單次采樣模式下出現(xiàn)的多種衰蕩類型進行了歸納分析, 發(fā)現(xiàn)可以通過單次采樣數(shù)據(jù)預(yù)測多次采樣的測量結(jié)果, 實驗結(jié)果與預(yù)期一致.解決了CRDS實驗“平均”和“擬合”的先后順序問題.在此基礎(chǔ)上, 利用優(yōu)衰蕩出現(xiàn)概率滿足二項分布模型的特性,建立了優(yōu)衰蕩出現(xiàn)頻率隨觸發(fā)閾值變化的概率模型, 用于選擇合適的觸發(fā)閾值.實驗表明提升觸發(fā)閾值可以有效地抑制多模衰蕩, 使測量靈敏度提升約一個數(shù)量級.隨著觸發(fā)閾值的提升, 通過優(yōu)衰蕩得到的Allan方差將趨于一個定值, 但是衰蕩過程獲取時間將逐漸延長.因此, 在CW-CRDS檢測中觸發(fā)閾值應(yīng)設(shè)置在保證全部衰蕩過程均為優(yōu)衰蕩的最小閾值處.之后, 采用擬合度判定法對實驗數(shù)據(jù)進行了篩選.最后給出了兩種方法的適用范圍, 擬合度判定法雖然簡單但局限性較大, 閾值選擇法可適用于腔誤調(diào)程度不嚴(yán)重的情況.

    1 引 言

    近年來, 腔衰蕩光譜技術(shù) (cavity ring-down spectroscopy, CRDS)在光物相互作用中等離子體分析[1,2]、溫室氣體測量[3,4]、污染氣體檢測[5]、分子光譜分析[6,7]以及反射率測量[8,9]等諸多方面得到應(yīng)用.CRDS是一種具有高靈敏度的直接吸收光譜技術(shù), 通過測量特定波長激光的衰蕩時間來實現(xiàn)弱吸收物質(zhì)的成分與濃度檢測[10].CRDS裝置的測量靈敏度與衰蕩時間的測量精度直接關(guān)聯(lián)[11], 衰蕩時間測量精度越高, 裝置測量靈敏度也越高.另一方面, 同樣的工作條件下, 激光激發(fā)衰蕩腔模式不同, 所對應(yīng)的衰蕩時間也存在明顯差異, 最終影響測量精度與準(zhǔn)確度.研究表明[12], 與單模衰蕩相比, 多模衰蕩激發(fā)下的CRDS測量靈敏度會大幅降低.

    自CRDS技術(shù)出現(xiàn)以來, 入射激光與衰蕩腔之間的模式匹配及衰蕩過程的篩選一直被認(rèn)為是一個重要的研究課題.在O’Keefe等[13]首次實現(xiàn)的腔衰蕩光譜測量中, 采用線寬較寬的脈沖激光器作為光源, 雖然簡化了實現(xiàn)激光-衰蕩腔頻率匹配的條件, 但是由于多個腔縱模的同時存在、脈沖間不同階數(shù)橫模的任意激發(fā)等因素而大幅降低了衰蕩時間的測量精度[14].Romanini等[15]將光源改為窄線寬的連續(xù)激光器, 通過壓電陶瓷(PZT)調(diào)節(jié)腔長實現(xiàn)了連續(xù)波腔衰蕩光譜技術(shù)(continuouswave cavity ring-down spectroscopy, CW-CRDS).這種結(jié)構(gòu)由于可以最大程度地避免多個腔縱模被同時激發(fā), 從而有效提高了CRDS的測量靈敏度及信噪比.與脈沖 CRDS 相同, CW-CRDS 也需要引入適當(dāng)?shù)臋M模匹配來抑制高階橫模的激發(fā).然而在實際應(yīng)用中, 由于加工及裝調(diào)誤差, 嚴(yán)格的橫模匹配難以實現(xiàn).Huang和Lehmann[16]對CW-CRDS中由于橫模匹配引入的噪聲進行了分析.研究表明, 基橫模與高階橫模同時激發(fā)會導(dǎo)致衰蕩曲線偏離理想的單指數(shù)衰減形式, 而更加符合雙指數(shù)模拍模型.其通過在衰蕩腔內(nèi)加入孔徑光闌有效抑制了高階橫模的產(chǎn)生, 使系統(tǒng)性能得到顯著提升.但是實驗中仍然記錄到了由于裝調(diào)偏差而產(chǎn)生的TEM01和TEM10模的衰蕩過程, 兩者衰蕩時間與基橫模(TEM00)不同, 從而影響了測量靈敏度.崔立紅等[17]對橫模匹配中的腔誤調(diào)進行了分析.腔誤調(diào)包括腔失調(diào) (misalignment)和腔失配 (mismatch), 其中腔失調(diào)是指入射光的光軸與衰蕩腔的光軸存在偏移或傾斜, 腔失配是指耦合進衰蕩腔的入射光腰斑與衰蕩腔基模腰斑之間存在尺寸或位置差異.該研究指出能量耦合效率主要受腔失配的影響, 而對衰蕩腔高階模激發(fā)來說, 腔失調(diào)的影響要大于腔失配.其提出可以采用法布里-珀羅干涉儀掃描法和多維象限儀探測法對特定激光器及衰蕩腔的腔失配及腔失調(diào)進行監(jiān)測, 通過PZT調(diào)節(jié)回路對衰蕩腔進行裝調(diào)校正, 減少多橫模衰蕩過程的發(fā)生, 從而提高CRDS裝置的測量靈敏度.但是這種方法的缺點在于需要使用額外的儀器(元器件)及復(fù)雜的反饋控制電路才能滿足較高精度裝調(diào)的需求.

    本文結(jié)合能量耦合及光腔響應(yīng)理論, 分析CW-CRDS實驗中存在的多(橫)模衰蕩現(xiàn)象, 提出了通過選擇合適的觸發(fā)閾值(閾值選擇)和分析衰蕩曲線擬合度(擬合度判定)實現(xiàn)CW-CRDS模式篩選的方法, 分析了兩種方法在痕量氣體檢測中的優(yōu)點和適用范圍.同時, 利用二項分布概率模型分析了在衰蕩曲線擬合前后進行數(shù)據(jù)平均的區(qū)別,并分別進行了實驗驗證, 回答了長期困擾衰蕩測量中“平均后擬合還是擬合后平均”的問題[18].

    2 理 論

    在CW-CRDS中, 入射光耦合進入衰蕩腔, 當(dāng)出射光強達到觸發(fā)閾值時, 關(guān)斷入射光, 衰蕩腔輸出形成一個衰蕩過程.理想情況下, 耦合過程中腔內(nèi)只有基橫模(以下簡稱為基模)被激發(fā), 探測器接收到出射光強隨時間的變化呈現(xiàn)單指數(shù)下降趨勢, 相應(yīng)地, 衰蕩時間呈正態(tài)分布, 計算其均值及標(biāo)準(zhǔn)差可以得到CRDS的測量靈敏度[11].但是, 實際上由于存在著腔失調(diào)與腔失配, 衰蕩腔內(nèi)除了基模外還會有高階橫模(以下簡稱為高階模)成分被激發(fā).由于空間分布不同導(dǎo)致腔內(nèi)損耗不同[16,19],高階模的衰蕩時間會遠(yuǎn)低于基模的衰蕩時間, 從而影響整體衰蕩時間的測量精度.實驗中可以利用基模與高階模的空間分布不同這一特點, 通過在衰蕩腔內(nèi)加入孔徑光闌, 來抑制高階模的產(chǎn)生.但是,腔內(nèi)孔徑光闌的引入同時也會造成部分基模衰蕩能量的損失, 進而導(dǎo)致信噪比的降低.因此, 有必要能夠找到某些非光闌方法來抑制甚至消除高階模對衰蕩時間測量精度帶來的影響.

    這里引入失調(diào)參量ξ和失配參量ο, 假設(shè)入射光的腰斑位置在z軸坐標(biāo)0點, 且在兩個方向初始相位均為0, 則有[20]

    式中以下標(biāo)L表示入射激光的參量, 下標(biāo)C表示衰蕩腔的參量;ε為入射光光軸與衰蕩腔光軸的偏移量;γ為入射光光軸與衰蕩腔光軸的傾角;w為光斑半徑;q為高斯光束q參數(shù).

    為簡化分析, 在分析腔失配時, 假設(shè)腰斑位置良好匹配, 僅光斑尺寸存在差異, 則(1)式可以改寫為

    其中θ為入射光的半發(fā)散角.

    在實際裝調(diào)中, 最常見的腔誤調(diào)是經(jīng)橫模匹配后, 入射激光腰斑與衰蕩腔基模腰斑之間存在的尺寸差異.通過引入裝調(diào)指示光, 可以使入射光光軸與衰蕩腔光軸基本交疊, 使得軸線偏移量ε近似為0.假設(shè)腔失調(diào)的傾角為入射光半發(fā)散角的1/10, 入射光束腰斑半徑wL與衰蕩腔基橫模腰斑半徑wC之差為Δw, 這里用Δw與wC的比值來表征腔誤調(diào)程度.在僅存在腔失配的情況下, 耦合進高階模的光功率Pmn占激光耦合進入衰蕩腔的總功率P0的比例[20]表示為

    同時存在腔失調(diào)和腔失配的情況下, 耦合進高階模的光功率Pmn占激光耦合進衰蕩腔的總功率P0的比例[20]表示為

    將(2)式代入(3)和(4)式, 可以分別得到兩種情況下高階模耦合功率占比隨腔誤調(diào)度(Δw/wC)的變化, 如圖1所示.

    圖1 腔誤調(diào)時高階模功率耦合占比情況Fig.1.Proportion of higher-order cavity mode excitement in a misadjusted CRDS system.

    從圖1可以看到, 當(dāng)衰蕩腔出現(xiàn)誤調(diào)時, 因部分激光能量的耦合使高階模受到激發(fā), 引發(fā)高階模衰蕩過程.在腔誤調(diào)程度較輕時, 耦合進入高階模的能量較少, 其輸出光強較弱, 出現(xiàn)高階模衰蕩的概率較低.隨著模式失配的進一步惡化, 耦合到高階模的能量迅速增加.在衰蕩曲線采集過程中, 這將導(dǎo)致兩方面的后果:一是從衰蕩腔輸出的高階模光束光強增加; 二是同一觸發(fā)閾值下, 采集到的衰蕩過程屬于高階模衰蕩的概率增加.由此可知, 在腔誤調(diào)程度較輕時, 耦合進衰蕩腔的激光能量大部分注入腔基模中, 高階模具有較少的能量, 此時探測器接收到的基模光強會明顯高于高階模光強.在這種情況下, 在獲取衰蕩過程中可以通過選擇合適的觸發(fā)閾值, 對不同模式的衰蕩過程進行篩選, 只保留基模衰蕩過程.

    理想情況下, 激光關(guān)斷是在瞬時完成, 衰蕩曲線應(yīng)呈單指數(shù)下降.但實際測量中, 受有限的關(guān)斷時間影響, 衰蕩曲線會偏離單指數(shù)趨勢, 偏離程度與光源關(guān)斷時間及理想衰蕩時間相關(guān).假設(shè)在tp時刻, 光源開啟充光; 在td時刻, 出射光強到達觸發(fā)閾值, 光源關(guān)斷, 則入射光場振幅Ei的時域變化可以寫為[21]

    其中E為光源光場振幅,η為關(guān)斷效率,ts為關(guān)斷時間,u(t)為階躍函數(shù).關(guān)斷后, 出射光場振幅Eo的時域表達式為[22]

    其中τ為衰蕩時間.不失一般性, 假設(shè)入射光場電場強度幅值為在零時刻發(fā)生關(guān)斷, 關(guān)斷時間為 200 ns, 如圖2 中插圖所示.基模衰蕩時間設(shè)為 31.5 μs, 由于損耗較大, 高階模的衰蕩時間較短, 設(shè)為 19 μs.探測過程中的高斯噪聲功率為入射光光強的0.025倍, 則探測器可能接收到的出射光強信號Io隨時間的變化如圖2所示, 其中基模衰蕩記為優(yōu)衰蕩(good decay), 高階模衰蕩記為劣衰蕩(bad decay).衰蕩曲線的單指數(shù)性用調(diào)整擬合度來表示, 相較于傳統(tǒng)擬合度, 調(diào)整擬合度剔除了變量個數(shù)對擬合度的影響, 其確定系數(shù)的表達式為

    其中SSres為樣本殘差平方和, 自由度為dfr;SStot為樣本離差的平方和, 自由度為dft.從圖2可以看到, 基模衰蕩曲線的調(diào)整擬合度要高于高階模衰蕩曲線的調(diào)整擬合度.因此可以認(rèn)為, 在測得衰蕩曲線后, 通過比較調(diào)整擬合度, 能夠?qū)Σ煌J降乃ナ庍^程進行甄別, 從而篩選出基模衰蕩過程.

    圖2 受有限的關(guān)斷時間影響的衰蕩曲線仿真Fig.2.Simulation of ring-down curves affected by the finite shutdown time.

    3 實 驗

    用于研究CW-CRDS非光闌模式篩選的實驗裝置為實驗室搭建的CW-CRDS裝置, 未安裝腔內(nèi)光闌, 如圖3所示.光源采用光纖耦合輸出的分布反饋式二極管激光器(DFB-LD, NEL), 中心波長為1654 nm.激光器控制器 (LDC-3724C, Lightwave)可以通過控制電流或溫度, 對DFB-LD輸出波長進行調(diào)節(jié), 電流設(shè)定為 130 mA, 溫度設(shè)定為 20 ℃.為保證入射激光耦合進衰蕩腔, 需要實現(xiàn)穩(wěn)定的頻率匹配, 通過信號發(fā)生器 (33600A, Aglient)在激光器控制器上加載一個三角波電流調(diào)制信號, 調(diào)制頻率為 10 Hz, 幅值為 26 mV.該幅值對激光頻率的調(diào)制略大于一個自由光譜范圍.激光光束經(jīng)準(zhǔn)直器、模式匹配透鏡組, 耦合進入衰蕩腔.入射光路中加入光隔離器(Thorlabs)來抑制光反饋.衰蕩腔的腔鏡采用兩片相同參數(shù)的平凹高反鏡, 曲率半徑為 500 mm, 在 1600—1700 nm 范圍內(nèi), 標(biāo)稱反射率優(yōu)于 0.9999, 腔鏡間隔為 340 mm.出射光經(jīng)聚焦透鏡后, 是由光電探測器(PDA10CS-ES,Thorlabs)接收、示波器 (MSO 4104C, Tektronix)進行采集記錄.當(dāng)出射光強達到觸發(fā)閾值時, 通過示波器向激光控制器發(fā)出觸發(fā)信號, 實現(xiàn)關(guān)斷.利用示波器對衰蕩曲線進行記錄, 用于后續(xù)數(shù)據(jù)處理.實驗氣體為標(biāo)準(zhǔn)氮氣(N2), 通過質(zhì)量流量控制器, 以 0.5 L/min 的速度輸入衰蕩腔.為減小水汽及空氣顆粒對實驗結(jié)果的影響, 在衰蕩腔輸入輸出氣路上, 均裝有干燥劑和顆粒過濾器.

    4 實驗結(jié)果與分析

    4.1 衰蕩過程類型分析

    圖3 CW-CRDS 實驗裝置圖Fig.3.Schematic of CW-CRDS experimental instrument.

    圖4 四次平均模式下的 (a)衰蕩曲線、(b), (c)衰蕩時間及 (d)擬合度分布Fig.4.(a) Typical decays, (b), (c) distributions of ring-down time, and (d) decay curve fitness in four times averaged mode sampling.

    利用數(shù)字存儲示波器對衰蕩過程進行記錄, 采樣模式為平均采樣模式, 平均次數(shù)為4次, 總共采集10000條衰蕩曲線, 實驗結(jié)果如圖4所示.在10000組實驗數(shù)據(jù)中, 根據(jù)階梯分布的衰蕩時間可以劃分為三種類型的衰蕩過程, 其典型的衰蕩曲線如圖4(a)所示.三類衰蕩過程所占比例不同, 其中絕大部分是衰蕩時間均值約為32 μs的衰蕩過程,比例為92.84%, 為基模衰蕩(優(yōu)衰蕩).另外兩類衰蕩過程的衰蕩時間較短, 推測其為高階模衰蕩, 稱之為劣衰蕩1和劣衰蕩2, 均值分別為28.6 μs和25.8 μs; 兩者所占比例較低, 分別為 6.86% 和 0.3%,如圖4(b)所示.通過衰蕩時間分布, 篩選出優(yōu)衰蕩過程, 其衰蕩時間呈正態(tài)分布, 如圖4(c) 所示.如果衰蕩過程只存在優(yōu)衰蕩, 可以得到較小的衰蕩時間標(biāo)準(zhǔn)差, 對應(yīng)的測量靈敏度較高.但是, 由于存在衰蕩時間相差較大的兩類劣衰蕩, 導(dǎo)致標(biāo)準(zhǔn)差急劇增大, 測量靈敏度大幅下降.對衰蕩曲線的擬合度進行計算, 發(fā)現(xiàn)三種類型衰蕩過程的衰蕩曲線擬合度存在嚴(yán)重交疊, 如圖4(d)所示.在這種情況下, 很難提出有效判據(jù)來實現(xiàn)模式篩選.另外, 假設(shè)這三類衰蕩過程分別對應(yīng)衰蕩腔內(nèi)同時存在的基模衰蕩、以及其他兩種不同階數(shù)的高階模衰蕩,那么受腔損耗影響, 這三種類型的衰蕩時間之間的差值應(yīng)該逐漸增大[19], 這與三類衰蕩時間均值分布大致成等差數(shù)列的實驗結(jié)果相矛盾, 即原假設(shè)不成立.

    為證實實驗中出現(xiàn)了多模衰蕩, 利用掃描頻譜法測量了衰蕩腔出射激光模式的頻率分布, 如圖5所示, 其中 A, B, C 為相鄰的三個強激發(fā)腔模, a, b是兩個弱激發(fā)腔模, 介于 A, B, C 腔模之間.測得A與B, B與C的模式頻率間隔分別為445.9 MHz和 439.2 MHz, A 與 a, B 與 b 的模式頻率間隔分別為 174.9 MHz和 171.4 MHz.與利用間距L的兩面曲率半徑為Rc的腔鏡構(gòu)建的對稱圓形孔徑穩(wěn)定球面腔諧振頻率公式為

    其中c為光速;m,n為橫模階次.計算得到的基橫模頻率ν00= 441.2 MHz 和高階模v10= 174.9 MHz基本一致.由此推測, 衰蕩腔內(nèi)主要有兩種模式被激發(fā), 即 TEM00模和 TEM10模.

    圖5 衰蕩腔內(nèi)存在多模激發(fā)Fig.5.Multimode excitation in ring-down cavity.

    采用上述CW-CRDS裝置再次進行實驗, 將示波器的采樣模式變?yōu)闊o平均過程的單次采樣模式, 采集10000條衰蕩曲線.在10000組實驗數(shù)據(jù)中, 只存在兩類衰蕩過程, 衰蕩時間均值為 31.5 μs的優(yōu)衰蕩和均值為19.1 μs的劣衰蕩, 前者占總數(shù)的98.97%, 其典型的衰蕩曲線如圖6所示.計算兩條衰蕩曲線的調(diào)整擬合優(yōu)度, 與理論仿真(圖2)的結(jié)果相近.

    圖6 單次采樣的典型衰蕩曲線Fig.6.Typical decays in single sampling.

    分析上述兩種不同采樣模式下的實驗結(jié)果, 認(rèn)為受腔誤調(diào)影響, 實驗所用的CW-CRDS裝置中有且僅有兩種類型的衰蕩過程, 這里分別被稱為“優(yōu)衰蕩”和“劣衰蕩”, 這兩個事件為互斥事件、且其事件出現(xiàn)概率之和為1.由10000組單次采樣實驗結(jié)果表明, 前者出現(xiàn)頻率(約為0.99)遠(yuǎn)高于后者(約為0.01), 這表明耦合進衰蕩腔能量絕大部分在基模上, 即CW-CRDS裝置的腔誤調(diào)程度較輕,可通過選擇合適的觸發(fā)閾值來實現(xiàn)模式篩選.

    在四次平均模式下, 這兩種類型的衰蕩過程可以有五種組合方式:四個優(yōu)衰蕩、三個優(yōu)衰蕩一個劣衰蕩、兩個優(yōu)衰蕩兩個劣衰蕩、一個優(yōu)衰蕩三個劣衰蕩、以及四個劣衰蕩.由單次采樣下優(yōu)劣衰蕩出現(xiàn)頻率推得上述五種組合出現(xiàn)概率分別為0.9606, 0.0388, 5.9 × 10—4, 4.0 × 10—6和 1.0 ×10—8.為方便表述, 這里定義平均模式采樣中的“優(yōu)衰蕩”是全部由單次采樣中的優(yōu)衰蕩組成, 其他組合方式均為“劣衰蕩”, 并根據(jù)所含單次采樣中的劣衰蕩數(shù)目進行排序.按上述概率對10000組四次平均模式下采集到的衰蕩過程進行估計, 則優(yōu)衰蕩、劣衰蕩1、劣衰蕩2、劣衰蕩3和劣衰蕩4的頻數(shù)約為 9606, 388, 6, 0 和 0.由此可推測在四次平均模式下采集10000組衰蕩過程, 看到的衰蕩過程類型可能為三類, 優(yōu)衰蕩、劣衰蕩1和劣衰蕩2.該推測結(jié)果與10000組四次平均模式下的實驗結(jié)果(觀測到三種類型衰蕩個數(shù)分別為9284, 686和30)接近.據(jù)此推測, CW-CRDS 實驗中, 以單次采樣模式進行衰蕩過程采集, 在一次采集中結(jié)果為優(yōu)衰蕩屬于單次伯努利實驗;N次衰蕩過程采集是相互獨立事件, 優(yōu)衰蕩出現(xiàn)次數(shù)滿足二項分布.在平均模式下進行衰蕩過程采集, 假設(shè)平均次數(shù)為N(N≥ 2),則可能存在的劣衰蕩的種類為N, 在已知優(yōu)衰蕩出現(xiàn)概率的基礎(chǔ)上, 利用二項式展開可以得到各類型衰蕩過程出現(xiàn)的概率.已知兩類基本衰蕩過程的衰蕩時間, 可以推測出不同采樣模式下不同類型衰蕩時間的數(shù)學(xué)期望.在同一采樣模式下, 多個類型衰蕩時間期望值之間呈等差數(shù)列, 這就解釋了四次平均模式采樣實驗中三類衰蕩時間大致呈等差數(shù)列的現(xiàn)象.另外, 只要在單次采樣模式下抑制或消除了劣衰蕩, 使其出現(xiàn)概率趨近于0, 則之后不論采用哪種模式進行采樣, 都能獲得較好的測量靈敏度.所以, 對于衰蕩過程均為優(yōu)衰蕩的情況, “平均后擬合”還是“擬合后平均”除因處理速度帶來的差異外, 無明顯測量精度差別.但當(dāng)存在劣衰蕩時,“平均后擬合”將造成劣衰蕩與優(yōu)衰蕩混疊, 增大模式篩選難度, 且平均次數(shù)越多, 篩選難度越大, 提取的衰蕩時間精度越差.這種情況下, 應(yīng)該先以單次采樣模式對衰蕩過程進行采集, 對衰蕩曲線進行擬合, 提取衰蕩時間, 經(jīng)篩選后再平均計算.

    4.2 閾值選擇模式篩選法

    通過理論分析可知, 提高觸發(fā)閾值可以對不同模式的衰蕩過程進行篩選.本文通過建立優(yōu)衰蕩出現(xiàn)頻率隨觸發(fā)閾值變化的概率模型, 對合適觸發(fā)閾值的選擇進行分析.根據(jù)單次采樣與平均采樣間關(guān)系的論述可以知道, 利用不同閾值單次采樣模式獲得的優(yōu)衰蕩出現(xiàn)頻率可以建立N次(N≥ 2)平均模式優(yōu)衰蕩出現(xiàn)頻率隨觸發(fā)閾值變化概率模型, 不失一般性, 本文以N= 4 為例.

    通過上文分析, 在單次采樣模式下優(yōu)衰蕩的出現(xiàn)次數(shù)fg滿足二項分布, 即fg~B(n,p).實驗中,為建立可靠的概率模型, 通過有限次測量過程中優(yōu)衰蕩的出現(xiàn)頻率去推斷總體概率p, 需要對實驗所需的衰蕩事件個數(shù), 即抽樣容量n進行推斷.二項分布概率數(shù)學(xué)模型表示為

    其中δ為上述推斷的精度, 1—α為置信度.

    當(dāng)n較大時, 利用棣莫弗-拉普拉斯中心極限定理, 將二項分布用其極限分布正態(tài)分布表示:

    由(8)和(9)式可以得到

    通過正態(tài)分布概率值表, 可以得到概率臨界值x0.則有

    當(dāng)p= 0.5 時, (11)式等號右邊有最大值, 則滿足推斷所需精度和置信度的抽樣容量為該最大值取整后加一.

    假設(shè)CW-CRDS裝置以驅(qū)動電流三角波調(diào)制方式實現(xiàn)波長調(diào)節(jié), 最終獲得1 Hz甲烷濃度測量速率.假定三角波頻率10 Hz, 相應(yīng)入射激光光譜掃描范圍略大于一個自由光譜范圍[23], 理論上1 s內(nèi)最多可有效觸發(fā)40個衰蕩事件.概率模型所需精度δ= 1/(2 × 40), 約為 0.013.取α= 0.1,則概率臨界值為1.65, 計算得到所需抽樣容量為4028.為便于分析表述, 本文選取的抽樣容量為4000, 也即需要獲取4000個衰蕩事件才能分析優(yōu)衰蕩出現(xiàn)頻率.

    分別選取 10.7, 14.7, 18.7 和 22.7 mV 作為觸發(fā)閾值, 以單次采樣模式, 分別采集4000條衰蕩曲線, 得到優(yōu)衰蕩出現(xiàn)頻率分別為0.9903, 0.9953,0.9990和1.利用二項分布分別計算得出在四次平均采樣模式下, 上述觸發(fā)閾值出現(xiàn)優(yōu)衰蕩出現(xiàn)頻率如圖7所示.

    以指數(shù)函數(shù)進行擬合, 則優(yōu)衰蕩出現(xiàn)頻率擬合曲線為

    該方程給出了在實驗條件下, 由單次采樣數(shù)據(jù)得到的四次平均采樣模式下優(yōu)衰蕩出現(xiàn)頻率隨觸發(fā)閾值變化的概率模型.為證明模型的正確性, 在四次平均采樣模式下, 分別以 10.7, 12.7, 16.7 和 20.7 mV為觸發(fā)閾值進行實驗, 在每個觸發(fā)閾值下分別進行4組實驗, 每組采集1000條衰蕩曲線, 得到優(yōu)衰蕩出現(xiàn)頻率.從圖7可以看到, 實驗測得結(jié)果與概率模型預(yù)測結(jié)果基本相符.

    從圖7還可以看出, 隨著觸發(fā)閾值提高, 優(yōu)衰蕩出現(xiàn)的概率趨近于1, 且出現(xiàn)頻率的波動區(qū)間逐漸收斂.利用該概率模型, 可以大致得到不同觸發(fā)閾值下, 優(yōu)衰蕩出現(xiàn)的概率, 為閾值選擇模式篩選提供參考.

    圖7 優(yōu)衰蕩出現(xiàn)頻率隨觸發(fā)閾值的變化Fig.7.Variance of the frequency of good decays with trigger thresholds.

    通過衰蕩時間長短, 對單次采樣模式下四種不同觸發(fā)閾值 (10.7, 14.7, 18.7 和 22.7 mV)得到的衰蕩過程進行篩選, 則篩選前后各自Allan方差曲線如圖8所示.由圖8可以看到, 篩選后的Allan方差比原始數(shù)據(jù)大約小一個數(shù)量級.高觸發(fā)閾值(22.7 mV)下衰蕩過程直接得到的Allan方差與其他觸發(fā)閾值下經(jīng)篩選后的數(shù)值相近.

    對單次/四次平均采樣模式下, 不同觸發(fā)閾值(10.7, 12.7, 14.7, 16.7, 18.7, 20.7, 22.7 mV)實驗中由40個衰蕩過程計算得到的Allan方差值進行對比, 結(jié)果如圖9所示.在單次采樣模式下, 示波器采集的一條衰蕩曲線代表一個衰蕩過程, Allan方差圖中橫坐標(biāo)衰蕩曲線個數(shù)為40; 在四次平均模式下, 示波器采集的一條衰蕩曲線是由四個衰蕩過程平均后得到的, 為保證衰蕩過程個數(shù)為40,Allan方差圖中橫坐標(biāo)衰蕩曲線個數(shù)為10.從圖9可知, 隨著觸發(fā)閾值的提高, 40個衰蕩過程的Allan方差值線性減小.當(dāng)進行數(shù)據(jù)篩選后,Allan方差值隨著觸發(fā)閾值逐漸趨近于一個定值,對于本實驗采用的CW-CRDS裝置, 該值為3 ×10—10cm—1.

    圖8 不同閾值單次采樣下的 Allan 方差圖 (a) 10.7 mV; (b) 14.7 mV; (c) 18.7 mV; (d) 22.7 mVFig.8.Allan deviations in single sampling with different trigger thresholds:(a) 10.7 mV; (b) 14.7 mV; (c) 18.7 mV; (d) 22.7 mV.

    圖9 不同閾值下 40 個衰蕩過程的 Allan 方差Fig.9.Allan deviations of 40 decays with different trigger thresholds.

    實驗中發(fā)現(xiàn), 衰蕩過程獲取速率與觸發(fā)閾值密切相關(guān), 觸發(fā)閾值越高, 獲取速率越慢.在不同觸發(fā)閾值下, 取重復(fù)測量10次、每次采樣128個衰蕩過程所需時間的平均值作為采樣耗時, 可以計算得到該觸發(fā)閾值下噪聲等效吸收系數(shù)(noise equivalentabsorption coefficient, NEA)以及 1 Hz測量速率下甲烷測量靈敏度, 如圖10所示.從圖10可以看到, 隨著觸發(fā)閾值的提升, CW-CRDS裝置的測量靈敏度迅速提升, 其中觸發(fā)閾值為20.7 mV和22.7 mV時, 測量靈敏度大幅提升.由于觸發(fā)閾值為20.7 mV時獲取有效衰蕩記錄的過程更快,則1 Hz測量速率下相較于閾值22.7 mV時的測量靈敏度更高.所以, 對于實驗用CW-CRDS裝置, 通過閾值選擇法進行模式篩選時, 可以將觸發(fā)閾值設(shè)定為20.7 mV.通過基于衰蕩時間的數(shù)據(jù)篩選實現(xiàn)劣衰蕩剔除后, 可以看到不同觸發(fā)閾值下,該裝置對于甲烷的測量靈敏度大致相同, 當(dāng)觸發(fā)閾值為14.7 mV時, 相同檢測速率下測量靈敏度數(shù)值最小.

    圖10 不同閾值的噪聲等效吸收系數(shù)及測量靈敏度Fig.10.Noise equivalent absorption coefficients and sensitivities under different trigger thresholds.

    4.3 擬合度判定模式篩選法

    除閾值選擇外, 也可以通過擬合度判定實現(xiàn)模式篩選.相較于平均模式, 單次采樣模式下兩類衰蕩過程調(diào)整擬合度分布的交疊較少, 以觸發(fā)閾值為10.7 mV的單次采樣為例, 其衰蕩時間及衰蕩曲線調(diào)整擬合優(yōu)度如圖11所示.盡管交疊較少,但兩類衰蕩過程擬合度之間交疊現(xiàn)象仍然存在.如果以衰蕩曲線的調(diào)整擬合度為判據(jù)進行模式篩選,為保證劣衰蕩全部被剔除, 需要以劣衰蕩曲線擬合度的最大值為標(biāo)準(zhǔn)進行選擇, 這將導(dǎo)致在總共4000個衰蕩過程中, 除39個劣衰蕩外, 還有140個優(yōu)衰蕩被誤剔除.這將影響CW-CRDS裝置對優(yōu)衰蕩的獲取效率, 降低裝置測量靈敏度.

    圖11 單次采樣模式下的衰蕩時間(a)及曲線調(diào)整擬合度(b)分布Fig.11.Distributions of ring-down times (a) and curve fitness (b) in single sampling.

    4.4 模式篩選數(shù)值方法的適用范圍

    利用上述數(shù)值方法對CW-CRDS中模式篩選,雖然無需引入影響能量傳遞的光闌結(jié)構(gòu), 但是也存在相應(yīng)的適用范圍.理論分析可知, 由于實際測量中存在關(guān)斷時間、電路響應(yīng)延遲等因素影響, 衰蕩曲線不再是理想的單指數(shù)函數(shù).衰蕩時間越短, 對衰蕩曲線的單指數(shù)函數(shù)擬合優(yōu)度的影響越大.所以從本質(zhì)上看, 利用擬合度判別和在低觸發(fā)閾值下的數(shù)據(jù)篩選, 都是利用劣衰蕩的衰蕩時間遠(yuǎn)小于優(yōu)衰蕩這一特性.但事實上, 實際應(yīng)用中衰蕩時間可能會因為腔內(nèi)吸收損耗的大幅增加而急劇縮短.以本文搭建的CW-CRDS裝置為例, 空腔衰蕩時間為31.5 μs, 由 CRDS 痕氣檢測的基本公式[3]可知, 對于甲烷 (其吸收截面為 1.5 × 10—20cm2/molecule)氣體濃度測量, 利用擬合度判別或數(shù)據(jù)篩選來提高測量靈敏度的方法, 僅適用于氣體濃度不大于2 ppmv 的場合.當(dāng)甲烷濃度大于 2 ppmv 時, 衰蕩時間將小于19 μs, 與空腔狀態(tài)下劣衰蕩的衰蕩時間相近, 該方法將不再適用.而對于硫化氫(H2S)氣體濃度測量, 由于其吸收截面為2.7 ×10—24cm2/molecule), 則上述最大量程可達到 1 ×104ppmv.可以看到, 采用這種方法進行模式篩選,在一定程度上限制了CW-CRDS裝置的最大量程.而通過提高觸發(fā)閾值進行模式篩選則不會受到待測氣體濃度影響, 其本質(zhì)更依賴于CW-CRDS裝置自身性能, 即正常情況下, 在衰蕩腔中耦合進基模的能量遠(yuǎn)大于耦合進高階模的能量.通過(4)式計算可知, 在不存在腔失調(diào)的情況下, 入射光束腰斑半徑與模式匹配目標(biāo)(衰蕩腔基模腰斑半徑)的偏差小于40%時, 或經(jīng)良好匹配(腰斑半徑偏差為5%), 腔失調(diào)光軸傾角小于4倍入射光束發(fā)散角時, 可通過提高閾值進行模式篩選.

    5 結(jié) 論

    本文通過分析連續(xù)波激光衰蕩過程存在腔誤調(diào)情況時不同衰蕩模式之間的能量耦合、基于光腔響應(yīng)理論的衰蕩過程與激光關(guān)斷時間之間的關(guān)系,分別提出了利用提高觸發(fā)閾值與利用衰蕩曲線調(diào)整擬合優(yōu)度來進行模式篩選的方法, 并且在實驗中對這兩種模式篩選方法進行了驗證.在頻譜掃描以及模式分析的基礎(chǔ)上, 解釋了單次采樣與平均采樣之間的關(guān)系, 提出了優(yōu)衰蕩出現(xiàn)頻率隨觸發(fā)閾值變化的概率模型, 為閾值選擇法進行模式篩選提供了參考.實驗結(jié)果表明, 篩選后計算得到的Allan方差比直接使用原始數(shù)據(jù)大約小一個數(shù)量級.也就是說, 通過提升閾值或利用衰蕩時間進行的模式篩選能有效地提高衰蕩時間的測量精度, 從而提高CRDS的測量靈敏度.隨著觸發(fā)閾值的提高, 由優(yōu)衰蕩得到的Allan方差將趨近于一個定值, 但是衰蕩過程的獲取速率將逐漸降低; 二者平衡可以實現(xiàn)快速高精度采樣, 在滿足模式篩選的最小觸發(fā)閾值下進行腔衰蕩實驗.擬合度判定法雖然方法簡單,能起到模式篩選的作用, 但效率較低, 并且會限制痕量氣體檢測中的最大量程.閾值選擇法不受腔內(nèi)吸收損耗變化影響, 可以在采集衰蕩過程階段實現(xiàn)模式篩選, 但不適用于腔誤調(diào)較大的情況.

    猜你喜歡
    基模高階閾值
    有限圖上高階Yamabe型方程的非平凡解
    高階各向異性Cahn-Hilliard-Navier-Stokes系統(tǒng)的弱解
    滾動軸承壽命高階計算與應(yīng)用
    哈爾濱軸承(2020年1期)2020-11-03 09:16:02
    小波閾值去噪在深小孔鉆削聲發(fā)射信號處理中的應(yīng)用
    基于自適應(yīng)閾值和連通域的隧道裂縫提取
    比值遙感蝕變信息提取及閾值確定(插圖)
    河北遙感(2017年2期)2017-08-07 14:49:00
    從基模理論談新媒體環(huán)境下網(wǎng)民媒介素養(yǎng)的提高
    采寫編(2017年2期)2017-06-29 11:28:36
    “基模導(dǎo)向”在初中數(shù)學(xué)教學(xué)中的應(yīng)用
    室內(nèi)表面平均氡析出率閾值探討
    革新實驗室:一種新的工作場所學(xué)習(xí)方法的基模
    亚洲欧美精品综合一区二区三区| 色综合欧美亚洲国产小说| 免费高清视频大片| 一区福利在线观看| 欧美日韩精品网址| 亚洲国产精品久久男人天堂| 99riav亚洲国产免费| 免费高清视频大片| 国产三级在线视频| 又黄又粗又硬又大视频| 欧美中文综合在线视频| 免费看十八禁软件| 天堂√8在线中文| ponron亚洲| 亚洲午夜理论影院| 在线视频色国产色| 丁香六月欧美| 精品午夜福利视频在线观看一区| 亚洲一码二码三码区别大吗| 91麻豆精品激情在线观看国产| 高清在线国产一区| 国产亚洲av嫩草精品影院| 90打野战视频偷拍视频| 精品人妻1区二区| 亚洲一区中文字幕在线| 黑人欧美特级aaaaaa片| 久久久精品国产亚洲av高清涩受| 欧美日韩乱码在线| 国产欧美日韩一区二区三| 色哟哟哟哟哟哟| 国产亚洲精品一区二区www| 亚洲国产中文字幕在线视频| 精品欧美国产一区二区三| 黑丝袜美女国产一区| 国产又色又爽无遮挡免费看| 精品一区二区三区视频在线观看免费| 久久青草综合色| 性欧美人与动物交配| 高清毛片免费观看视频网站| 啦啦啦 在线观看视频| 精品久久久久久,| 国产精品亚洲一级av第二区| 一级毛片精品| 日本撒尿小便嘘嘘汇集6| 人成视频在线观看免费观看| 欧美激情 高清一区二区三区| 嫩草影视91久久| 97人妻天天添夜夜摸| 可以在线观看毛片的网站| 国产亚洲欧美精品永久| 欧美乱妇无乱码| 欧美日韩亚洲综合一区二区三区_| 国产欧美日韩一区二区三| 亚洲精品在线观看二区| 侵犯人妻中文字幕一二三四区| 亚洲色图综合在线观看| 精品午夜福利视频在线观看一区| 精品欧美国产一区二区三| 亚洲av熟女| 成人免费观看视频高清| 国产主播在线观看一区二区| 97碰自拍视频| 国产亚洲精品久久久久久毛片| 久久精品国产综合久久久| 国产人伦9x9x在线观看| 极品人妻少妇av视频| 亚洲国产欧美一区二区综合| 琪琪午夜伦伦电影理论片6080| 成年女人毛片免费观看观看9| a在线观看视频网站| 91成年电影在线观看| 久久婷婷人人爽人人干人人爱 | 亚洲av成人不卡在线观看播放网| 亚洲一区中文字幕在线| 黄色毛片三级朝国网站| 国产精品久久电影中文字幕| 国产精品美女特级片免费视频播放器 | 久久国产亚洲av麻豆专区| 熟女少妇亚洲综合色aaa.| 午夜老司机福利片| 一本久久中文字幕| 三级毛片av免费| 久久精品成人免费网站| 久久中文字幕人妻熟女| 久久久久久大精品| 久久天躁狠狠躁夜夜2o2o| aaaaa片日本免费| 午夜免费鲁丝| svipshipincom国产片| 成人国语在线视频| 精品熟女少妇八av免费久了| 免费人成视频x8x8入口观看| 国产亚洲欧美在线一区二区| 国产欧美日韩一区二区三| 亚洲成人精品中文字幕电影| 少妇粗大呻吟视频| 亚洲欧美激情在线| 日韩精品中文字幕看吧| 黄色毛片三级朝国网站| 久久人妻av系列| 精品午夜福利视频在线观看一区| 美女 人体艺术 gogo| 一区二区三区精品91| 久久精品国产亚洲av高清一级| 久久久久久久午夜电影| 国产精品二区激情视频| 久久狼人影院| 国产欧美日韩精品亚洲av| 久久久精品欧美日韩精品| 日韩高清综合在线| 国产欧美日韩综合在线一区二区| 亚洲avbb在线观看| 99久久99久久久精品蜜桃| 日韩有码中文字幕| 久久精品国产亚洲av高清一级| 亚洲人成电影观看| 久久久久久亚洲精品国产蜜桃av| 国产熟女午夜一区二区三区| 这个男人来自地球电影免费观看| 免费女性裸体啪啪无遮挡网站| 黄色视频不卡| 久久中文看片网| 极品教师在线免费播放| 久久久久九九精品影院| 国产av又大| 女人精品久久久久毛片| 久久热在线av| 免费在线观看视频国产中文字幕亚洲| 国产精品亚洲美女久久久| 一级毛片高清免费大全| 淫妇啪啪啪对白视频| 曰老女人黄片| 亚洲欧美激情在线| 一级a爱片免费观看的视频| 黄色女人牲交| 90打野战视频偷拍视频| 久久久国产精品麻豆| 久久精品影院6| 国产精品一区二区三区四区久久 | 淫妇啪啪啪对白视频| 黄色丝袜av网址大全| 99久久99久久久精品蜜桃| 精品国产超薄肉色丝袜足j| 亚洲色图av天堂| 欧美不卡视频在线免费观看 | 久久久水蜜桃国产精品网| 国产成年人精品一区二区| 国产精品,欧美在线| 高清黄色对白视频在线免费看| 午夜老司机福利片| 99在线视频只有这里精品首页| 一区二区三区精品91| www.www免费av| 欧美成人一区二区免费高清观看 | 欧美成狂野欧美在线观看| 日本a在线网址| 日本免费a在线| 日本免费a在线| 精品免费久久久久久久清纯| 久久中文字幕人妻熟女| 久久人妻福利社区极品人妻图片| 国产成人精品在线电影| 久久久久久国产a免费观看| 美女高潮到喷水免费观看| 亚洲精品美女久久av网站| 日本精品一区二区三区蜜桃| 亚洲成人国产一区在线观看| 久久中文看片网| 日韩欧美三级三区| 亚洲成人国产一区在线观看| 18禁国产床啪视频网站| 黄片小视频在线播放| 美女高潮喷水抽搐中文字幕| 日韩欧美三级三区| 免费看十八禁软件| 午夜免费激情av| 日韩高清综合在线| 成人手机av| 精品欧美国产一区二区三| 后天国语完整版免费观看| 91国产中文字幕| 男人操女人黄网站| 99香蕉大伊视频| 18禁黄网站禁片午夜丰满| 国产极品粉嫩免费观看在线| 精品一区二区三区四区五区乱码| 久久久精品国产亚洲av高清涩受| 日本撒尿小便嘘嘘汇集6| 免费不卡黄色视频| 欧美黄色淫秽网站| 欧美一级毛片孕妇| 91麻豆av在线| 多毛熟女@视频| 久久热在线av| 国产激情欧美一区二区| 久久 成人 亚洲| 色播在线永久视频| 99国产精品免费福利视频| 久久国产精品男人的天堂亚洲| avwww免费| 欧美激情 高清一区二区三区| 亚洲成av片中文字幕在线观看| 99国产精品免费福利视频| 国产一区二区激情短视频| 天天添夜夜摸| 亚洲七黄色美女视频| 久久久久久久久久久久大奶| 99riav亚洲国产免费| 国产精华一区二区三区| 久久九九热精品免费| 制服诱惑二区| 18禁黄网站禁片午夜丰满| 亚洲成国产人片在线观看| 午夜福利成人在线免费观看| 午夜久久久久精精品| 亚洲 欧美一区二区三区| 在线观看舔阴道视频| 亚洲最大成人中文| 欧美成人午夜精品| 免费在线观看影片大全网站| 91字幕亚洲| 精品久久久久久,| 婷婷六月久久综合丁香| 免费在线观看完整版高清| av有码第一页| 18禁美女被吸乳视频| 久久天躁狠狠躁夜夜2o2o| tocl精华| 免费在线观看亚洲国产| 校园春色视频在线观看| av天堂久久9| 日日夜夜操网爽| 免费在线观看视频国产中文字幕亚洲| 一级a爱视频在线免费观看| 在线观看www视频免费| 韩国精品一区二区三区| 欧美在线黄色| 岛国视频午夜一区免费看| 成人免费观看视频高清| 亚洲av熟女| 99在线视频只有这里精品首页| 99久久精品国产亚洲精品| 成熟少妇高潮喷水视频| 女性被躁到高潮视频| avwww免费| 久久久久国内视频| 首页视频小说图片口味搜索| 欧美性长视频在线观看| 免费少妇av软件| 搡老妇女老女人老熟妇| 久久精品91无色码中文字幕| 别揉我奶头~嗯~啊~动态视频| 国产97色在线日韩免费| 伦理电影免费视频| 男男h啪啪无遮挡| 手机成人av网站| 97人妻精品一区二区三区麻豆 | 免费观看精品视频网站| 久久精品亚洲精品国产色婷小说| 国产精品98久久久久久宅男小说| 99久久久亚洲精品蜜臀av| 美女大奶头视频| 免费在线观看完整版高清| 成在线人永久免费视频| 色播亚洲综合网| 久99久视频精品免费| 色av中文字幕| 啦啦啦韩国在线观看视频| 搞女人的毛片| 老司机靠b影院| 美女大奶头视频| 国产精品一区二区三区四区久久 | 午夜精品国产一区二区电影| 国产xxxxx性猛交| АⅤ资源中文在线天堂| 一区在线观看完整版| 久久影院123| 亚洲国产日韩欧美精品在线观看 | 国内毛片毛片毛片毛片毛片| 免费看十八禁软件| 精品国产一区二区久久| 久久久久精品国产欧美久久久| 麻豆一二三区av精品| 变态另类成人亚洲欧美熟女 | 国产精华一区二区三区| 女人被狂操c到高潮| 精品久久久久久久久久免费视频| 免费在线观看视频国产中文字幕亚洲| 国产一级毛片七仙女欲春2 | 久久久久久国产a免费观看| 在线观看免费午夜福利视频| 午夜久久久在线观看| 精品久久久久久久人妻蜜臀av | 涩涩av久久男人的天堂| 99热只有精品国产| 一夜夜www| 免费在线观看黄色视频的| 精品久久久久久成人av| 麻豆一二三区av精品| 亚洲专区国产一区二区| 涩涩av久久男人的天堂| 十八禁网站免费在线| 母亲3免费完整高清在线观看| 极品人妻少妇av视频| 成人18禁高潮啪啪吃奶动态图| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品一区av在线观看| 99热只有精品国产| 国产成人精品无人区| 国产麻豆成人av免费视频| 亚洲第一电影网av| 国产免费男女视频| 91字幕亚洲| 人人妻人人爽人人添夜夜欢视频| 18禁美女被吸乳视频| 午夜老司机福利片| 精品久久久久久久久久免费视频| 亚洲午夜精品一区,二区,三区| 成年人黄色毛片网站| 麻豆国产av国片精品| 久久伊人香网站| 极品人妻少妇av视频| 国产日韩一区二区三区精品不卡| 国产激情欧美一区二区| 久久精品91无色码中文字幕| 日韩一卡2卡3卡4卡2021年| 一卡2卡三卡四卡精品乱码亚洲| 日韩av在线大香蕉| bbb黄色大片| 国产在线精品亚洲第一网站| 国产精品国产高清国产av| 国产亚洲精品第一综合不卡| 在线天堂中文资源库| 十八禁人妻一区二区| 亚洲第一欧美日韩一区二区三区| 黄色毛片三级朝国网站| 日本黄色视频三级网站网址| 久久精品aⅴ一区二区三区四区| 久久人妻福利社区极品人妻图片| 18禁黄网站禁片午夜丰满| 黑人欧美特级aaaaaa片| 老熟妇乱子伦视频在线观看| 日韩成人在线观看一区二区三区| 亚洲无线在线观看| 99香蕉大伊视频| 国产成人av激情在线播放| 一级a爱视频在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产三级黄色录像| tocl精华| 黄片播放在线免费| 国产男靠女视频免费网站| 大型av网站在线播放| 亚洲欧美激情综合另类| 国产1区2区3区精品| 成人av一区二区三区在线看| 久久中文看片网| 在线观看午夜福利视频| 国产精品久久电影中文字幕| 757午夜福利合集在线观看| 色在线成人网| 国产精品av久久久久免费| 国产欧美日韩综合在线一区二区| 50天的宝宝边吃奶边哭怎么回事| 亚洲专区国产一区二区| 日本 av在线| 久久亚洲真实| 法律面前人人平等表现在哪些方面| 欧美性长视频在线观看| 久久久国产欧美日韩av| 极品教师在线免费播放| 性少妇av在线| 亚洲国产日韩欧美精品在线观看 | 少妇裸体淫交视频免费看高清 | 黑人欧美特级aaaaaa片| 亚洲三区欧美一区| 亚洲av成人av| 成人三级做爰电影| 久久精品亚洲熟妇少妇任你| 国产精品久久久av美女十八| 中文字幕另类日韩欧美亚洲嫩草| 桃色一区二区三区在线观看| 国产熟女xx| 91成人精品电影| 欧美激情高清一区二区三区| 日韩大尺度精品在线看网址 | 成人18禁高潮啪啪吃奶动态图| 在线观看一区二区三区| 18禁黄网站禁片午夜丰满| 曰老女人黄片| 色综合站精品国产| 91大片在线观看| 国产主播在线观看一区二区| 精品一区二区三区四区五区乱码| 女人被躁到高潮嗷嗷叫费观| 久久久久国内视频| 欧美日本视频| 美女午夜性视频免费| 一进一出抽搐动态| 在线观看免费视频日本深夜| 亚洲av熟女| 亚洲午夜理论影院| 欧美黑人欧美精品刺激| 级片在线观看| 国产精品久久久av美女十八| 国产成人av教育| 亚洲成人久久性| 在线观看免费日韩欧美大片| 亚洲色图av天堂| 国产黄a三级三级三级人| 高清在线国产一区| 最好的美女福利视频网| 国产av又大| 9热在线视频观看99| 亚洲精品粉嫩美女一区| 中文字幕色久视频| 天天躁夜夜躁狠狠躁躁| 丝袜人妻中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲激情在线av| 99久久国产精品久久久| 久久国产精品影院| 国产午夜精品久久久久久| 国产单亲对白刺激| ponron亚洲| 一级,二级,三级黄色视频| 久久精品亚洲精品国产色婷小说| 国产成人精品在线电影| 99国产精品99久久久久| АⅤ资源中文在线天堂| 色尼玛亚洲综合影院| 国产不卡一卡二| 精品一区二区三区四区五区乱码| 妹子高潮喷水视频| 桃红色精品国产亚洲av| 久久国产精品影院| 夜夜爽天天搞| 国产精品乱码一区二三区的特点 | av电影中文网址| 国产一区二区三区综合在线观看| 亚洲国产中文字幕在线视频| 黄色 视频免费看| 亚洲五月天丁香| 超碰成人久久| 一二三四社区在线视频社区8| 亚洲第一欧美日韩一区二区三区| 99国产极品粉嫩在线观看| 亚洲成av人片免费观看| 欧美精品亚洲一区二区| 亚洲aⅴ乱码一区二区在线播放 | 日韩欧美一区二区三区在线观看| www国产在线视频色| 国产精品综合久久久久久久免费 | svipshipincom国产片| 亚洲少妇的诱惑av| 好男人在线观看高清免费视频 | 午夜精品在线福利| 操美女的视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 禁无遮挡网站| 国产精品一区二区精品视频观看| 国产精品爽爽va在线观看网站 | 亚洲欧美精品综合久久99| 午夜免费成人在线视频| 热99re8久久精品国产| 亚洲av成人不卡在线观看播放网| 久久久精品国产亚洲av高清涩受| 黄色 视频免费看| 亚洲精品久久国产高清桃花| 一级黄色大片毛片| 亚洲欧美一区二区三区黑人| 欧美乱色亚洲激情| 国产欧美日韩精品亚洲av| 天堂影院成人在线观看| 黄色女人牲交| 国产精品 欧美亚洲| 啦啦啦免费观看视频1| av片东京热男人的天堂| 欧美日韩黄片免| 99riav亚洲国产免费| 欧美不卡视频在线免费观看 | 国产成人免费无遮挡视频| 亚洲第一青青草原| 欧美国产日韩亚洲一区| 亚洲国产看品久久| 人妻久久中文字幕网| 国产不卡一卡二| 亚洲国产中文字幕在线视频| 久久亚洲精品不卡| 亚洲av电影不卡..在线观看| 久久亚洲精品不卡| 制服人妻中文乱码| 日本精品一区二区三区蜜桃| 中文字幕av电影在线播放| 精品人妻在线不人妻| 亚洲第一电影网av| 国产91精品成人一区二区三区| 欧美大码av| 黄片小视频在线播放| 欧美久久黑人一区二区| 久久天堂一区二区三区四区| 国产人伦9x9x在线观看| 99精品欧美一区二区三区四区| 久久国产精品人妻蜜桃| 日本精品一区二区三区蜜桃| 久久人人爽av亚洲精品天堂| 一二三四社区在线视频社区8| 怎么达到女性高潮| 亚洲精品美女久久久久99蜜臀| xxx96com| 十八禁网站免费在线| 亚洲欧美日韩高清在线视频| 十八禁网站免费在线| 久久人人爽av亚洲精品天堂| 美女高潮喷水抽搐中文字幕| 黄片播放在线免费| 91精品三级在线观看| 精品一区二区三区av网在线观看| 午夜福利,免费看| 国产成人av激情在线播放| www.999成人在线观看| 啦啦啦观看免费观看视频高清 | 最新美女视频免费是黄的| 精品乱码久久久久久99久播| 久久午夜综合久久蜜桃| 久久中文看片网| 国内精品久久久久精免费| 国产午夜精品久久久久久| 99国产精品99久久久久| 亚洲av电影在线进入| 一本综合久久免费| 日本a在线网址| 国产精品国产高清国产av| 电影成人av| 亚洲欧美日韩另类电影网站| 一边摸一边抽搐一进一小说| 午夜福利高清视频| 亚洲国产精品久久男人天堂| 淫秽高清视频在线观看| 精品不卡国产一区二区三区| 日本撒尿小便嘘嘘汇集6| 在线观看一区二区三区| 看片在线看免费视频| 亚洲五月天丁香| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品国产色婷婷电影| 国产黄a三级三级三级人| 午夜日韩欧美国产| 免费看美女性在线毛片视频| 男男h啪啪无遮挡| 久热爱精品视频在线9| 日韩三级视频一区二区三区| 精品久久久久久久人妻蜜臀av | 欧美在线一区亚洲| 女同久久另类99精品国产91| 精品欧美国产一区二区三| 国产高清视频在线播放一区| 亚洲欧美激情综合另类| 在线观看一区二区三区| 国产在线精品亚洲第一网站| 久久久国产精品麻豆| 欧美激情久久久久久爽电影 | 999久久久精品免费观看国产| 久久久精品国产亚洲av高清涩受| 亚洲第一av免费看| 国产精品一区二区免费欧美| 国产av精品麻豆| 成年版毛片免费区| 久久久国产成人免费| 亚洲国产精品sss在线观看| 伦理电影免费视频| 欧美激情极品国产一区二区三区| 国产激情欧美一区二区| 亚洲精品国产精品久久久不卡| 91在线观看av| 亚洲七黄色美女视频| 老熟妇乱子伦视频在线观看| 日韩高清综合在线| 女人被躁到高潮嗷嗷叫费观| 精品无人区乱码1区二区| 国产在线观看jvid| 成人18禁在线播放| 长腿黑丝高跟| 在线观看舔阴道视频| 国产成人精品无人区| 日韩中文字幕欧美一区二区| 可以在线观看的亚洲视频| 亚洲成人精品中文字幕电影| 熟女少妇亚洲综合色aaa.| 国产欧美日韩一区二区三| 亚洲 欧美一区二区三区| 极品教师在线免费播放| 欧美不卡视频在线免费观看 | 国产在线观看jvid| 高清在线国产一区| 国产黄a三级三级三级人| 无限看片的www在线观看| 国产精品久久视频播放| 成人av一区二区三区在线看| 久久久国产成人精品二区| 午夜免费观看网址| 在线国产一区二区在线| 国产伦人伦偷精品视频| 亚洲熟女毛片儿| 成人三级做爰电影| 一个人免费在线观看的高清视频| 激情视频va一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 90打野战视频偷拍视频| 999久久久精品免费观看国产| 国产成人啪精品午夜网站| 叶爱在线成人免费视频播放| 成人av一区二区三区在线看| 亚洲精品国产色婷婷电影| 亚洲七黄色美女视频|