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

    基于等效電路模型的鈣鈦礦太陽電池效率損失機(jī)理分析*

    2021-05-14 02:45:12徐婷王子帥李炫華沙威
    物理學(xué)報(bào) 2021年9期
    關(guān)鍵詞:傳輸層等效電路太陽電池

    徐婷 王子帥 李炫華 沙威?

    1) (浙江大學(xué)信息與電子工程學(xué)院, 現(xiàn)代光學(xué)儀器國家重點(diǎn)實(shí)驗(yàn)室, 杭州 310027)

    2) (香港大學(xué)電機(jī)電子工程系, 香港 999077)

    3) (西北工業(yè)大學(xué)材料學(xué)院, 凝固技術(shù)國家重點(diǎn)實(shí)驗(yàn)室, 西安 710072)

    理解并量化影響鈣鈦礦太陽電池效率的因素, 對研發(fā)高性能器件尤為重要.目前, 太陽電池普遍認(rèn)可的三大損失為光學(xué)損失、歐姆損失和非輻射復(fù)合損失.其中, 非輻射復(fù)合包括體復(fù)合和表面復(fù)合, 已被證明是制約電池效率提升的決定性因素.本文提出了一種分析電池伏安特性曲線的等效電路模型, 能對上述損失機(jī)制進(jìn)行全面描述, 并通過與漂移-擴(kuò)散模型及實(shí)驗(yàn)結(jié)果的對比, 證實(shí)了電路模型的可靠性, 擬合誤差在2%以內(nèi).根據(jù)該模型, 可以準(zhǔn)確判斷電池內(nèi)的主導(dǎo)復(fù)合機(jī)制, 并可從實(shí)際電池伏安曲線中提取不同效率損失對應(yīng)的物理參數(shù), 繪制電壓掃描過程中各機(jī)制隨電壓的演化曲線, 從而理解效率損失的物理機(jī)理.該模型從電路角度分析了不同損失機(jī)制對電池工作特性的影響, 有助于定位提高效率的關(guān)鍵點(diǎn), 是一個較全面的鈣鈦礦太陽電池仿真分析工具.

    1 引 言

    近年來, 隨著鈣鈦礦太陽電池研究的深入, 電池內(nèi)光電轉(zhuǎn)換效率的損失機(jī)制分析引起了學(xué)術(shù)界和產(chǎn)業(yè)界的廣泛關(guān)注.在理想的單節(jié)鈣鈦礦太陽電池中, 電子和空穴僅僅通過輻射復(fù)合發(fā)出光子, 其理論效率極限, 即Shockley-Queisser 極限, 約為31%[1].而對于實(shí)際太陽電池, 目前已實(shí)驗(yàn)認(rèn)證的最高光電轉(zhuǎn)換效率為25.5%[2], 距理論極限值仍相差甚遠(yuǎn).究其原因, 一方面, 鈣鈦礦中自由載流子的輻射復(fù)合較弱[3]; 另一方面, 電池內(nèi)各種損失機(jī)制的存在, 極大影響了載流子的產(chǎn)生、輸運(yùn)與收集,最終導(dǎo)致光電轉(zhuǎn)換效率降低.

    基于細(xì)致平衡理論, 我們之前的工作報(bào)道了對于實(shí)際的鈣鈦礦太陽電池主要有三種損失機(jī)制限制了電池的效率: 一是光學(xué)損失, 二是體缺陷輔助非輻射(SRH)復(fù)合損失, 三是歐姆損失[4].考慮到實(shí)際電池會受到器件陷光結(jié)構(gòu)和光子循環(huán)效應(yīng)的影響, 根據(jù)提出的修正細(xì)致平衡模型, 量化了影響電池工作的主要損失機(jī)制, 并成功將理論和實(shí)驗(yàn)的伏安曲線(J-V曲線)擬合誤差控制在4%以內(nèi).但該模型沒有考慮或區(qū)分其他類型的非輻射復(fù)合對器件效率的影響.對于鈣鈦礦太陽電池, 非輻射復(fù)合已經(jīng)被證明是主要復(fù)合機(jī)制[5,6].目前, 廣泛認(rèn)為鈣鈦礦太陽電池內(nèi)的非輻射復(fù)合途徑主要有體SRH 復(fù)合、體俄歇復(fù)合和表面復(fù)合.其中, 由于鈣鈦礦材料具有俄歇復(fù)合率低的特點(diǎn)[7], 對電池建模時常常忽略俄歇復(fù)合.體復(fù)合主要與鈣鈦礦光吸收層內(nèi)的固有點(diǎn)缺陷、雜質(zhì)等有關(guān), 而表面復(fù)合主要與鈣鈦礦層的表面缺陷[8]相關(guān).因此判斷電池中的主導(dǎo)非輻射復(fù)合機(jī)制, 分析并量化體復(fù)合與表面復(fù)合對特性曲線的影響, 對電池的效率和穩(wěn)定性的提高具有重要意義[9,10].

    現(xiàn)有的分析非輻射復(fù)合類型的方法主要分為測試方法和仿真方法.測試方法主要有理想因子判定法[11]和微擾測量法, 如阻抗譜分析法[12?14]、瞬態(tài)光電壓測試法[15,16]與光譜測量法[17]等.但理想因子法無法分析表面復(fù)合的影響, 微擾法難以量化正常工作時電池內(nèi)的體復(fù)合損失和表面復(fù)合損失,并且需要昂貴的設(shè)備支持.而仿真方法可基于漂移-擴(kuò)散模型, 結(jié)合體復(fù)合公式和表面復(fù)合公式, 以仿真太陽電池J-V曲線.以往的漂移-擴(kuò)散模型[18]對光學(xué)損失的探討較少且缺乏與細(xì)致平衡模型之間的聯(lián)系, 直到2017 年, Ren 等[19]將考慮陷光結(jié)構(gòu)和輻射角度限制的Roosbroeck-Shockley 公式,以及選擇性電極邊界條件應(yīng)用于漂移-擴(kuò)散方程中,才得出了漂移-擴(kuò)散模型與細(xì)致平衡模型的等價條件.不幸的是, 漂移-擴(kuò)散模型由于過強(qiáng)的非線性特征, 不適用于擬合實(shí)驗(yàn)曲線, 難以提取物理參數(shù)從而理解實(shí)際電池的工作機(jī)理.

    為解決上述問題, 本文基于修正細(xì)致平衡模型[4], 同時全面考慮光學(xué)損失、歐姆損失、體復(fù)合損失和表面復(fù)合損失, 提出了新的等效電路模型, 以描述鈣鈦礦太陽電池在實(shí)際情況下的工作特性, 分析效率損失機(jī)理.同時, 為驗(yàn)證模型的可靠性, 本文將電路模型與漂移-擴(kuò)散仿真的太陽電池J-V曲線、以及實(shí)驗(yàn)J-V曲線進(jìn)行擬合并對比, 通過特征參數(shù)(串聯(lián)電阻、并聯(lián)電阻、體復(fù)合系數(shù)和表面復(fù)合系數(shù))的提取與數(shù)值比較, 評估等效電路模型區(qū)分非輻射復(fù)合類型的能力, 以及量化各損失機(jī)制的效果.

    2 理 論

    2.1 等效電路模型

    光照下, 理想的光伏電池可看作一個理想二極管和一個恒流源組成的并聯(lián)電路, 恒流源產(chǎn)生的電流Jph即為光生電流.基于修正細(xì)致平衡模型[4], 描述實(shí)際鈣鈦礦太陽電池伏安曲線的改進(jìn)等效電路模型如圖1 所示.

    圖1 鈣鈦礦太陽電池的等效電路模型圖Fig.1.Equivalent circuit model of perovskite solar cells.

    在制造過程中, 實(shí)際的太陽電池會出現(xiàn)陷阱和孔洞等缺陷, 這些缺陷造成的漏電流損失用并聯(lián)電阻Rsh來表示.正負(fù)電極、載流子傳輸層、鈣鈦礦層與載流子傳輸層界面產(chǎn)生的歐姆損失的總效果可用串聯(lián)電阻Rs刻畫.理想二極管Dr代表鈣鈦礦的輻射復(fù)合, 而非輻射復(fù)合過程則由二極管Dnr描述.

    模型對應(yīng)的伏安特性方程為

    其中,V為太陽電池的光生電壓,Jr為由光子循環(huán)引起的輻射復(fù)合電流,Jbulk為考慮體復(fù)合的非輻射復(fù)合電流,Jsurf為表面復(fù)合電流,Jsh為并聯(lián)電阻電流(漏電流).

    光電流由以下方程給出:

    其中,c0代表光速,Γ是AM 1.5 G 的太陽光譜,λ代表波長,q是基本電荷.吸收率α等于鈣鈦礦活性層吸收的功率與入射太陽光的功率之比, 其值取決于鈣鈦礦層的厚度、選取的相關(guān)材料的折射率和器件陷光結(jié)構(gòu)的設(shè)計(jì)等.光電流可通過數(shù)值求解麥克斯韋方程得到.

    基于細(xì)致平衡的輻射電流表示為

    其中,Γ0是太陽電池T=300 K 時的黑體輻射譜,kB是玻爾茲曼常數(shù).

    忽略俄歇復(fù)合, 鈣鈦礦太陽電池非輻射復(fù)合的主導(dǎo)機(jī)制為體復(fù)合和表面復(fù)合.因此, 體復(fù)合電流和表面復(fù)合電流可分別表述為

    其中,Jbulk代表體復(fù)合電流,Jsurf代表表面復(fù)合電流,γbulk和γsurf分別為體復(fù)合系數(shù)和表面復(fù)合系數(shù),L是鈣鈦礦活性層的厚度,Lsurf是傳輸層與鈣鈦礦層界面的有效厚度,ni是活性材料的本征載流子濃度, 而是空穴傳輸層/鈣鈦礦層界面處鈣鈦礦側(cè)的平衡多子空穴濃度(或電子傳輸層/鈣鈦礦層界面處鈣鈦礦側(cè)的平衡多子電子濃度).(4)式和(5)式的具體推導(dǎo)見附錄A.

    漏電流Jsh由(6)式描述:

    基于該改進(jìn)的等效電路模型, 通過對實(shí)際J-V曲線擬合, 提取出串聯(lián)電阻Rs、并聯(lián)電阻Rsh、非輻射體復(fù)合系數(shù)γbulk和非輻射表面復(fù)合系數(shù)γsurf四個參數(shù), 進(jìn)而分析造成鈣鈦礦太陽電池效率損失的主要因素, 以及各損失機(jī)制(串聯(lián)電阻損耗、并聯(lián)電阻損耗、體復(fù)合損耗和表面復(fù)合損耗)對效率損失的貢獻(xiàn)比例.除此之外, 通過Jbulk-V,Jsurf-V和Jsh-V曲線的對比, 可以直觀分析隨著電壓變化,體復(fù)合、表面復(fù)合以及電阻損耗的演進(jìn)情況, 從而實(shí)現(xiàn)從電路角度分析鈣鈦礦太陽電池的工作機(jī)理,理解損失機(jī)制的目的, 為進(jìn)一步提高鈣鈦礦太陽電池的光電轉(zhuǎn)換效率提供更精確的方向.

    2.2 仿真方法

    為了驗(yàn)證等效電路模型分辨非輻射復(fù)合機(jī)制的可靠性與準(zhǔn)確性, 本文將根據(jù)漂移-擴(kuò)散模型[17]仿真得到的三條伏安曲線作為參考曲線.其中, 第一條J-V曲線代表實(shí)際太陽電池工作時非輻射復(fù)合機(jī)制僅存在體復(fù)合, 第二條J-V曲線描述了表面復(fù)合為主導(dǎo)非輻射復(fù)合類型時太陽電池的工作特性, 而第三條曲線則代表沒有非輻射復(fù)合但改變兩側(cè)傳輸層遷移率時對應(yīng)的太陽電池J-V曲線.

    描述鈣鈦礦太陽電池器件工作特性的漂移-擴(kuò)散模型主要由泊松方程、漂移-擴(kuò)散方程和電流連續(xù)性方程給出.泊松方程如下式所示:

    電子和空穴的漂移-擴(kuò)散方程為

    電流連續(xù)性方程為

    其中,ψ為電勢,q為基本電荷,n和p分別為電子密度和空穴密度,Jn和Jp分別為電子電流密度和空穴電流密度,μn,μp,Dn,Dp分別為電子遷移率、空穴遷移率、電子擴(kuò)散系數(shù)和空穴擴(kuò)散系數(shù).G為產(chǎn)生率,R代表復(fù)合率.

    本文考慮的復(fù)合機(jī)制主要有輻射復(fù)合和非輻射復(fù)合.非輻射復(fù)合又主要考慮體復(fù)合與活性層-載流子傳輸層界面處的表面復(fù)合.故輻射復(fù)合率Rrad、體復(fù)合率Rbulk和表面復(fù)合率Rsurf分別為

    其中,ni為鈣鈦礦活性層本征載流子濃度;krad為輻射復(fù)合系數(shù);τn,τp分別為電子、空穴的體復(fù)合壽命;τsurfn,τsurfp分別為電子、空穴的表面復(fù)合壽命;nt,pt分別代表陷阱電子、空穴濃度;n+,p?分別代表界面兩側(cè)的電子、空穴密度.

    為減小電極表面復(fù)合對J-V曲線的影響, 選用以下選擇性接觸邊界條件:

    其中Snc,Sna,Spc,Spa分別為陰極電子、陽極電子、陰極空穴和陽極空穴的表面復(fù)合速率.n0c,n0a,p0c,p0a為邊界處陰極電子、陽極電子、陰極空穴和陽極空穴的濃度.

    3 仿真結(jié)果分析

    根據(jù)2.2 節(jié)中的漂移-擴(kuò)散模型, 假設(shè)鈣鈦礦太陽電池內(nèi)的非輻射復(fù)合機(jī)制只有體復(fù)合, 或只有表面復(fù)合, 或不考慮非輻射復(fù)合但改變活性層兩側(cè)傳輸層的遷移率大小.若假設(shè)僅存在非輻射體復(fù)合, 即 令τbulk=100 ns ,τsurf=Inf[20], 傳輸層多子遷移率μmaj=20 cm2/(V·s) , 傳輸層少子遷移率μmin≈0 cm2/(V·s)[21].若假設(shè)僅存在表面復(fù)合, 則令τbulk=Inf ,τsurf=1 ns ,μmaj=20 cm2/(V·s) ,μmin≈0 cm2/(V·s).若假設(shè)不存在非輻射復(fù)合并且傳輸層遷移率低, 則令τbulk=Inf ,τsurf=Inf ,μmaj=1×10?2cm2/(V·s),μmin=1×10?5cm2/(V·s).在無初始偏置電壓的情況下, 用 0.1 V/s 的速度從短路( 0 V )掃描到稍大于開路(Voc), 得到對應(yīng)的實(shí)際鈣鈦礦太陽電池仿真J-V特性曲線, 如圖2 中紅色點(diǎn)劃線所示.漂移-擴(kuò)散模型中使用的其他仿真參數(shù)主要有: 器件厚度ETL 200 nm/Absorber 500 nm/HTL 200 nm[22]; 電子傳輸層、空穴傳輸層以及活性層的相對介電常數(shù)分別為εrETL=4 ,εrHTL=4 ,εrAbsorber= 31[23]; 設(shè)活性層電子和空穴的遷移率相等, 為μn=μp=20 c m2/(V·s)[21]; 內(nèi)建電場為Vbi=0.8 V[24]; 禁帶寬度為Eg=1.6 eV[25,26];鈣鈦礦層的態(tài)密度為Nc=Nv= 1 019/cm3[22]; 鈣鈦礦層和傳輸層之間的有效界面厚度為Lsurf=0.02 nm[20].故模型中的鈣鈦礦層本征載流子濃度為

    觀察圖2, 根據(jù)2.1 節(jié)中的等效電路模型, 分別對只含體復(fù)合(圖2(a)中紅色點(diǎn)劃線)、表面復(fù)合占優(yōu)(圖2(b) 中紅色點(diǎn)劃線)以及不考慮非輻射復(fù)合但改變傳輸層(圖2(c)中的紅色點(diǎn)劃線)的驗(yàn)證曲線進(jìn)行擬合, 得到對應(yīng)的黑色擬合曲線.發(fā)現(xiàn)三種情況下均有較好的擬合結(jié)果, 其中只含體復(fù)合時曲線的相對擬合誤差為1.01%, 表面復(fù)合主導(dǎo)時的相對擬合誤差僅為0.2%, 而改變傳輸層時相對擬合誤差也僅為0.31%.極小的擬合誤差, 證明了該等效電路模型可以較為準(zhǔn)確地描述鈣鈦礦太陽電池的伏安特性.

    同時, 提取擬合參數(shù)并進(jìn)行光電轉(zhuǎn)換效率分析, 得到的結(jié)果如表1 所列.其中, 經(jīng)漂移-擴(kuò)散模型仿真, 空穴傳輸層/鈣鈦礦界面處、靠近鈣鈦礦側(cè)的平衡多子空穴濃度均為=9.29×1010/cm3,同樣靠近電子傳輸層/鈣鈦礦界面處鈣鈦礦側(cè)的平衡多子電子濃度均為=9.29×1010/cm3.

    由于本文的漂移-擴(kuò)散模型與等效電路模型具有統(tǒng)一性, 因此對于同一個鈣鈦礦太陽電池參考曲線來說, 其體非輻射復(fù)合系數(shù)和表面復(fù)合系數(shù)應(yīng)分別與體載流子壽命和表面載流子壽命相關(guān).根據(jù)附錄A 可知,τ=1/γ, 故計(jì)算兩個模型仿真得到的非輻射復(fù)合參數(shù)(表2), 以比較兩個模型間非輻射復(fù)合系數(shù)和載流子壽命的關(guān)系.

    圖2 不同缺陷類型和傳輸層遷移率對應(yīng)的鈣鈦礦太陽電池 J -V 曲線圖 (a) 非輻射復(fù)合機(jī)制僅考慮體復(fù)合; (b) 表面復(fù)合為主導(dǎo)非輻射復(fù)合機(jī)制; (c) 不考慮非輻射復(fù)合且改變傳輸層遷移率.其中, 紅色點(diǎn)劃線為漂移-擴(kuò)散模型仿真得到的 J -V 曲線, 而黑色實(shí)線為經(jīng)等效電路模型擬合得到的 J -V 曲線Fig.2.The J -V curves of perovskite solar cells with different non-radiative recombination types and different transport layers:(a) Only bulk recombination is considered; (b) surface recombination is the dominant non-radiative recombination mechanism;(c) without non-radiative recombination and the mobility of transport layers is changed.The red-dot lines represent J -V curves that are simulated by drift-diffusion model, and the curves fitted by equivalent circuit model are shown in the dark solid lines.

    表1 不同情況下鈣鈦礦太陽電池J -V 曲線對應(yīng)的特征參數(shù)表Table 1.Parameters retrieved from the J -V curves of different cases.

    表2 不同情況下經(jīng)等效電路模型和漂移-擴(kuò)散模型仿真得到的非輻射復(fù)合參數(shù)表Table 2.Nonradiative recombination parameters retrieved from different cases by equivalent circuit model and drift-diffusion model.

    僅考慮體復(fù)合時, 漂移-擴(kuò)散模型設(shè)置τbulk=100 ns ,τsurf=Inf , 故其對應(yīng)的體復(fù)合系數(shù)計(jì)算值應(yīng)為 1 07s?1, 而等效電路模型的擬合值是2.07×106s?1, 說明該情況下體復(fù)合系數(shù)的擬合誤差在可接受范圍內(nèi)(約為計(jì)算值的 1 /5 ).對應(yīng)的表面復(fù)合系數(shù)γsurf計(jì)算值應(yīng)為無限小, 但實(shí)際等效電路模型的擬合值是 3.48×105s?1.此時, 根據(jù)表面復(fù)合系數(shù)擬合值相對體復(fù)合系數(shù)擬合值顯著小(相差一個數(shù)量級), 可推斷該鈣鈦礦太陽電池內(nèi)表面復(fù)合較弱.僅考慮表面復(fù)合時, 設(shè)置τbulk=Inf ,τsurf=1 ns ,其體復(fù)合系數(shù)的計(jì)算值應(yīng)為無限小, 表面復(fù)合系數(shù)的計(jì)算值應(yīng)為 1 09s?1, 而等效電路模型的擬合值為γbulk= 1.30×107s?1,γsurf=1.95×109s?1.同樣, 該情況下表面復(fù)合系數(shù)的擬合誤差在允許范圍內(nèi)(約為計(jì)算值的2 倍), 而體復(fù)合系數(shù)擬合值相比表面復(fù)合系數(shù)擬合值顯著小, 可認(rèn)為鈣鈦礦中體復(fù)合較弱.當(dāng)改變傳輸層遷移率時, 體復(fù)合系數(shù)和表面復(fù)合系數(shù)的計(jì)算值應(yīng)均為無限小, 而對應(yīng)的等效電路擬合值分別為γbulk=8.75×104s?1,γsurf=0.86 s?1.對比前兩個例子, 可以明顯看出該情況下, 體復(fù)合系數(shù)擬合值和表面復(fù)合系數(shù)擬合值都相對較小, 即非輻射復(fù)合系數(shù)擬合值與漂移-擴(kuò)散模型的壽命設(shè)定值之間有較好的對應(yīng)關(guān)系.綜上, 在已知Lsurf和的情況下, 根據(jù)J-V曲線擬合得到的γbulk和γsurf的大小比較, 可以粗略判斷該電池內(nèi)的非輻射復(fù)合情況.然而, 單純比較γbulk和γsurf的數(shù)值大小, 難以直觀理解鈣鈦礦太陽電池中的非輻射復(fù)合機(jī)理, 因此需要更形象的比較方法.

    除此之外, 分析并比較表1 中的結(jié)果, 發(fā)現(xiàn)三種情況下的短路電流大小相近.僅考慮表面復(fù)合時的開路電壓最小, 不考慮非輻射復(fù)合但改變傳輸層遷移率時的效率最高但填充因子最小.另外, 僅考慮體復(fù)合或僅考慮表面復(fù)合時的Rs都很小,Rsh都很大.而降低傳輸層遷移率時,Rs相對以上兩種情況最大,Rsh最小.同樣, 單純觀察Rs和Rsh的數(shù)值大小, 也難以分析歐姆損失對工作特性的影響.故根據(jù)(1)式作出電池總電流和各子電流隨電壓變化的曲線, 如圖3 所示.

    為進(jìn)一步分析串聯(lián)電阻、并聯(lián)電阻、體復(fù)合和表面復(fù)合對鈣鈦礦太陽電池光電轉(zhuǎn)換效率的影響, 繪制效率損失示意圖, 如圖4 所示(繪制方法見附錄B).

    觀察圖3 和圖4 可知, 僅考慮體復(fù)合和僅考慮表面復(fù)合時的電阻電流都極小, 并且串并聯(lián)電阻造成的效率損失均為 0.眾所周知, 串聯(lián)電阻損失的主要來源為載流子傳輸層及其界面的歐姆損耗.對應(yīng)于該兩種情況的漂移-擴(kuò)散模型參數(shù)為μmaj=20 cm2/(V·s),μmin≈0 cm2/(V·s) , 即此時鈣鈦礦活性層兩側(cè)均為高電導(dǎo)多子傳輸層(少子阻擋層).故該電池的歐姆損耗應(yīng)極小, 也就是串聯(lián)電阻損失應(yīng)極小.該分析與串聯(lián)電阻造成的效率損失約為0%相統(tǒng)一.并且當(dāng)降低傳輸層遷移率后μmaj=1×10?2cm2/(V·s), 因?yàn)閭鬏攲幼杩乖龃? 其串聯(lián)電阻損失增加為42.16%.另外, 擴(kuò)散-漂移模型不能刻畫因電池內(nèi)缺陷和針孔造成的漏電流損失(漏電流橫向傳輸, 不被上下電極收集).因此由等效電路得出的并聯(lián)阻抗損失幾乎為0.綜上, 該等效電路模型可以較好地描述歐姆電阻對伏安特性的影響.

    圖3 根據(jù)(1)式分解的不同情況下的鈣鈦礦太陽電池電流組成示意圖 (a), (d) 僅考慮體復(fù)合; (b), (e) 非輻射復(fù)合以表面復(fù)合為主; (c), (f)不考慮非輻射復(fù)合但改變傳輸層.其中 J 代表鈣鈦礦太陽電池的總電流, J bulk 為體復(fù)合電流, J surf 為表面復(fù)合電流, J sh 為電阻電流Fig.3.Decompositions of the total current density of perovskite solar cells according to Eq.(1): (a), (d) Only bulk recombination is considered; (b), (e) only surface recombination is considered; (c), (f) without non-radiative recombination and with different transport layers.J represents the total current, Jbulk represents the bulk recombination current and Jsurf represents the surface recombination current.J sh represents the resistance current.

    對于圖3(a), 在掃描電壓從 0 V 增加至Vbi=0.8 V 的過程中, 體復(fù)合電流Jbulk和表面復(fù)合電流Jsurf均變化不顯著.隨后, 在電壓從Vbi繼續(xù)向Voc掃描的過程中,Jbulk和Jsurf均以指數(shù)規(guī)律增大, 但Jsurf開始指數(shù)增長的電壓起始點(diǎn)明顯滯后于Jbulk的起始增長點(diǎn), 導(dǎo)致在該過程中Jbulk始終大于Jsurf.該結(jié)果與漂移-擴(kuò)散模型的預(yù)設(shè)條件τbulk=100 ns ,τsurf=Inf基本一致.同時, 根據(jù)圖4, 計(jì)算的體復(fù)合造成的效率損失占91.18%, 遠(yuǎn)高于表面復(fù)合造成的損失8.82%.說明在該鈣鈦礦太陽電池工作時,兩種非輻射復(fù)合機(jī)制中體復(fù)合占優(yōu), 而表面復(fù)合幾乎可忽略.對于圖3(b), 在掃描電壓從 0 V 增加至Vbi=0.8 V 的過程中,Jbulk和Jsurf同樣變化不顯著.而在電壓從Vbi增加至Voc的過程中,Jbulk和Jsurf均從Vbi開始以指數(shù)規(guī)律增大, 但Jsurf的增長速度明顯快于Jbulk.同時, 根據(jù)圖4, 體復(fù)合造成的效率損失占總損失的22.41%, 表面復(fù)合占77.59%.表明該鈣鈦礦太陽電池工作時, 表面復(fù)合是其主要的非輻射復(fù)合機(jī)制, 是效率損失的決定性因素.該結(jié)論同樣與漂移-擴(kuò)散模型的預(yù)設(shè)條件τbulk=Inf ,τsurf=1 ns 一致.對于圖3(c),Jbulk仍從Vbi開始隨電壓以指數(shù)規(guī)律增大, 但Jsurf始終約為 0.根據(jù)圖4, 該情況下體復(fù)合造成了57.28%的效率損失,而表面復(fù)合沒有造成損失.對應(yīng)于漂移-擴(kuò)散模型,仿真中未設(shè)置非輻射復(fù)合而僅改變了傳輸層遷移率大小, 也就代表此時電池內(nèi)幾乎不含非輻射復(fù)合, 但因?yàn)殡娐纺P驮诜禽椛鋸?fù)合為 0 的極端條件下存在誤差, 所以擬合的結(jié)果仍有體復(fù)合存在.特別地, 在該電池中電阻造成了額外的損耗, 較好地符合了模型對應(yīng)的物理機(jī)理.綜上, 本文提出的等效電路模型可以有效區(qū)分鈣鈦礦太陽電池工作時的主導(dǎo)非輻射復(fù)合機(jī)制, 并且能夠量化不同非輻射復(fù)合類型, 串并聯(lián)電阻對于光電轉(zhuǎn)換效率的影響,進(jìn)而分析影響電池伏安曲線變化的可能原因.

    圖4 不同情況下鈣鈦礦太陽電池的效率損失示意圖Fig.4.Efficiency loss of perovskite solar cells in different cases.

    4 實(shí)驗(yàn)結(jié)果分析

    接下來進(jìn)一步驗(yàn)證本文提出的等效電路模型對實(shí)際鈣鈦礦太陽電池伏安特性的建模能力, 評估其量化伏安特性損失機(jī)制的有效性.基于文獻(xiàn)[27],對在不同晶界處理工程下鈣鈦礦太陽電池的測試數(shù)據(jù)進(jìn)行分析.根據(jù)等效電路模型, 擬合提取的特征參數(shù)如表3 所列, 其中Control 代表未引入路易斯堿或酸官能團(tuán)的MAPbI3鈣鈦礦太陽電池,DTS 代表引入DTS[28]的MAPbI3鈣鈦礦太陽電池, 而DR3T 則為引入基于BDT 的DR3TBDTT(簡稱DR3T)[29]的MAPbI3鈣鈦礦太陽電池.值得注意的是, 在等效電路模型的仿真中, 由于鈣鈦礦太陽電池的Lsurf和未知, 表3 引入U(xiǎn)surf=以描述表面復(fù)合對太陽電池的影響.

    根據(jù)表3 可以發(fā)現(xiàn), 引入DTS 之后, 相比Control, 鈣鈦礦太陽電池的體復(fù)合系數(shù)顯著降低,開路電壓增大, 表明DTS 具有良好的MAPbI3活性層晶界鈍化作用, 可以有效提高鈣鈦礦太陽電池的工作性能.同樣, 引入DR3T 后, 體復(fù)合系數(shù)進(jìn)一步降低, 開路電壓增大, 表明DR3T 的晶界鈍化作用優(yōu)于DTS; 但由于表面復(fù)合的增加, 最終效率提升不明顯.

    圖5 為不同晶界處理情況下, 鈣鈦礦太陽電池伏安特性的等效電路模型擬合結(jié)果, 而圖6 為對應(yīng)的效率損失示意圖.其中圖5(a)為未進(jìn)行晶界修飾時的鈣鈦礦太陽電池等效電路擬合曲線和各子電流的曲線, 其總電流-電壓的理論曲線和實(shí)驗(yàn)曲線的相對擬合誤差為1.08%.圖5(b)為引入DTS時的電流-電壓曲線, 其理論和實(shí)驗(yàn)曲線的相對擬合誤差為0.70%.圖5(c)為引入DR3T 時的曲線,其相對擬合誤差為0.95%.極小的相對誤差, 再一次表明該等效電路模型可以很好地描述實(shí)際鈣鈦礦太陽電池的電流密度-電壓曲線.

    表3 不同情況下鈣鈦礦太陽電池J -V 曲線對應(yīng)的特征參數(shù)表Table 3.Parameters retrieved from the J -V curves of different cases.

    圖5 根據(jù)(1)式分解的不同情況下的鈣鈦礦太陽電池電流組成示意圖 (a) 未進(jìn)行鈣鈦礦層晶界修飾的鈣鈦礦太陽電池器件;(b) 鈣鈦礦層引入DTS 的太陽電池器件; (c)鈣鈦礦層引入DR3T 的器件.其中 J theoretical 代表等效電路模型擬合得到的鈣鈦礦太陽電池的總電流, J bulk 為其體復(fù)合電流, J surf 為表面復(fù)合電流, J experimental 為實(shí)驗(yàn)測得的電流曲線; 插圖表示漏電流 J sh 隨電壓的變化Fig.5.Decompositions of the total current density of perovskite solar cells according to Eq.(1): (a) Devices based on the control MAPbI 3 films; (b) devices based on the DTS passivated MAPbI 3 films; (c) devices based on the DR3T passivated MAPbI 3 films.Jtheoretical represents the total theoretical current, Jbulk represents the bulk recombination current, Jsurf represents the surface recombination current and Jexperimental represents the experimental current.The insets show the bias voltage dependence of J sh.

    圖6 不同界面工程處理下鈣鈦礦太陽電池的效率損失示意圖Fig.6.Efficiency loss of perovskite solar cells with different grain boundaries.

    根據(jù)圖5 和圖6, 未進(jìn)行晶界修飾的鈣鈦礦太陽電池中, 體復(fù)合為主導(dǎo)非輻射復(fù)合機(jī)制, 造成了83.3%的效率損失.引入DTS 后, 體復(fù)合損失減小為48.06%, 而表面復(fù)合損失增大為14.40% (說明判斷不同情況下鈣鈦礦太陽電池的非輻射復(fù)合特性時, 單看非輻射復(fù)合系數(shù)的大小是不準(zhǔn)確的, 需要進(jìn)行效率損失計(jì)算).而在引入DR3T 后, 表面復(fù)合進(jìn)一步增大為38.03%, 體復(fù)合進(jìn)一步減小為19.65%, 表面復(fù)合因此成為該鈣鈦礦太陽電池的主導(dǎo)非輻射復(fù)合機(jī)制.結(jié)合文獻(xiàn)[27] 中的分析可知, DTS 或DR3T 分子與鈣鈦礦之間的相互作用是引起缺陷鈍化、促進(jìn)載流子在晶界間傳輸?shù)脑?因此在鈣鈦礦層引入DTS 和DR3T 能有效地減小鈣鈦礦層的體復(fù)合損失.同時, 由于DTS 對MAPbI3有更大的能極差, 故DTS 在一定程度上會阻礙載流子的傳輸, 即增大串聯(lián)電阻損失.另外,引入DR3T 會使鈣鈦礦層中缺陷態(tài)變得更淺, 在有利于電子和空穴的收集的同時, 也更易產(chǎn)生較大的表面復(fù)合損失.綜上可知, 本等效電路模型的擬合結(jié)果與文獻(xiàn)的分析基本一致.因此, 本文提出的等效電路模型可以準(zhǔn)確判斷實(shí)際鈣鈦礦太陽電池中的主導(dǎo)非輻射復(fù)合機(jī)制, 并且通過各子電流的對比和效率損失的計(jì)算, 可以較好地理解電池的工作機(jī)理, 分析電池優(yōu)化策略的實(shí)用性.

    5 結(jié) 論

    本文提出了一種描述鈣鈦礦太陽電池伏安特性的改進(jìn)等效電路模型, 該模型考慮了光子循環(huán)、陷光結(jié)構(gòu)、非輻射復(fù)合途徑即體復(fù)合和表面復(fù)合、以及串并聯(lián)電阻損耗, 是一個較全面的仿真模型.且通過擬合結(jié)合體復(fù)合和表面復(fù)合公式的漂移-擴(kuò)散模型、仿真得到的J-V曲線、以及不同晶界鈍化處理下的實(shí)際J-V曲線, 驗(yàn)證了該等效電路模型描述鈣鈦礦太陽電池伏安特性的準(zhǔn)確性, 其擬合誤差在2%以內(nèi).基于該等效電路模型, 通過擬合J-V曲線, 可以提取出描述體復(fù)合的體復(fù)合系數(shù)、描述表面復(fù)合的表面復(fù)合系數(shù)、描述歐姆損耗的串聯(lián)電阻和并聯(lián)電阻, 即量化影響電池效率的損耗因素.除此之外, 還可繪制體復(fù)合伏安曲線、表面復(fù)合伏安曲線和電阻伏安曲線, 研究電壓掃描過程中各損耗機(jī)理對效率的影響, 有助于判斷電池內(nèi)的主導(dǎo)損失機(jī)制和弄清工作機(jī)理, 以準(zhǔn)確定位效率提升的重點(diǎn).

    為了逼近Shockley-Queisser 理論極限, 通過對只含體復(fù)合的電池J-V曲線、表面復(fù)合主導(dǎo)的J-V曲線以及不考慮非輻射復(fù)合但改變傳輸層的J-V曲線的仿真和分析, 發(fā)現(xiàn)選用合適的載流子阻擋層可以有效減小串聯(lián)電阻損耗, 甚至將串聯(lián)電阻對效率的影響忽略不計(jì).同時, 優(yōu)化光學(xué)設(shè)計(jì), 提高鈣鈦礦活性層的質(zhì)量及鈍化界面缺陷, 可以有效減小非輻射復(fù)合損耗和并聯(lián)電阻損耗, 獲得更高的光電轉(zhuǎn)換效率和填充因子[30].此外, 本模型在擬合參數(shù)的唯一性上需要進(jìn)一步地改善, 可考慮同時擬合不同光強(qiáng)下的伏安曲線.并且本模型沒有考慮影響鈣鈦礦太陽電池的反常遲滯現(xiàn)象, 因而無法解釋掃描方向不同時離子移動對J-V曲線的影響[31],同樣需要在將來的工作中完善.

    附錄A

    已知體SRH 復(fù)合率在n≈p,n?nt和p?pt情況下的計(jì)算公式為[32]

    其中n為電子濃度,p為空穴濃度,ni為鈣鈦礦活性層本征載流子濃度,pt為陷阱空穴濃度,nt為陷阱電子濃度,τn為電子的體復(fù)合壽命,τp為空穴的體復(fù)合壽命.

    非平衡載流子濃度乘積為

    其中,EFn為電子的準(zhǔn)費(fèi)米能級,EFp為空穴的準(zhǔn)費(fèi)米能級.如果載流子遷移率足夠大, (EFn?EFp)/q為鈣鈦礦體兩端施加的電壓Vap.

    因此體復(fù)合電流應(yīng)為

    其中q為基本電荷,L為鈣鈦礦層厚度,γbulk為體復(fù)合系數(shù),kB為玻爾茲曼常數(shù),T為溫度.

    對于空穴傳輸層/鈣鈦礦層界面, 已知表面SRH 復(fù)合

    率在p??n+,n+?nt和p??pt情況下的計(jì)算公式為

    其中n+為界面處靠近鈣鈦礦側(cè)的非平衡少子電子濃度,p?為界面處靠近傳輸層側(cè)的非平衡多子空穴濃度,τsurfn為電子的表面復(fù)合壽命,τsurfp為空穴的表面復(fù)合壽命.

    根據(jù)玻爾茲曼統(tǒng)計(jì), 少子濃度n+為

    同理, 在電子傳輸層/鈣鈦礦界面處, 靠近鈣鈦礦側(cè)的非平衡少子空穴濃度為

    其中為電子傳輸層/鈣鈦礦層界面處、靠近鈣鈦礦側(cè)的平衡多子電子濃度和平衡少子空穴濃度.

    故空穴傳輸層/鈣鈦礦界面處的表面復(fù)合電流和電子傳輸層/鈣鈦礦界面處的表面復(fù)合電流分別為

    其中Lsurfh為空穴傳輸層/鈣鈦礦層的界面層厚度,Lsurfe為電子傳輸層/鈣鈦礦層的界面層厚度.

    假設(shè)傳輸層/鈣鈦礦層界面的有效厚度為Lsurf,γsurfn≈γsurfp及故有

    附錄B

    基于理想因子法判斷主導(dǎo)非輻射復(fù)合機(jī)制時, 為了避免寄生電阻的影響, Tvingstedt 和Deibel[11]選用開路電壓Voc作為光強(qiáng)的函數(shù), 以更為準(zhǔn)確地求解理想因子.同樣, 考慮到在開路電壓點(diǎn), 實(shí)際鈣鈦礦太陽電池內(nèi)部無電流的情況, 可單獨(dú)考慮非輻射體復(fù)合、非輻射表面復(fù)合、串聯(lián)電阻和并聯(lián)電阻對該點(diǎn)功率大小的影響, 并與該點(diǎn)處的理想功率值比較.得到四組功率和一組理想功率后, 即可根據(jù)功率的損失比例, 量化以上四種因素對鈣鈦礦太陽電池的效率損失.

    以第4 節(jié)實(shí)驗(yàn)結(jié)果分析中的DR3T 情況為例.已知等效電路模型擬合得到的開路電壓為Voc=1.12 V.首先令體復(fù)合系數(shù)γbulk=0 s?1, 表面復(fù)合系數(shù)γsurf=0 s?1, 串聯(lián)電阻Rs= 0 ?·cm2, 并聯(lián)電阻Rsh=Inf ?·cm2, 作出J-V曲線如圖B1中黑色曲線所示.并計(jì)算當(dāng)V=1.12 V 時, 該鈣鈦礦太陽電池對應(yīng)的理想功率Pid=18.18 mW/cm2.第二步, 令γbulk= 7.17×105s?1,γsurf=0 s?1,Rs=0 ?·cm2,Rsh=Inf ?·cm2, 作J-V曲線如圖B1 中紅色曲線所示, 并計(jì)算當(dāng)V= 1.12 V 時, 該鈣鈦礦太陽電池的對應(yīng)功率Pbulk=12.36 mW/cm2.第三步, 令γbulk=0 s?1,Usurf=1.96×10?6nm·cm3·s?1,Rs= 0 ?·cm2,Rsh=Inf ?·cm2, 作J-V曲線如圖B1 中藍(lán)色曲線所示, 并計(jì)算當(dāng)V=1.12 V時, 該鈣鈦礦太陽電池的對應(yīng)功率為Psurf=6.93 mW/cm2.第四步, 令γbulk=0 s?1,Usurf=0 nm·cm3·s?1,Rs=4.20?·cm2,Rsh=Inf ?·cm2, 作J-V曲線如圖B1 中粉色曲線所示, 并計(jì)算當(dāng)V=1.12 V 時, 該鈣鈦礦太陽電池的對應(yīng)功率為PRs=6.42 m W/cm2.最后, 令γbulk=0 s?1,Usurf=0 nm·cm3·s?1,Rs=0 ?·cm2,Rsh=1.63×103?·cm2,作J-V曲線如圖B1 中青色曲線所示, 并計(jì)算當(dāng)V=1.12 V時, 該鈣鈦礦太陽電池的對應(yīng)功率PRsh=17.41 mW/cm2.此時根據(jù)以上各功率與理想功率的關(guān)系, 分別量化體復(fù)合系數(shù)、表面復(fù)合系數(shù)、串聯(lián)電阻以及并聯(lián)電阻對效率損失的影響.

    圖 B1 量化鈣鈦礦太陽電池效率損失的方法示意圖Fig.B1.The method of quantifying efficiency loss of perovskite solar cells.

    猜你喜歡
    傳輸層等效電路太陽電池
    磁致伸縮振動能量收集器的全耦合非線性等效電路模型
    基于撕裂法的變壓器繞組集總參數(shù)等效電路頻率響應(yīng)計(jì)算方法
    基于Python語言的網(wǎng)絡(luò)傳輸層UDP協(xié)議攻擊性行為研究
    ZnO電子傳輸層在有機(jī)無機(jī)雜化鈣鈦礦太陽能電池中的應(yīng)用
    物聯(lián)網(wǎng)無線通信傳輸層動態(tài)通道保障機(jī)制
    基于物聯(lián)網(wǎng)GIS的消防智能巡檢系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)
    幾種新型鈣鈦礦太陽電池的概述
    光對聚合物太陽電池的影響
    柔性砷化鎵太陽電池
    CIGS薄膜太陽電池柔性化
    天天躁狠狠躁夜夜躁狠狠躁| 啦啦啦在线观看免费高清www| 18在线观看网站| 欧美97在线视频| 国产精品国产三级国产专区5o| 少妇人妻精品综合一区二区| 9热在线视频观看99| 老司机影院成人| 亚洲自偷自拍图片 自拍| 成人漫画全彩无遮挡| 日本猛色少妇xxxxx猛交久久| 美女扒开内裤让男人捅视频| 国产亚洲午夜精品一区二区久久| av又黄又爽大尺度在线免费看| 国产极品粉嫩免费观看在线| 69精品国产乱码久久久| 免费在线观看黄色视频的| av片东京热男人的天堂| 大香蕉久久成人网| 美女主播在线视频| 亚洲美女视频黄频| 中国三级夫妇交换| 日韩大片免费观看网站| 黑人猛操日本美女一级片| 成人18禁高潮啪啪吃奶动态图| 综合色丁香网| 亚洲欧美一区二区三区久久| 久久久久网色| 成年av动漫网址| 一级爰片在线观看| 国产精品麻豆人妻色哟哟久久| 男女午夜视频在线观看| 无遮挡黄片免费观看| 美女中出高潮动态图| 国产探花极品一区二区| 你懂的网址亚洲精品在线观看| 国产黄频视频在线观看| 成人国产麻豆网| 久久影院123| 午夜福利视频精品| 精品第一国产精品| 一区二区三区乱码不卡18| 热re99久久精品国产66热6| 亚洲伊人色综图| 少妇被粗大的猛进出69影院| 18禁国产床啪视频网站| 免费观看人在逋| 精品福利永久在线观看| 婷婷色综合大香蕉| 男女边吃奶边做爰视频| av国产精品久久久久影院| 欧美乱码精品一区二区三区| 夜夜骑夜夜射夜夜干| 国产精品久久久久久久久免| 亚洲精品国产区一区二| 亚洲第一av免费看| h视频一区二区三区| 热99久久久久精品小说推荐| 久久久久久久久久久免费av| 中文天堂在线官网| 哪个播放器可以免费观看大片| 日本午夜av视频| 嫩草影视91久久| 免费在线观看黄色视频的| 欧美日韩视频精品一区| 日韩人妻精品一区2区三区| 天天添夜夜摸| 久久久久久久久久久免费av| 最近最新中文字幕大全免费视频 | 熟女少妇亚洲综合色aaa.| 午夜激情久久久久久久| 亚洲精品,欧美精品| 日本欧美视频一区| 赤兔流量卡办理| 国产在线免费精品| 欧美日韩成人在线一区二区| 美女主播在线视频| 操出白浆在线播放| 女性生殖器流出的白浆| 国产精品麻豆人妻色哟哟久久| 精品人妻熟女毛片av久久网站| 欧美变态另类bdsm刘玥| 你懂的网址亚洲精品在线观看| 欧美亚洲 丝袜 人妻 在线| 免费观看av网站的网址| 日韩大码丰满熟妇| 免费观看a级毛片全部| bbb黄色大片| 精品一区二区三区av网在线观看 | 欧美 亚洲 国产 日韩一| 女的被弄到高潮叫床怎么办| 免费少妇av软件| 国产精品二区激情视频| 国产一区二区三区综合在线观看| 欧美日韩一区二区视频在线观看视频在线| 久久精品久久久久久噜噜老黄| 我的亚洲天堂| 另类精品久久| 国产av一区二区精品久久| 老司机影院毛片| av.在线天堂| 欧美日韩成人在线一区二区| 美女扒开内裤让男人捅视频| 热99国产精品久久久久久7| 亚洲美女搞黄在线观看| 国产一卡二卡三卡精品 | 国产97色在线日韩免费| 久久精品国产综合久久久| 久久久久久久国产电影| 两个人免费观看高清视频| 国产成人精品在线电影| 成人毛片60女人毛片免费| 波多野结衣一区麻豆| 最新在线观看一区二区三区 | 大码成人一级视频| 久久久久久久精品精品| 欧美精品亚洲一区二区| 亚洲国产看品久久| 国产伦人伦偷精品视频| 亚洲伊人久久精品综合| 天美传媒精品一区二区| 麻豆精品久久久久久蜜桃| 国产淫语在线视频| 国产熟女欧美一区二区| 亚洲美女黄色视频免费看| 夫妻午夜视频| 一级a爱视频在线免费观看| 日韩精品免费视频一区二区三区| 免费观看av网站的网址| 亚洲国产av影院在线观看| 丝袜美足系列| 久久久久久久国产电影| 亚洲精品美女久久av网站| 免费在线观看黄色视频的| 黄色毛片三级朝国网站| 美女福利国产在线| 美女视频免费永久观看网站| 日韩欧美一区视频在线观看| 90打野战视频偷拍视频| 国产女主播在线喷水免费视频网站| 成人国产av品久久久| 国产片内射在线| 丰满乱子伦码专区| 日韩欧美精品免费久久| 午夜精品国产一区二区电影| 中文字幕人妻丝袜一区二区 | 高清欧美精品videossex| 国产精品久久久久久精品古装| 国产野战对白在线观看| 日韩熟女老妇一区二区性免费视频| 香蕉丝袜av| 别揉我奶头~嗯~啊~动态视频 | 男女之事视频高清在线观看 | 男的添女的下面高潮视频| 国产97色在线日韩免费| 亚洲成色77777| 国产极品天堂在线| 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 99re6热这里在线精品视频| 汤姆久久久久久久影院中文字幕| 亚洲一码二码三码区别大吗| 国产97色在线日韩免费| 久久精品亚洲熟妇少妇任你| 国产又爽黄色视频| 欧美日韩一区二区视频在线观看视频在线| 欧美人与性动交α欧美精品济南到| 日本猛色少妇xxxxx猛交久久| 在线 av 中文字幕| 精品人妻一区二区三区麻豆| 亚洲成人av在线免费| 日本午夜av视频| 9191精品国产免费久久| 亚洲少妇的诱惑av| 纵有疾风起免费观看全集完整版| 男女边吃奶边做爰视频| 国产日韩欧美视频二区| 妹子高潮喷水视频| 啦啦啦 在线观看视频| 亚洲一卡2卡3卡4卡5卡精品中文| 男男h啪啪无遮挡| 国产乱人偷精品视频| 国产成人精品在线电影| 国产成人免费观看mmmm| 不卡av一区二区三区| 国产一区有黄有色的免费视频| 黄频高清免费视频| √禁漫天堂资源中文www| 秋霞在线观看毛片| 午夜福利乱码中文字幕| 亚洲精品视频女| 久久久久人妻精品一区果冻| 亚洲熟女毛片儿| 久久久久精品人妻al黑| 制服诱惑二区| 亚洲精品av麻豆狂野| 韩国精品一区二区三区| 亚洲成人av在线免费| 成人国产av品久久久| 丝袜美足系列| 午夜精品国产一区二区电影| 成人免费观看视频高清| 亚洲欧美成人综合另类久久久| 精品一区二区三区av网在线观看 | 国产老妇伦熟女老妇高清| 久久精品熟女亚洲av麻豆精品| 交换朋友夫妻互换小说| 日韩av免费高清视频| 在线观看免费高清a一片| 久久人人97超碰香蕉20202| 欧美激情极品国产一区二区三区| 免费在线观看黄色视频的| 亚洲精品av麻豆狂野| 纯流量卡能插随身wifi吗| 国产成人精品在线电影| 精品久久蜜臀av无| 欧美乱码精品一区二区三区| 在线观看一区二区三区激情| 老鸭窝网址在线观看| 精品一区在线观看国产| 久久精品久久久久久久性| 黄色视频不卡| 国产成人免费观看mmmm| 日韩免费高清中文字幕av| 丝袜喷水一区| 天堂中文最新版在线下载| 国产精品秋霞免费鲁丝片| 热re99久久国产66热| 国产深夜福利视频在线观看| 色网站视频免费| 免费黄频网站在线观看国产| av有码第一页| 国产欧美日韩综合在线一区二区| 国产精品欧美亚洲77777| 国产精品嫩草影院av在线观看| 亚洲国产看品久久| 亚洲av日韩精品久久久久久密 | 国产深夜福利视频在线观看| 80岁老熟妇乱子伦牲交| 街头女战士在线观看网站| 日日爽夜夜爽网站| 亚洲国产欧美日韩在线播放| 热99久久久久精品小说推荐| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产精品一区二区三区在线| 久久久国产一区二区| 精品少妇久久久久久888优播| 一区二区三区激情视频| 熟妇人妻不卡中文字幕| 99国产综合亚洲精品| 欧美日韩成人在线一区二区| 青春草亚洲视频在线观看| 一边摸一边做爽爽视频免费| 女人高潮潮喷娇喘18禁视频| 免费高清在线观看日韩| 少妇被粗大猛烈的视频| av卡一久久| 97在线人人人人妻| 国产在线视频一区二区| 日韩制服骚丝袜av| 超碰成人久久| 狠狠精品人妻久久久久久综合| av视频免费观看在线观看| a级毛片在线看网站| 久久久精品94久久精品| 亚洲精品第二区| 色94色欧美一区二区| 国产成人精品无人区| avwww免费| tube8黄色片| 亚洲精品久久久久久婷婷小说| 99久久综合免费| 国产精品久久久av美女十八| 久久久久久久久久久免费av| 赤兔流量卡办理| 日韩伦理黄色片| 国产成人免费观看mmmm| 国产精品一二三区在线看| 欧美最新免费一区二区三区| 欧美少妇被猛烈插入视频| 国产欧美日韩一区二区三区在线| 中文天堂在线官网| 一边摸一边抽搐一进一出视频| 最近的中文字幕免费完整| 在线天堂中文资源库| 亚洲一码二码三码区别大吗| 自拍欧美九色日韩亚洲蝌蚪91| 美国免费a级毛片| 亚洲av成人不卡在线观看播放网 | 老司机在亚洲福利影院| 国产毛片在线视频| 国产极品粉嫩免费观看在线| 看免费av毛片| 日本91视频免费播放| 成年人免费黄色播放视频| 欧美日韩一级在线毛片| 亚洲 欧美一区二区三区| 黄网站色视频无遮挡免费观看| 十八禁网站网址无遮挡| 视频在线观看一区二区三区| 好男人视频免费观看在线| 亚洲成人国产一区在线观看 | 午夜免费鲁丝| 18禁观看日本| 男女免费视频国产| 多毛熟女@视频| 免费日韩欧美在线观看| 成年动漫av网址| 亚洲精品视频女| 校园人妻丝袜中文字幕| 久久久久久久国产电影| 亚洲欧美精品自产自拍| 黑人猛操日本美女一级片| 另类精品久久| 国产伦人伦偷精品视频| 夫妻午夜视频| 久久久久久免费高清国产稀缺| 亚洲精品美女久久久久99蜜臀 | 各种免费的搞黄视频| 亚洲欧美一区二区三区国产| 久久鲁丝午夜福利片| 久久久久国产一级毛片高清牌| 国产精品一区二区精品视频观看| 国产午夜精品一二区理论片| 日日撸夜夜添| 国产日韩欧美视频二区| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲美女视频黄频| 一级毛片我不卡| 日韩一区二区视频免费看| 黄片无遮挡物在线观看| 色视频在线一区二区三区| 嫩草影视91久久| 男女之事视频高清在线观看 | 成人国产av品久久久| kizo精华| 国产一区二区三区av在线| 国产免费视频播放在线视频| 日韩一区二区视频免费看| 国产国语露脸激情在线看| 老司机在亚洲福利影院| 国产av精品麻豆| 国产精品一区二区在线不卡| 飞空精品影院首页| 交换朋友夫妻互换小说| 日韩欧美精品免费久久| 国产一卡二卡三卡精品 | 国产成人av激情在线播放| 国产精品麻豆人妻色哟哟久久| 在线观看免费日韩欧美大片| 亚洲国产欧美网| 国产成人系列免费观看| 超色免费av| 亚洲国产欧美在线一区| 美女主播在线视频| 久久久久精品性色| 看免费成人av毛片| 国产一区亚洲一区在线观看| 欧美 亚洲 国产 日韩一| 久久久精品免费免费高清| 视频在线观看一区二区三区| 久久青草综合色| 日本欧美视频一区| 在线看a的网站| 国产一区二区 视频在线| 熟妇人妻不卡中文字幕| 纵有疾风起免费观看全集完整版| 我的亚洲天堂| 熟妇人妻不卡中文字幕| 男女免费视频国产| 天堂8中文在线网| 日韩视频在线欧美| 丝袜人妻中文字幕| 午夜激情av网站| 亚洲国产欧美一区二区综合| 一边摸一边做爽爽视频免费| 黄片无遮挡物在线观看| 女的被弄到高潮叫床怎么办| 精品国产一区二区久久| 国产乱人偷精品视频| 99久久精品国产亚洲精品| 十八禁高潮呻吟视频| 操美女的视频在线观看| 日本欧美视频一区| 欧美另类一区| 欧美老熟妇乱子伦牲交| 国产99久久九九免费精品| 亚洲国产日韩一区二区| 久久韩国三级中文字幕| 欧美日韩亚洲综合一区二区三区_| 最近手机中文字幕大全| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 亚洲人成电影观看| 亚洲国产欧美日韩在线播放| 国产麻豆69| 亚洲图色成人| 中国三级夫妇交换| 欧美在线黄色| 久久ye,这里只有精品| 欧美变态另类bdsm刘玥| 综合色丁香网| videosex国产| 亚洲欧美清纯卡通| 满18在线观看网站| 99热全是精品| 99精品久久久久人妻精品| 久久久久久人妻| 99久久精品国产亚洲精品| 久久久国产欧美日韩av| 一级片'在线观看视频| 一区福利在线观看| 9色porny在线观看| 国产精品.久久久| av卡一久久| 美女脱内裤让男人舔精品视频| 超碰成人久久| 亚洲精品日韩在线中文字幕| 超碰97精品在线观看| 日本午夜av视频| www.自偷自拍.com| 中文字幕色久视频| 久久精品aⅴ一区二区三区四区| 多毛熟女@视频| 国产伦理片在线播放av一区| 高清不卡的av网站| 肉色欧美久久久久久久蜜桃| 一二三四在线观看免费中文在| 大码成人一级视频| 80岁老熟妇乱子伦牲交| 日本av手机在线免费观看| 捣出白浆h1v1| 女人爽到高潮嗷嗷叫在线视频| 天堂中文最新版在线下载| 丝袜脚勾引网站| 麻豆av在线久日| 亚洲精品中文字幕在线视频| 少妇人妻久久综合中文| 日韩精品免费视频一区二区三区| 国产精品女同一区二区软件| 午夜福利免费观看在线| www日本在线高清视频| 天天躁夜夜躁狠狠久久av| 亚洲欧美激情在线| 日韩精品有码人妻一区| 搡老岳熟女国产| 久久久精品94久久精品| 成年人免费黄色播放视频| 久久精品国产综合久久久| 国产成人精品久久久久久| 亚洲伊人久久精品综合| 啦啦啦啦在线视频资源| 亚洲中文av在线| 中文天堂在线官网| 成人漫画全彩无遮挡| 国产精品香港三级国产av潘金莲 | 久久狼人影院| 如日韩欧美国产精品一区二区三区| 一本大道久久a久久精品| 岛国毛片在线播放| 亚洲情色 制服丝袜| 人妻 亚洲 视频| 欧美日韩亚洲国产一区二区在线观看 | 最近中文字幕2019免费版| 国产精品久久久久成人av| 人成视频在线观看免费观看| 国产精品一二三区在线看| 久久精品久久久久久噜噜老黄| 久久久国产欧美日韩av| www.精华液| 99国产精品免费福利视频| 卡戴珊不雅视频在线播放| 男女之事视频高清在线观看 | 久久人人97超碰香蕉20202| 一本—道久久a久久精品蜜桃钙片| 美女午夜性视频免费| 精品一区在线观看国产| 五月开心婷婷网| 欧美国产精品va在线观看不卡| 18在线观看网站| 18禁观看日本| 大片电影免费在线观看免费| 99国产精品免费福利视频| 亚洲国产最新在线播放| 美女大奶头黄色视频| 老汉色∧v一级毛片| 桃花免费在线播放| 欧美变态另类bdsm刘玥| 久久国产亚洲av麻豆专区| 精品免费久久久久久久清纯 | 欧美激情 高清一区二区三区| 国产av国产精品国产| 免费观看av网站的网址| 国产精品国产三级国产专区5o| 人人妻人人澡人人看| 亚洲欧美色中文字幕在线| 成人影院久久| 久久国产精品大桥未久av| 黄色视频在线播放观看不卡| 少妇被粗大猛烈的视频| 男人添女人高潮全过程视频| 久久99热这里只频精品6学生| 欧美激情高清一区二区三区 | 国产麻豆69| 日韩一区二区视频免费看| 我的亚洲天堂| 天堂中文最新版在线下载| 国产乱来视频区| 国产麻豆69| 一区二区日韩欧美中文字幕| 免费观看性生交大片5| 天堂中文最新版在线下载| 天堂俺去俺来也www色官网| 国产成人91sexporn| 精品人妻在线不人妻| 免费观看性生交大片5| 亚洲欧美一区二区三区黑人| 婷婷色av中文字幕| 国产精品亚洲av一区麻豆 | 大片免费播放器 马上看| 纵有疾风起免费观看全集完整版| 18禁动态无遮挡网站| 爱豆传媒免费全集在线观看| 一本一本久久a久久精品综合妖精| 丰满少妇做爰视频| 欧美亚洲 丝袜 人妻 在线| 国产乱人偷精品视频| 国产亚洲最大av| 欧美最新免费一区二区三区| 亚洲五月色婷婷综合| 亚洲欧洲国产日韩| 国产av一区二区精品久久| 啦啦啦视频在线资源免费观看| 亚洲欧美一区二区三区黑人| 亚洲国产欧美日韩在线播放| 欧美日韩亚洲国产一区二区在线观看 | 精品一区二区三卡| 日韩成人av中文字幕在线观看| 成年美女黄网站色视频大全免费| 国产精品久久久久久精品古装| 日韩中文字幕视频在线看片| 亚洲成av片中文字幕在线观看| 欧美 亚洲 国产 日韩一| 制服诱惑二区| 欧美97在线视频| av在线播放精品| 欧美精品一区二区大全| 免费黄网站久久成人精品| 欧美中文综合在线视频| 午夜福利免费观看在线| 一级黄片播放器| 性少妇av在线| 人妻一区二区av| 亚洲四区av| 午夜av观看不卡| 欧美精品高潮呻吟av久久| 精品少妇久久久久久888优播| 满18在线观看网站| 久久久久久久国产电影| www.熟女人妻精品国产| 女性生殖器流出的白浆| 亚洲精品成人av观看孕妇| 色婷婷av一区二区三区视频| 国产毛片在线视频| 老司机亚洲免费影院| 精品国产乱码久久久久久小说| 久久久精品免费免费高清| 99九九在线精品视频| 久久性视频一级片| 丝袜喷水一区| av国产精品久久久久影院| 老鸭窝网址在线观看| 亚洲美女搞黄在线观看| 国产精品三级大全| 国产高清不卡午夜福利| 欧美成人午夜精品| 久久久久视频综合| 老汉色av国产亚洲站长工具| 午夜激情av网站| 国产成人91sexporn| 国产毛片在线视频| 中文乱码字字幕精品一区二区三区| 国产极品粉嫩免费观看在线| 蜜桃在线观看..| 一本大道久久a久久精品| 美女主播在线视频| 少妇被粗大的猛进出69影院| 亚洲综合精品二区| 国产一区二区三区av在线| 日韩电影二区| 色播在线永久视频| 中文欧美无线码| 国产xxxxx性猛交| 男男h啪啪无遮挡| 久久99精品国语久久久| 欧美在线黄色| 色播在线永久视频| 国产精品无大码| 一区二区日韩欧美中文字幕| 一级毛片我不卡| 欧美人与性动交α欧美精品济南到| 日韩一区二区视频免费看| 一级毛片我不卡| 99香蕉大伊视频| 亚洲精华国产精华液的使用体验| 超碰成人久久| 亚洲欧美成人综合另类久久久| 久久毛片免费看一区二区三区| 亚洲自偷自拍图片 自拍| 观看av在线不卡| 赤兔流量卡办理| 日本欧美国产在线视频| 久久久精品国产亚洲av高清涩受| 男女床上黄色一级片免费看|