亚洲人天堂I色综合久久久久网I超碰精品在线观看I激情伊人I欧美天堂久久I九九热在线观看I一区二区精品在线I日韩美女黄色片

技術文章

Technical articles

當前位置:首頁技術文章理解iR補償 ——溶液電阻補償

理解iR補償 ——溶液電阻補償

更新時間:2020-03-26點擊次數:10913

引言

在Gamry聽到的一些常見的技術問題都與iR補償有關–

  • 未補償的iR是哪來的?
  • 我需要在我的實驗中用iR補償嗎?
  • 我應該如何設置iR補償參數?
  • 為什么iR補償在我的系統中不能運行呢?

在這個技術報告中,我們將嘗試回答這些問題,讓您對iR補償有一個基本的了解。

背景

一些初級背景信息描述了iR誤差大體的情況。后面的內容主要集中于介紹用于測量和修正iR誤差的“電流截斷”iR補償方法。正反饋iR補償有提到,但是沒有具體討論。

該技術報告假定您對電化學工作站操作有一定了解。如果沒有,您可以瀏覽我們的Primer on Potentiostats。有經驗的電化學工作站使用者應該跳過初級教程,繼續閱讀本文。

如果您對電化學阻抗的基本原理有一些了解的話,也有幫助。在我們的網頁上有Primer on Electrochemical Impedance可以查閱。請特別注意一下,典型化學過程是如何被映射到電路元件中的。

什么時候我需要iR補償?

關于這個問題,我們將在此給出一個大概的回答。更加完整的答案需要知道下面討論的信息和一些關于您測試系統的信息才能得出。

一般而言,當如下列舉中的一個或多個條件實現時,我們就需要iR補償:

  • 您正在做一個終產生數值結果的量化測試,例如腐蝕速率,平衡常數或速率常數

  • 您電解池中的溶液不導電

  • 您使用的電流非常得高

  • 您的電解池幾何不夠理想

遺憾的是,這些標準都是主觀的。舉個例子,0.5 M KCl 水溶液在電化學分析應用中可能被認為導電性很好,而電鍍應用中導電性卻很差。

一個常用的簡單的經驗法則:

  • 記錄一下有iR補償和沒有iR補償的初始數據曲線

如果當iR補償運行后曲線的形狀明顯得變化了,那么補償是必需的。

iR補償常會給數據增加額外噪聲,所以曲線上增加的噪聲不被考慮為一個明顯的變化。

iR誤差哪來的?

讓我們來看看一個典型的三電極電化學測試電解池。我們在電解池里標注了一些參考點,這個電解池將在該技術報告的剩余部分多次引用。

Counter electrode output    at the potentiostat

參考點:

A: 對電極在電化學工作站處的輸出端

B: 對電極的金屬表面

C: 對電極的電解質表面

D: 魯金毛細管的電解質

E: 工作電極的電解質表面

F: 工作電極的金屬表面

G: 參比電極的電解質表面

H: 參比電極在電化學工作站處的輸入端

I: 工作電極在電化學工作站處的輸出端

我們也可以把電解池理解為一個(簡化的)電子元件網絡,如下所示:

a network of electronic components

幾乎所有的電化學工作站都能夠較好地控制和測量H和I點間的電勢。

遺憾的是,我們真的想要控制和測量的其實是E和F點間的電勢。這是我們想要知道的通過電化學界面的電勢差。

Vimportant = Vf — Ve

為了討論這個問題,I點與F點等值。G等于H加上由工作電極/參比電極電勢差引起的一個常數補償電壓,也叫做開路電壓Voc。因為參比中沒有電流通過,Rbrigde上的電壓降為0,所以G與D等值。從

Vmeasured = VI — VH

我們可以得到

Vmeasured = Vf — Vd + Voc

我們離得到反應電勢僅有一步之遙了。E點等于D加上討厭的電阻Ru

這怎么關聯起來呢?通過歐姆定律:

Ve — Vd = Icell * Ru

所以,把上式代入后得到:

Vmeasured = Vf — Ve + Icell * Ru + Voc

重新排列結果后得到:

Vimportant = Vf — Ve = Vmeasured — Icell * Ru — Voc

我們可以測量電壓Vmeasured。我們測量電流Icell。我們能查閱或測量開路電壓Voc。不過不知道Ru,我們不能得到Vimportant

而這個問題是iR補償的核心。iR代表Icell﹡Ru

幸運地是,Ru可以被測量和修正。

什么類型的電解池/系統我們需要擔心呢?

多數情形,電化學家們能夠操縱他們的實驗,所以iR降不是一個問題。一個簡單的方法就是加入一些不會反應的鹽,酸或者能增加電解質導電性的基底(支撐電解液)。

導電性增加 –--> Ru降低 –> iR誤差降低

另一個方法是去減小魯金毛細管和工作電極間的距離。可以設計這樣的電解池以使距離非常地小。

顯然,如果真這么簡單的話,我們就不會有那么多的問題了。

增加支撐電解質將影響電化學性能,盡管離子不直接參與反應。它將改變雙電層(Cf)的成分。這可能會影響反應物和產物的溶解度或結構。也可能改變表面原子層的結構。在很多的例子中,不改變電解質可能很重要。例如,腐蝕化學家想要研究他或她的電解質的腐蝕性,而不是電解質加鹽時的情形。

同樣的,電解池設計方法常常不是一個好的解決方案。很多電解池設計被物理上局限了。一個更微妙問題是,把參比電極和工作電極放的太近會改變表面的電流密度,導致改變我們正嘗試要做的測量結果。

 No surface effect-current desisty lines have room to adjust to probe

 

Reducted current density effect 

所以,您需要擔心iR誤差嗎?

如果您的電解質不導電或者您的參比探針與反應表面離得非常遠,那么是的。通常您能夠很容易地測量iR。

iR誤差如何測量?

下面的電路中有一些關于測量iR和Ru的重要線索。它代表了一個電化學電解池電學行為的常見簡化模型。

Look at Rfaradaic. It has a capacitor, Cfaradaic in parallel with it. Ru doesn't

看看Rfaradaic,它有一個電容Cfaradaic與其并聯。而Ru沒有。

由此表明,交流實驗可以區分兩種電阻。高頻信號直接通過Cf沒有電壓降,而它們在通過Ru時被迫有Icell﹡Ru的電壓降,就像低頻信號。

的確,如果您記錄了這個電解池的電化學阻抗譜,它的Bode表現看起來就像下圖:

record the Electrochemical Impedance Spectrum of this cell

低頻時Cf近乎為開路,測得的阻抗是Ru和Rf的總和。高頻時Cf相當于短路,測得的阻抗是Ru

所以如果您需要擔心Ru,您可以測量它然后乘以電解池電流。

產生的誤差電壓就是沒有補償的iR,Ve-Vd。如果iR小于幾個毫伏,不需要再擔心它了。

舉個例子,從阻抗譜假設Ru = 100歐姆,假定Icell = 10 μA。

iR = 100﹡10﹡10-6 = 1 mV

對于大多數的電化學現象,這是一個小誤差。

從另一個角度考慮 – 如果Rf>>Ru, 則不用考慮Ru

用直流技術測量Ru

用交流阻抗測量Ru理論上是個很好的方法,但是有些時候我們需要一個更快捷更便宜的方法。而且常常,我們想要一邊測試一邊做些其他事,比如記錄電流vs.電壓曲線。

幸運的是,可以用一個等效直流方法。事實上,對它更好的形容是“瞬態”技術。它被稱作電流截斷iR補償。

同樣,考慮用簡單的Randle電解池模型模擬具有溶液電阻的電化學反應。

Gamry的電化學工作站含有一個電路,它可以快速關閉(切斷)流經電解池的電流,等候一個較短的時間(10微秒到30000微秒),然后再把電流打開。很多其他商業儀器也有相同的功能。

為了做電流切斷測量,電解池電壓(Vmeasured)在電流切斷的緊接前和緊接后測得。理想地,測得的電解池電壓將如下圖所示:

Suppose we are measuring 1.0 V while current is flowing

假定電流流動時,我們在測1.0 V。在時間0時,我們非常快地切斷電流!電壓馬上下降一個Ru的電壓降。然后電壓開始慢慢地降低。

為什么?

電壓的減小是由感應電流電容器(Rfaradaic)慢放電引起的。這個現象時間較長時變得重要。短時間段內,電容器可以保持電壓到Vmeasured – Vu,這是我們的Vimportant。而這是我們想要得到的。

照例,事情不會如此簡單。三個問題妨礙了這個簡單模型:采樣速率,輸出電容和噪聲。

采樣速率

理想化電流截斷波形伴隨的一個問題是采樣的速率。在上圖中,采樣是2微秒(非常快)。如果我們假定衰變曲線是直線然后反向外推到切斷時間,我們可以明顯調慢采樣。

我們用我們一直用的相同Randle電解池試試。采樣在1毫秒和2毫秒,反向外推至0毫秒開關關閉的時候。您可以得到類似下圖的結果:

The estimated Vu is calculated

Vu可以通過Vu=V­1+(V1-V2)估算。

在這個例子中,V1=0.671 V,V2=0.481 V,而通過直線外推,Vu估值在0.862 V。

不過等等,Vu常常非常接近0.938 V。我們被誤導選擇了太慢的時基。您能夠在圖上看到關閉跡線的曲率。當然它是彎曲的-它是指數式衰減。

順便提一下,這些都是從我們原始Randel電解池模型得來的真實數值,用Mathcad計算然后用Excel畫圖。

所以,我們離題太遠了。毫秒的時間設置對于該電解池的iR測量太慢了。應該使用其他更快但又不至于太快的時基。

您如何辨別?

在數學上,衰減時間常數是Rfaradaic﹡Cfaradaic。對于這個電解池,

taufaradaic=3000 Ohm﹡1μF=3 ms

如果您粗略地知道這些值,您可以挑選一個短的時間tau為RC/10。或者您可以通過減小采樣速率直到數值穩定來慢慢地接近正確答案。

不過采樣速率變短存在著另一個問題 – 恒電位輸出和纜線電容。

纜線電容

回想原始電解池模型:

那個看起來無害的電容Ccable能夠引起很多的問題。如果您有典型的屏蔽電纜,Ccable的值可以是每英尺50 pF。對于一個5英尺長的纜線就是250 pF。另外還有一個差不多100 pF大小的給開關本身,電路板和驅動放大器的電容。

我們可以用這個電路作為一個模型:

The cable capacitance forms an RC section with Ru and Rsolution

纜線電容產生一個與Ru和Rsolution在一起的RC部分。這意味著通過Ru的電壓沒有無限快地消失。

為了這個討論,我們必須假設對電極電容很大而且在這些時標充當短路。幸運的是,這個假設很合理。

假定您決定將iR采樣設在50微秒和100微秒。這些點在下圖中顯示為方塊。用這兩種測量,iR估算很明顯又再次非常不精確。您必須等到纜線電容瞬變消失才能進行測量。

helpful to see this on a log(time) scale so you can see the cell cable & Faradaic capacitors being discharged

用時間指數比例查看這一現象可能會有幫助,所以您可以查看電解池電纜和法拉第電容在放電。

a positive noise contribution to V1 and a negative noise contribution to V2.

您需要找到兩個限制放電曲線間的時間范圍。

纜線電容必須充分放電,但是法拉第電容必須仍然在近似線性區域。如果您的法拉第電容不太大,比纜線電容大很多,電流截斷iR補償將會有用。

噪聲

當電流截斷iR補償被用在實際系統上時,噪聲可能是一個主要問題。

基本上,電流截斷iR補償是一個微差測量。記住我們用于估算Vu的方程:

Vu = V1 + (V1-V2)

差分項V1-V2對噪聲非常敏感。

在適于電流截斷iR補償的系統中,V1和V2的差值很小,從幾個毫伏到幾百個毫伏。

假定V1有一個正噪聲貢獻,V2有一個負噪聲貢獻。那么平均噪聲是0,但是Vu的誤差是兩倍那么大!

The control loop algorithm treats the potentiostat as a loop within a loop

您可能會說:“那又怎樣。只要開一些濾波器。噪聲就會沒了。”

但是我們正嘗試測量一個快速(10-1000秒)的現象。這種情況下我們不能增加一個5 Hz的濾波器,如果它沒有*地消失的話,整個瞬變將會失真。

下面列出的所有都能幫助減小噪聲:

  • 使用法拉第籠將外來噪聲隔絕在測量以外。

  • 使用信號平均法以使在保留真實值的同時噪聲項被平均掉。

  • 如果您的噪聲源頻率是已知的,使用同步采樣法以便所有噪聲形成的誤差在同一個方向。

  • 后,如果噪聲還是太大,不要用外推法。就取個平均值吧,比如:


Vu = (V1 + V2) / 2

當您嘗試測量低電流時,噪聲情況更差。在這個情況下,當電流截斷開關打開時,參比電極和工作電極在更高電流時有更多噪聲。

如何修正?

至此,我們只談論了iR的測量。如果我們知道Vu的值,我們可以減去Vmeasured值,得到Vimportant。這被稱為后處理修正。

后處理修正會遇到的一個問題是我們在施加真實電勢前不能預測其數值。這在掃描電位的實驗中是特別有問題的。在這些實驗中,電位掃描速率不是常數,掃描極限可能非常不準確。

我們會傾向于連續地進行恒電位測量和Vu修正。畢竟當您想施加1 V到電解池上,您想要Vimportant = 1 V,而不是Vmeasured = 1 V。

當您使用電化學工作站時,情況變得簡單了。電化學工作站不需要對Vu進行修正,因為它的工作是控制電流而不是電位。它還需要測量Vu

盡管不是有效但有一個簡單的用電流截斷自動修正iR的方法是讓電化學工作站在施加信號上增加其對Vu的合適評估值。這一方法可以用下列公式表示,其中方框中的數值代表測量點:

Vapplied[i] = Vrequested[i] + Vu[i-1]

初值沒有修正。然后從第一個數據點測得的誤差將被疊加到第二個數據點的施加電壓上。繼續這樣做,隨著數據的積累修正變得更加準確。

需要注意的是這是一個動力學修正。Ru可以在實驗中改變,系統將會自動地補償這個改變。

控制回路算法

將誤差電位直接反饋進第二個數據點不是運行修正的復雜的方法。一個更好的了解反饋機制的方法是把iR修正看成一個控制回路。

控制回路算法把電化學工作站當成一個回路中的回路。內部回路是電化學工作站本身,它測量Vmeasured然后用反饋機制控制它。這個回路純粹由模擬電子學制成,如下所示:

The iR correction also happens in a control loop outside of the potentiostat loop

同樣,我們減少了一些與這部分討論不相關的組件。

電化學工作站是一個控制回路。它測量Vmeasured,并與Vapplied作比較,對對電極電壓做修正直到兩者間的電勢差為0。所有這些是連續發生的。

iR修正也常發生在電化學工作站回路以外的控制回路。如下圖所示:

 block that produces Vactual

外回路與內回路做著相同的工作,不過它在電腦上數字化地完成修正。它的工作是查看Vimportant = Vapplied,而現在內回路的工作是查看Vmeasured = Vactual。Vactua是來自外回路的新值。

我們也展示了一個有趣的產生Vactual的模塊。Vactual被稱為增益模塊或控制器模塊,取決于您跟誰說。它的輸出值由下式給出:

Simplified Potentiostat with PFIR D/A Converter

它已知將工程師們控制成PID回路。每一個增益是單獨控制的。通過調整增益,我們能夠得到iR補償,比簡單反饋算法運行得更加好。事實上,我們想要用積分控制做iR補償。

在Gamry,我們通常在我們的直流軟件中使用控制回路算法。它用一些我們實驗控制語言里得到的線來設置。說明,就像下面這樣:

if (IRToggle)Pstat.SetIruptMode (IruptClfg, EuExtrap, IruptTime, POTEN.Eoc (), 0.8)Pstat.SetVchFilter (100000.0)elsePstat.SetIruptMode (IruptOff)Pstat.SetVchFilter (5.0)

0.8是控制回路增益,Ki,Kp定在1.0,Kd是0。注意電壓通道濾波器(VchFilter)被設為通過100 kHz信號。如果我們不用iR補償,濾波器被設為過濾盡可能多的環境噪聲。

我們可以通過修改包括控制回路模式,Vu計算,電流截斷計時和增益等的參數以適應反應和測試中的電解池動力學。

電流截斷iR補償的優勢

電流截斷iR補償與其他iR補償方法相比有一些優勢,包括:

  • 不需要提前知道Ru

  • Ru可以在實驗過程中改變,而不制造補償中的誤差

  • 補償與用于測量電流的電流量程無關,所以在自動量程實驗中運行

  • 掃描參數如斜升極限值和掃描速率都自動修正了

電流截斷iR補償的局限性

電流截斷iR補償在一些電化學系統上運行良好,但是當應用到其他系統時無法正常工作。

這些失敗在技術現實局限性方面一般是可說明的。

局限性包括:

  • 需要一個大的法拉第電容

  • 每點時間局限性

  • Rfaradaic應大于Ru

  • Ru的值必須小于一個極限值

下面將詳細解釋這些局限性。

需要一個大的法拉第電容

如上所述,當電流被斷時Cfaradaic保持“直流”電位。如果法拉第電容缺失或者太小,電流截斷通常會驅動系統到一個大的電位和電流。這個問題明顯的癥狀是測得電流比預期電解池電流要高很多倍。過載指示也可能會出現。

Gamry電化學工作站的電流截斷在法拉第電容大于20 μF時效果好。對于一個“裸金屬”電極,您可以估算電容為20 μF/cm2,所以電極面積必須大于等于1 cm2。如果電極覆蓋有任何形式絕緣涂層的電極時,我們不建議您使用電流截斷iR補償。

這個要求通常將電流截斷iR補償限制在腐蝕測試與電池和燃料電池研究。電流截斷對于常用于物理電化學電解池的電極尺寸效果不太好。

每點時間局限性

電流截斷iR補償假設您多數時間施加的是直流電位和電流。分斷時間應該比測量數據曲線中每一個數據點所需的時間小很多。

默認情況下,Gamry軟件自動地為分斷選擇一個總的電流截斷時間和采樣時間。每當電流量程變化時,這些時間會被調整,更長的分斷時間和更慢的采樣速率用于更敏感的電流量程。分斷時間的范圍通常為10微秒到64微秒。

Gamry建議當每個數據點時間等于或大于1秒時,您僅施加電流截斷iR補償。

如果您正在掃描電位,這將掃描速率限制為等于或小于5 mV/sec。

Rfaradaic應該大于Ru

Rfaradaic和Ru間的比值也有限制。因為相同的電流流經兩個電阻,這也是通過電化學界面的電位和誤差電位間比值的一個限制。

通常關于這個比值的一個更嚴重的局限性出現在電解池。大多數電化學電解池在電極表面擁有一個非均勻的電流分布。工作電極的一些部分有比其他部分更多的電流。在這樣的情況下,簡單的Randle電解池模型不適用。界面不能用單個電位來表述。

除非您的電解池有一個專為均勻電流而設計的幾何,否則我們相信您應該將Ru保持等于或小于Rfaradaic的十分之一。如果比值大于1/10,我們相信您系統上得到的任何量化結果都將有誤差。需要注意的是這只是一個“直覺”近似。我們不能保證這個近似值將適用于您的電化學系統。

Ru的值不太大

Ru的值也有一個基線,其與Rfaradaic的值無關。我們的經驗顯示當Ru超過一些上*誤差會發生。對于大多數的Gamry系統,這個極限大約是10 kΩ。

正反饋iR補償

電流截斷iR補償只有當您在測試慢現象時是有用的,例如腐蝕反應或能源儲存設備表征。當需要非常快的測量時,它不能用。快速實驗的一個例子是用于測量化學熱力學的1000 V/sec循環伏安法。

幸運的是,有一個iR補償技術可以用于快速系統。它被稱為正反饋iR補償。

這個技術可以被認為是電化學工作站的一個額外模擬反饋路徑。所有有用的電化學工作站測量電解池電流。當電化學工作站里的正反饋激活了,一部分的電流信號以額外電壓輸入的方式反饋回來。

下圖是一個電化學工作站的高度簡化示意圖。在這個示意圖的右下部分,電流通過Rm上的電壓降來測量。這個電壓降被放大去產生一個電壓信號,叫做Isig。這個示意圖中,Isig在滿度電流是3伏。

在示意圖的右上部分,我們展示了一個標為PFIR DAIC的模塊。這是正反饋IR數模轉換器。它的輸出是一個電壓,是Isig已知的部分。當正反饋iR被激活,這個電壓以額外電壓輸入的方式被施加到控制放大器上。

注意:在這個示意圖和如下的討論中,PFIR DAC輸出在被用作反饋源之前不按比例縮小。這可能對于全部PFIR的實現不對--包括一些Gamry的實現。

含PFIR D/A轉換器的簡化電化學工作站

the PFIR DAC output is not scaled prior to being used as a source of feedbac

一些簡單的數學推導可獲得:

Isig = 3.0 × Icell/IFS = Icell ×Re

其中,Re是等效電流測量電阻,由下式得:

Re = 3.0 V / IFS

Re告訴我們能夠在任意電流量程上補償的Ru的大值。

在正反饋iR補償中,您需要在修正前知道Ru。Gamry的PHE軟件包含了一個簡便的先于實驗運行Ru測量方法。這個方法使用一個簡化的EIS形式,假設高頻的電解池阻抗等于Ru

一旦您給Ru輸入一個值并選擇正反饋iR修正,Gamry軟件將PFIR DAC設為輸出一部分Isig,Isig等于Ru對Re的比值。在這一設置下,電壓反饋是:

PFIR out = Ru/Re × Isig = Ru/Re × Icell × Re = Ru × Icell

我們通過Ru上的電壓來增加施加在電解池上的電勢。數值的分辨率由PFIR DAC分辨率控制。用一個14位DAC(數模轉換器),分辨率是Re/16384。

舉一個例子,我們來看看3 mA電流量程。在這個量程,Re是1000 Ω。用14位DAC的正反饋修正有一個1000/16834的分辨率或者每比特0.061 Ω。

正反饋iR補償的優勢

正反饋iR補償相比其他的iR補償方法有一些優勢,包括:

  • 非常快速的實驗可用

  • 掃描參數,例如斜升極限和掃描速率,被修正

正反饋iR補償的局限性

相比于其他補償方法,正反饋iR補償有一些局限性:

  • 需要預先知道Ru的值

  • 如果Ru在實驗過程中變化了,會有誤差

  • 在實驗過程中電流流程必須不變

  • 正反饋能導致電化學工作站振蕩

四虎国产成人精品免费一女五男 | 亚洲综合在线成人 | 色综合久久久久无码专区 | 久久伊人草 | 国产午夜精品一区二区三区欧美 | 亚洲黄色在线观看 | 亚洲黄色在线看 | 日韩精品久久久久久久酒店 | 国产精品久久久久久妇女 | 国产精九九网站漫画 | 游戏涩涩免费网站 | 男女互操视频 | 亚州综合视频 | 国产人澡人澡澡澡人碰视频 | 日韩欧美性视频 | 大波大乳videos巨大 | 天天摸天天操 | 人妻天天爽夜夜爽一区二区三区 | 色伊人| 日韩一级片免费在线观看 | 99激情网 | 国产午夜在线播放 | 屁屁影院一区二区三区 | 91插插插插插插插 | 久久久久久99精品久久久 | 国产精品成人国产乱一区 | 潘金莲三级80分钟 | 国产成人精品免高潮在线观看 | 国产精品爽爽爽 | 亚洲一区二区三区婷婷 | 日韩免费一级片 | 高跟鞋和丝袜猛烈xxxxxx | 久久婷婷伊人 | 青草青草久热 | 免费成人黄色网址 | 欧美激情999| 特黄老太婆aa毛毛片 | 亚洲电影在线看 | 国产精品偷拍 | 97av视频在线| 欧美激情视频一区二区三区 | 国产精品免费观看视频 | 亚洲中文字幕视频一区 | 久草国产在线 | 香蕉av网站 | 不卡一区二区三区四区 | 日韩午夜av | 四虎在线免费视频 | 女人18毛片毛片毛片毛片区二 | av黄页| www.涩涩爱| 色成人免费网站 | 一级性生活毛片 | 欧美在线视频一区二区三区 | 在线a视频 | 日韩爱爱片 | 欧洲美女粗暴牲交免费观看 | 青青草视频观看 | 国产91精品久久久久久久网曝门 | 久草av在线播放 | 成人免费毛片色戒 | av青青| 亚洲永久精品国产 | 大尺度摸揉捏胸床戏视频 | 岛国中文字幕 | 亚洲精品aⅴ中文字幕乱码 亚洲狠狠婷婷综合久久久久图片 | www.99爱| 日本久久久久 | 葵司有码中文字幕二三区 | 国产一区二区网址 | 加勒比色综合 | 免费在线观看的av | 国产精品福利一区二区三区 | 欧美黄色小说 | 国产免费叼嘿网站免费 | 999毛片 | 蜜乳av一区二区 | 精品一区91 | 特黄做受又粗又大又硬老头 | av波多野吉衣 | 国产亚洲精品电影 | 日本精品99 | 日韩精品免费一区二区在线观看 | 久久伊人一区二区 | 五月天激情视频在线观看 | 日韩一级黄 | 一本一道久久a久久精品蜜桃 | 一个色综合网站 | 91视频综合网 | 黑人巨大精品欧美 | 久久久久国产精品午夜一区 | 成人短视频在线播放 | 亚洲av少妇一区二区在线观看 | 国产一区二区在线视频 | 国产伦精品一区二区三区妓女下载 | 亚洲视频一区 | 五月婷婷激情网 | 手机看片国产 | 色哟哟精品一区 | 黄色福利| 日本欧美www| 一本大道久久a久久精二百 中文字幕日韩精品亚洲一区小树林 | 草草在线免费视频 | 欧美性猛交ⅹ乱大交3 | 好了av在线 | 国产精品白嫩极品美女视频 | 天堂av在线免费观看 | 国产三级小视频 | 在线观看三区 | 成年男女免费视频 | 人妻av中文系列 | 色肉色伦交av色肉色伦 | 国模在线 | 青青草在线播放 | 黄色一级免费观看 | 久久综合国产精品 | 欧美成人精品一区二区三区在线观看 | 国产精品久久久精品 | 538任你躁在线精品免费 | 欧美成人免费观看视频 | 成年人在线观看视频免费 | 日本一区不卡在线观看 | 亚洲美女福利视频 | av观看国产 | 天天色天天插 | 999视频在线播放 | 中文字幕超清在线观看 | 久久久免费电影 | 91av高清| 亚洲天堂成人在线 | 欧美熟妇毛茸茸 | 日本va在线| 春色av | 亚洲精品高清视频 | 日韩综合在线 | 96av在线视频 | 中文字幕免费看 | 韩国三级hd中文字幕叫床浴室 | 男女日批在线观看 | 人人草人人插 | 亚洲一卡二卡三卡四卡 | 国产一区二区三区在线播放无 | 黄大色黄大片女爽一次 | 久草网在线视频 | 少妇被躁爽到高潮 | 国产无遮挡裸体免费视频 | 91日韩中文字幕 | 久久在线免费视频 | 亚洲av最新在线网址 | 亚洲乱码中文字幕 | 东京热无码av一区二区 | 国产黄色一区二区 | www.黄色免费 | 91高清视频在线观看 | 亚洲第一页色 | 秋霞精品| 美女视频久久久 | 亚洲黄色免费观看 | 丁香激情综合 | 国产一区二区三区乱码 | 草草影院第一页 | 中文字幕在线播放第一页 | 91大神精品在线| 肉大捧一进一出免费视频 | 亚洲伦理在线视频 | 99自拍视频| 久久久综合网 | 巨乳女教师的诱惑 | 日本中文在线观看 | 免费观看av毛片 | 老熟女毛茸茸 | 男女裸体无遮挡做爰 | 亚洲自拍偷拍一区 | 成人免费在线播放视频 | 黑人玩弄人妻一区二 | 秋霞免费av | 秋霞综合网 | 波多野结衣不卡视频 | 三级三级久久三级久久 | 日韩国产在线一区 | 亚洲一区有码 | 欧美大片aaa | 亚洲欧美日韩精品久久亚洲区 | 色妺妺视频网 | 一区二区三区 日韩 | 国产午夜无码视频在线观看 | 天堂网a| 日韩三级免费看 | 国产三区精品 | 日韩在线精品强乱中文字幕 | 亚洲免费视频一区 | 国产欧美又粗又猛又爽 | 国产日韩欧美中文字幕 | 久久受 | 亚洲综合成人网 | 法国空姐 在线 | 久久久久亚洲av片无码下载蜜桃 | 两个人做羞羞的视频 | 国产xxxx做受视频 | 欧美性受视频 | 亚洲综合免费视频 | 欧美日韩生活片 | 91麻豆精品国产 | 性生活av | 777精品伊人久久久久大香线蕉 | 色av网 | 男生操女生逼逼 | 久久久国| 免费国产一区 | 极品少妇xxxx| 久久久成人av | 中文字幕在线视频一区二区 | 国产视频在线观看一区二区 | 日韩成人免费视频 | 日韩在观看线 | 中文字幕av专区 | 欧美又粗又长又爽做受 | 久久sese| 黄色的网站免费看 | 国产精品一品二区三区的使用体验 | 午夜精品久久久久久久99老熟妇 | 国产综合久久 | 亚洲AV无码阿娇国产精品 | 精品一区二区久久久久蜜桃 | 日本亚洲色图 | 成人午夜激情视频 | 日本国产一级片 | 欧美在线一级片 | 精品一区二区在线播放 | 天天爽天天爽天天爽 | 中文字幕 国产 | 国产亚洲电影 | 国产日日夜夜 | 日本不卡影院 | 在线观看国产麻豆 | 一级黄色美女视频 | 久操视频网 | 夜夜嗨网站 | 三年大片在线观看 | 日韩一区二区三区三四区视频在线观看 | 中文字幕在线一区二区三区 | 天天噜夜夜噜 | 在线免费视频一区 | 成人免费一区二区三区在线观看 | 亚洲熟妇丰满大屁股熟妇 | www九色| 日韩精品在线一区 | 国产日产欧美一区二区三区 | 中文字幕一区视频 | 国产午夜精品一区二区三区视频 | 国产精品jizz | 美女视频在线免费观看 | 一区二区三区四区av | 夜夜骑天天干 | 亚洲成人动漫在线观看 | 天堂av√ | 30一40一50女人毛片 | 665566综合网| 性欧美长视频 | av网站国产| 亚洲人体一区 | 欧美成人怡红院 | 18av在线视频 | 精品欧美乱码久久久久久1区2区 | 欧美成人不卡视频 | 四房婷婷 | 成人av电影免费观看 | 久久久免费高清视频 | 日韩精品av一区二区三区 | 18禁网站免费无遮挡无码中文 | 成年人91视频| 欧美老肥妇做爰bbww | 欧美高清hd19 | 红桃视频91| 中文字幕电影一区二区 | 国产成人短视频 | 日韩黄色av网站 | 精品一区二区三区毛片 | 成人精品二区 | a级黄色网址 | 欧美日韩一级在线观看 | 欧美视频在线观看一区二区三区 | 中文在线观看免费 | 丝袜视频在线 | 96av视频| 黑人精品无码一区二区三区 | 亚洲一区二区三区四区av | 黄色片在线观看视频 | 日本大胆欧美人术艺术 | 亚洲中文字幕在线观看 | 在线看片一区二区 | 一区二区精品在线 | 亚洲永久免费视频 | 爱啪啪av | 亚洲免费片 | 欧美一级在线观看 | 五月婷婷六月合 | 国产精品videossex国产高清 | 免费观看av | 欧美一级淫片aaaaaa | 麻豆入口 | 青青草97| 97人人在线视频 | 高潮在线视频 | 亚洲第一男人天堂 | 午夜寻花| 日本一区二区三区在线播放 | 一级淫片免费 | 少妇太紧太爽又黄又硬又爽 | 久久综合久久久 | 夜夜操夜夜 | 香蕉视频911| 丁香六月欧美 | 麻豆视频传媒入口 | 久久久国产亚洲 | 中文字幕成人av | 久久久久久久美女 | 日韩欧美偷拍 | 中文字幕欧美人妻精品一区蜜臀 | 李宗瑞91在线正在播放 | 好吊色欧美一区二区三区视频 | 福利资源导航 | 国产毛片毛片毛片 | av手机免费看 | 女人18毛片水真多 | 婷婷伊人网 | 久操网在线 | 亚洲免费一区 | 不卡的一区二区 | 国产女人在线视频 | 亚欧美在线观看 | 国产精品国语自产拍在线观看 | 国产精品无码久久久久高潮 | 免费a网址 | 一区二区在线免费视频 | 国产人妻大战黑人20p | 青青草草视频 | 麻豆视频一区二区 | 日韩一区二区精品 | 中文字幕免费高清在线观看 | 又黄又色的网站 | 麻豆av一区二区三区在线观看 | 天堂一二三区 | 69色|