一種基于fpga的電力諧波分析儀的制作方法

            文檔序號:6016115閱讀:280來源:國知局
            專利名稱:一種基于fpga的電力諧波分析儀的制作方法
            技術(shù)領(lǐng)域
            本發(fā)明屬于電力電子技術(shù)領(lǐng)域,特別涉及一種基于FPGA的電力諧波分析儀。
            背景技術(shù)
            電力系統(tǒng)的諧波分析在電能質(zhì)量分析,電力系統(tǒng)控制中具有重要意義,諧波的快速準(zhǔn)確的檢測是諧波抑制和無功補償?shù)幕A(chǔ)。目前,國內(nèi)外對于電力系統(tǒng)諧波檢測的方法很多,國際國內(nèi)針對電力系統(tǒng)諧波檢測和諧波分析的方法主要有以下幾種(1)模擬濾波器法早期的諧波檢測方法都是基于頻域理論,即采用模擬濾波原理。模擬濾波器有兩種,一是通過濾波器濾除基波電流分量,得到諧波電流分量;二是帶通濾波器得出基波分量,再與被檢測電流相減后得到諧波電流分量,其原理和電路結(jié)構(gòu)簡單,造價低,能濾除一些固有頻率的諧波。但這種檢測方法有其自身的缺點誤差大,實時性差,電網(wǎng)頻率變化時尤其明顯,對電路元件參數(shù)十分敏感,參數(shù)變化時檢測效果明顯變差,難以獲得理想的幅頻和相頻特性。(2)基于瞬時無功功率理論的諧波檢測法根據(jù)該理論,可以得到瞬時有功功率ρ和瞬時無功功率q,ρ和q中都含有直流分量和交流分量。由此可得被檢測電流的基波分量,將基波分量與總電流相減即得相應(yīng)的諧波電流。因為該方法忽略了零序分量,且對于不對稱系統(tǒng),瞬時無功的平均分量不等于三相的平均無功。所以該方法只適用于三相電壓正弦對稱情況下的三相電路諧波和基波無功電流的檢測。(3)基于神經(jīng)網(wǎng)絡(luò)的諧波檢測法隨著電網(wǎng)諧波檢測要求的速度、精度和實時性不斷提高,諧波檢測理論的不斷深入,人工神經(jīng)網(wǎng)絡(luò)等新型的方法正在被運用到諧波檢測中。人工神經(jīng)網(wǎng)絡(luò)包括自適應(yīng)線性人工神經(jīng)網(wǎng)絡(luò),基于多層前饋的神經(jīng)網(wǎng)絡(luò)和基于多層反饋的神經(jīng)網(wǎng)絡(luò)。用人工神經(jīng)網(wǎng)絡(luò)實現(xiàn)諧波與無功電流檢測不僅對周期性變化的電流具有很好的跟蹤性能而且對各種非周期變化的電流也能進(jìn)行快速跟蹤,對高頻隨機干擾有良好的識別能力。但是神經(jīng)網(wǎng)絡(luò)工程實際還有很多問題,例如,沒有規(guī)范的人工神經(jīng)網(wǎng)絡(luò)構(gòu)造方法,需要大量的訓(xùn)練樣本,如何確定需要的樣本數(shù)沒有規(guī)范方法,人工神經(jīng)網(wǎng)絡(luò)的精度對樣本有很大的依賴性等。(4)基于小波分析的諧波檢測法小波分析是一種新興起的信號處理方法,是信號時域分析的重要工具,它克服了傅立葉分析在頻域完全局部化而在時域完全無局部化的缺點,在頻域和時域同時具有良好的局部性,對信號的奇異點(如突變點、不連續(xù)點)也非常敏感,因而能算出某一特定時間的頻率分布,并將各種不同頻率組成的頻譜信號分解成不同頻率的信號塊,尤其適合非平穩(wěn)突變信號的分析與處理。通過對含有諧波的電流信號進(jìn)行正交小波分解,得到電流信號在各個頻段的分解結(jié)果,可較準(zhǔn)確地求出基波電流,進(jìn)而求得諧波,從而實時地跟蹤電力系統(tǒng)諧波的變化。雖然小波變換得到結(jié)果的解析度是好的,但處理時間和內(nèi)存使用的要求很
            尚ο(5)傅里葉變換法傅里葉變換法是在電力系統(tǒng)諧波分析中應(yīng)用非常廣泛的一種方法。傅里葉變換方法具有完備、正交等許多優(yōu)點,是一種經(jīng)典的信號分析方法。傅里葉變換的實質(zhì)是將分析信號分解為多個頻率的正弦分量之和的形式。在實際分析時,需要用到計算機進(jìn)行諧波分析。這就需要信號在時域和頻域上是離散的,并且都是有限長。此時的傅里葉變換,稱之為離散傅里葉變換(DFT)。為了提高運算速度,一般在進(jìn)行諧波分析時主要是應(yīng)用快速傅里葉變換算法(FFT)。實際分析的過程中,由于柵欄效應(yīng)和泄漏現(xiàn)象的影響,F(xiàn)FT會產(chǎn)生很大的誤差。這主要是因為將信號截斷為有限長以及非同步采樣。針對這兩個問題,可以采取不同的措施來減小分析誤差。通過加窗函數(shù)的方法來緩解泄漏現(xiàn)象,但是要選擇合適的窗函數(shù)才會達(dá)到目的。通過插值算法來應(yīng)對柵欄效應(yīng),這樣可以分析出被漏掉的頻率成份。加窗插值FFT中的窗函數(shù)已經(jīng)發(fā)展到矩形窗(泄漏現(xiàn)象比較明顯),Hamming窗、Harming窗、 Blackman窗等。文獻(xiàn)[25]中,為了限制泄漏的影響,經(jīng)典的Harming窗取代了矩形窗。雖然窗函數(shù)在一定程度上能夠減少泄漏的影響,并且DFT的準(zhǔn)確性已經(jīng)有所改善,但是由于寬大的主瓣,窗函數(shù)也會模糊DFT的結(jié)果。然后運用插值算法,有效地較少柵欄效應(yīng)。本發(fā)明依據(jù)傅里葉變換法來實現(xiàn)諧波的檢測。目前,在國內(nèi)外的研究中,主要通過在三種硬件平臺上來實現(xiàn)諧波的檢測(1)以虛擬儀器為核心這種方式是以通用的PC機為主要的開發(fā)平臺,配合虛擬儀器軟件,數(shù)據(jù)采集卡以及傳感器電路實現(xiàn)。例如,美國國家儀器公司的虛擬儀器軟件LabVIEW具有強大的模塊節(jié)點和計算能力,網(wǎng)絡(luò)功能。廈門大學(xué)的陳栩等人提出了一種采用基于LabVIEW的電能質(zhì)量監(jiān)控系統(tǒng)的設(shè)計方法,它使用FFT算法作為諧波檢測算法。這種實現(xiàn)方法靈活,擺脫了硬件上的束縛和干擾,充分利用了虛擬軟件的功能,提高了測量的精度和系統(tǒng)開發(fā)的效率。但是這種諧波分析儀不能隨身攜帶,并且需要一臺專用的PC機開發(fā)平臺,實現(xiàn)成本較高。(2)以通用處理器為核心這種方式采用通用的處理器(ARM,DSP等)作為控制核心和算法實現(xiàn)單元,并輔以外圍的AD轉(zhuǎn)換電路,信號調(diào)理電路等外圍電路?,F(xiàn)在有很多的電能監(jiān)控設(shè)備或者諧波分析儀都采用具有很強控制功能的DSP或者DSP和ARM雙CPU作為核心。這種方式既能很好滿足檢測算法中復(fù)雜運算的需求,也能夠滿足系統(tǒng)控制的需要,并且具有很好的網(wǎng)絡(luò)功能,能夠和本地的應(yīng)用端網(wǎng)絡(luò)連在一起。上海交通大學(xué)的王麗穎等人使用DSP和ARM雙CPU結(jié)構(gòu)實現(xiàn)了電能質(zhì)量監(jiān)控系統(tǒng)。采用雙CPU的工作模式,這種方式雖然可以很好的實現(xiàn)各種功能,DSP和ARM都存在開發(fā)周期長,系統(tǒng)可擴展性不高的缺點,并且這種雙CPU結(jié)構(gòu)電路的硬件電路復(fù)雜,價格昂貴,設(shè)備開發(fā)和調(diào)試十分困難,后期的升級和維護也十分不易。(3)以 FPGA 為核心這種方式采用FPGA來實現(xiàn)諧波檢測算法,并且加上外圍采集電路等實現(xiàn)。武漢理工大學(xué)的潘磊等人研究了一種基于FPGA的電力諧波檢測系統(tǒng)。使用FPGA實現(xiàn)了 FFT處理單元。FFT處理單元采用了 FPGA內(nèi)嵌的硬件乘法器和加法器做并行運算,實現(xiàn)了硬件的并行處理。但這種檢測系統(tǒng)由于采用硬件實現(xiàn)邏輯,系統(tǒng)可重構(gòu)性較弱,耦合性較高,對于需要添加新的功能模塊或?qū)δ骋荒K進(jìn)行修改時,系統(tǒng)的修改較大,同時需要對系統(tǒng)進(jìn)行重建;另外,該系統(tǒng)的可讀性較差,各模塊間的數(shù)據(jù)流向以及數(shù)據(jù)通訊方式很難直觀的被程序開發(fā)人員讀取,為后續(xù)的維護造成了很大的困難。山東大學(xué)的梁海霞等人研究了一種基于FPGA的諧波檢測裝置,提出利用SOPC技術(shù)解決方案,在FPGA中內(nèi)嵌Mos II軟核作為控制器,并實現(xiàn)了硬件的FFT算法。但該裝置沒有具體的系統(tǒng)實現(xiàn),只是進(jìn)行了功能的仿真,并且該裝置的FFT的處理精度為1 位, 由于采用同步采樣,在實際系統(tǒng)設(shè)計實現(xiàn)中,容易造成采樣不準(zhǔn)確的情況。

            發(fā)明內(nèi)容
            本發(fā)明針對現(xiàn)有技術(shù)存在的問題,提供一種基于SOPC技術(shù),應(yīng)用FPGA平臺的電力諧波分析儀,采用電流電壓互感器作為前端傳感器,使用高速的16位AD芯片設(shè)計采樣電路采集電流電壓信號,運用Verilog HDL硬件描述語言設(shè)計了采樣電路AD7606的驅(qū)動程序, 將采集到的電壓電流信號進(jìn)行緩存,然后進(jìn)行前期的諧波算法預(yù)處理。結(jié)合SOPC Builder, 在FPGA上嵌入Nios II軟核作為系統(tǒng)控制單元,設(shè)計了一個顯示IP核,將這個IP核掛載在NiosII軟核處理器上,控制檢測數(shù)據(jù)在液晶屏上的顯示輸出。完成系統(tǒng)構(gòu)建,并編寫C 程序進(jìn)行控制算法和顯示控制,最后在友晶公司的DE2-70開發(fā)板上進(jìn)行驗證。本發(fā)明的一種基于FPGA的電力諧波分析儀,包括有硬件部分和軟件部分。本發(fā)明的一種基于FPGA的電力諧波分析儀的硬件部分包括有傳感器電路和AD采樣電路。傳感器電路主要負(fù)責(zé)將電力系統(tǒng)中電流電壓信號轉(zhuǎn)換為-5V +5V之間的電壓信號,并完成信號的調(diào)理。AD采樣電路主要負(fù)責(zé)將模擬的諧波信號轉(zhuǎn)換為16位的數(shù)字信號, 提供給FPGA做進(jìn)一步的處理。本發(fā)明的一種基于FPGA的電力諧波分析儀,主要包括以下幾個模塊數(shù)據(jù)采集模塊、數(shù)據(jù)緩存模塊、數(shù)據(jù)類型轉(zhuǎn)換模塊、加窗模塊、FFT運算模塊、運算結(jié)果存儲和中斷產(chǎn)生模塊、Nios II軟核處理器模塊、LTM顯示模塊。本發(fā)明的數(shù)據(jù)采集模塊主要負(fù)責(zé)驅(qū)動AD7606進(jìn)行8路信號,1024點的同步采樣, 然后將數(shù)據(jù)讀取到FPGA片內(nèi),提供給下一級的數(shù)據(jù)緩存模塊。模塊通過對輸入時鐘計數(shù), 控制AD7606進(jìn)行采樣速率為10. 24K的等間隔采樣。本發(fā)明的數(shù)據(jù)緩存模塊主要完成采樣數(shù)據(jù)的存儲,控制,以及輸出。從ADC數(shù)據(jù)采樣模塊輸出的采樣數(shù)據(jù)是按照1到8路數(shù)據(jù)一次輸出,系統(tǒng)需要將這8路數(shù)據(jù)分別存儲在片內(nèi)RAM中,當(dāng)有輸出觸發(fā)信號時,按設(shè)定順序依次輸出1路采樣數(shù)據(jù)。本發(fā)明的數(shù)據(jù)類型轉(zhuǎn)化模塊主要負(fù)責(zé)將數(shù)據(jù)類型為有符號整型的采樣數(shù)據(jù)轉(zhuǎn)換為浮點類型數(shù)據(jù),然后將轉(zhuǎn)換為的浮點數(shù)據(jù)提供給加窗模塊。本發(fā)明的加窗模塊主要負(fù)責(zé)將給每一路的采樣數(shù)據(jù)加上IOM寬度的Blackman 窗。當(dāng)有效信號到達(dá)時,該模塊首先產(chǎn)生Blackman窗數(shù)據(jù),然后將Blackman窗數(shù)據(jù)與采樣數(shù)據(jù)相乘,并將計算結(jié)果提供給下一級的FFT運算模塊。
            本發(fā)明的FFT運算模塊主要實現(xiàn)10M點的浮點FFT運算。該模塊主要負(fù)責(zé)將加 窗后的數(shù)據(jù)進(jìn)行FFT運算,它決定了本發(fā)明整個系統(tǒng)的實時性。本發(fā)明的運算結(jié)果存儲和中斷產(chǎn)生模塊主要負(fù)責(zé)緩存FFT運算結(jié)果,并通知Nios II軟核處理器系統(tǒng)讀取運算結(jié)果。這個模塊產(chǎn)生緩存的控制信號,數(shù)據(jù)存儲地址,以及時序 控制等,存儲完成后產(chǎn)生1個時鐘的高電平,作為中斷信號,與Mos II軟核處理器系統(tǒng)交互。本發(fā)明的Nios II軟核處理器模塊主要完成總體的核心控制,各個模塊的協(xié)調(diào)作 用,負(fù)責(zé)讀取FFT運算結(jié)果,然后在軟件中對采樣數(shù)據(jù)進(jìn)行插值運算,計算出諧波的相關(guān)參 量,調(diào)用LTM顯示IP核,將諧波參量顯示在LTM液晶屏上。當(dāng)收到FFT運算模塊完成計算, 并完成運算結(jié)果的緩存后,產(chǎn)生一個中斷信號,通知Nios II軟核處理器系統(tǒng)讀取運算結(jié) 果。當(dāng)完成1路諧波計算后,該模塊發(fā)出命令,讀取下一路數(shù)據(jù)。直到完成所有電壓電流采 樣數(shù)據(jù)后,再發(fā)出命令開始下一輪的采樣。針對現(xiàn)有技術(shù)存在的問題,本發(fā)明的一種基于FPGA的電力諧波分析儀,前端的諧 波采集和處理使用硬件實現(xiàn),充分發(fā)揮了硬件加速的優(yōu)勢;利用FPGA中的Mos II軟核處 理器對處理后的數(shù)據(jù)做進(jìn)一步運算和控制,該過程在FPGA片內(nèi)總線上完成,解決了數(shù)據(jù)傳 輸?shù)钠款i,充分發(fā)揮了硬件設(shè)計的高速性和Nios II軟核處理器控制的靈活性。在lOOMHz 時鐘下,前端的處理模塊可以在77US內(nèi)完成一次10M點浮點FFT運算,Nios II軟核處理 器能在78ms完成1路諧波信號的插值和諧波參數(shù)計算,在實際應(yīng)用中,系統(tǒng)能夠滿足實時 檢測的需求。本系統(tǒng)的AD采用兩片AD模塊級聯(lián)的方式,具有16位采集精度,具有較高的 精度。對系統(tǒng)進(jìn)行的相關(guān)分析和測試表明,系統(tǒng)運行穩(wěn)定,檢測精度較高,處理速度快。


            圖1本發(fā)明的系統(tǒng)總體構(gòu)成示意圖;圖2本發(fā)明的系統(tǒng)總體流程圖;圖3本發(fā)明實施例使用的Nios II軟核處理器模塊結(jié)構(gòu)框圖;圖4本發(fā)明實施例使用的電壓采集電路圖;圖5本發(fā)明實施例使用的電流采集電路圖;圖6本發(fā)明實施例使用的單片AD7606-4電路原理圖;圖7本發(fā)明實施例使用的REF5025的電路原理圖;圖8本發(fā)明實施例使用的數(shù)據(jù)采集模塊圖;圖9本發(fā)明實施例使用的數(shù)據(jù)采集模塊設(shè)計流程圖;圖10是本發(fā)明實施例使用的AD采樣驅(qū)動單元狀態(tài)機流程圖;圖11本發(fā)明實施例使用的數(shù)據(jù)采集模塊的時序仿真圖;圖12本發(fā)明實施例使用的數(shù)據(jù)緩存模塊圖;圖13本發(fā)明實施例的寫數(shù)據(jù)時序仿真圖;圖14本發(fā)明實施例使用的數(shù)據(jù)類型轉(zhuǎn)換模塊圖;圖15本發(fā)明實施例使用的數(shù)據(jù)類型轉(zhuǎn)換模塊時序仿真圖;圖16本發(fā)明實施例使用的加窗模塊結(jié)構(gòu)圖;圖17本發(fā)明實施例使用的控制單元狀態(tài)機流程圖18本發(fā)明實施例使用的復(fù)數(shù)邏輯單元ALUl結(jié)構(gòu)圖;圖19本發(fā)明實施例使用的復(fù)數(shù)邏輯單元ALU2結(jié)構(gòu)圖;圖20本發(fā)明實施例使用的運算結(jié)果存儲模塊圖;圖21本發(fā)明實施例使用的運算結(jié)果存儲模塊的時序圖;圖22本發(fā)明實施例使用的SOPC Builder的配置界面圖;圖23本發(fā)明實施例的總體控制單元的流程圖;圖M對本發(fā)明實施例FFT運算單元進(jìn)行測試的單頻正弦信號頻譜圖;圖25對本發(fā)明實施例FFT運算單元進(jìn)行測試的單頻正弦信號仿真誤差圖;圖沈?qū)Ρ景l(fā)明實施例FFT運算單元進(jìn)行測試的雙頻正弦信號頻譜圖;圖27對本發(fā)明實施例FFT運算單元進(jìn)行測試的雙頻正弦信號電路仿真誤差圖;圖觀對本發(fā)明實施例FFT運算單元進(jìn)行測試的三角信號頻譜圖;圖四對本發(fā)明實施例FFT運算單元進(jìn)行測試的三角信號電路仿真誤差圖;圖30對本發(fā)明實施例系統(tǒng)實時行進(jìn)行測試的一路諧波計算所需時間圖。
            具體實施例方式下面結(jié)合附圖對本發(fā)明內(nèi)容作進(jìn)一步描述。本發(fā)明實施例是在Altera公司的FPGA開發(fā)板DE2-70上實現(xiàn)的。開發(fā)板預(yù)留了兩個擴展接口。在本系統(tǒng)中,這兩個擴展接口分別用于顯示模塊LTM和AD7606的采集電路連接。DE2-70開發(fā)板上的FPGA芯片是Altera公司Cyclone II系列EP2C70F896C6N。本系統(tǒng)總體構(gòu)成示意圖如圖1所示。系統(tǒng)的總體流程圖如圖2所示。本系統(tǒng)中的Mos II軟核處理器模塊結(jié)構(gòu)框圖如圖3所示。本發(fā)明的基于FPGA的電力諧波分析儀中的傳感器電路,主要包括電壓和電流信號采集部分。在電壓信號采集部分中,電壓的輸入范圍為士 IOOOVdf 士 1000V范圍內(nèi)的交流電壓轉(zhuǎn)換為士5V范圍內(nèi)的交流電壓信號。在電流信號采集部分,電流的輸入范圍為士20A,將士20A范圍內(nèi)的交流電流轉(zhuǎn)換為士5V范圍內(nèi)的交流電壓信號。本發(fā)明采用電壓互感器PT,電流互感器CT作為電壓采集部分和電流采集部分的傳感器。本發(fā)明選用“兵字”高精度的電壓互感器和電流互感器TV1013和TA1420。這兩種傳感器的精度都一樣,工作頻率范圍20Hz 20kHz,使用運算放大器的情況下,互感器信號輸出非線性度不大于0. 1%,相位誤差小于5’,滿足系統(tǒng)的精度和測量范圍的要求。(1)電壓信號采集電路電壓信號采集電路的輸入為士 1000V范圍內(nèi)的交流電壓,具體的設(shè)計電路如圖4 所示。電壓互感器TV1013的電流型電壓互感器,最大輸入電壓為1000V,互感器線圈匝數(shù)比為1 1,額定電流2mA。運算放大器采用了 TI的0PA2227該放大器具有高精度、低噪聲、 低功耗、低偏置的優(yōu)點,能夠滿足要求。由電壓互感器的原理可以知道,副邊的電壓輸出為U。ut = UinZR1^R2,因為電壓互感器的額定電流為2·,!^的阻值設(shè)定為500K,I 2的阻值設(shè)定為2. 5K,能夠滿足采樣電路輸入士 1000V交流電時,輸出士5V的交流電壓信號。電壓信號采集電路的輸出將送到AD7606 的模擬輸入端。(2)電流信號采集電路
            電流信號采集電路的輸入為士20A范圍內(nèi)的交流電流,具體的電路如圖5所示。電流互感器互感器TA1420-04線圈匝數(shù)比為4000 1,原邊額定輸入電流20A,次級的額定輸出電流為5mA。信號調(diào)理電路中,放大器依然采用0PA2227。通過電流互感器的參數(shù)可以知道,副邊電壓輸出為U。ut = Iin/4000*R,電流互感器的額定輸出電流為5mA,R的阻值設(shè)定為1K,能夠滿足采樣電路輸入士20A時,輸出士5V 的交流電壓信號。同樣地,將電流信號采集電路的輸出送到AD7606的模擬輸入端。本發(fā)明的AD采樣芯片選用了 ADI公司的AD7606-4,具體的電路原理如圖6所示。 AD7606-4具有16位數(shù)據(jù)精度,4路同步采樣輸入,同時所有通道均能以高達(dá)200kSPS的吞吐速率采樣,能夠滿足系統(tǒng)對精度和采樣速度的要求。AD7606-4片內(nèi)資源非常豐富,可以最大限度簡化采樣系統(tǒng)外圍電路的設(shè)計。它采用5V單電源供電,可以處理士 IOV或士5V范圍內(nèi)的雙極性輸入信號輸入,片內(nèi)箝位保護電路可以耐受最高達(dá)士 16. 5V的電壓,模擬輸入阻抗均為1ΜΩ。同時,片內(nèi)還集成了二階抗混疊濾波器,抗混疊濾波器的3dB截止頻率為22KHz。因此,在AD采樣電路之前省略抗混疊濾波器的設(shè)計。因為AD7606只有4路采樣通道,為了滿足3相電采樣的需求,系統(tǒng)使用兩片 AD7606-4并聯(lián),使系統(tǒng)可以實現(xiàn)多達(dá)8路信號的同步采樣。在本發(fā)明中,選用士5V輸入范圍,片外參考電壓方式。與FPGA的通信時,使用并行通信方式。發(fā)明中使用片外的基準(zhǔn)電壓REF5025作為AD7606的參考電壓,REF5025具有低噪聲、極低漂移、高精度的優(yōu)點,能夠使系統(tǒng)精度更加穩(wěn)定。圖7是REF5025的電路原理圖,輸入5V電壓,輸出2. 5V的標(biāo)準(zhǔn)參考電壓。 實施例的邏輯電路模塊數(shù)據(jù)采集模塊主要完成對芯片AD7606的驅(qū)動,并將采樣數(shù)據(jù)讀取到FPGA片內(nèi),提供給下一級的數(shù)據(jù)緩存模塊。當(dāng)數(shù)據(jù)采集模塊收到Mos II軟核處理器系統(tǒng)的啟動采樣信號時,數(shù)據(jù)采集模塊開始啟動新的一輪采樣。在一輪采樣中,數(shù)據(jù)采樣模塊將等間隔采樣 IOM次,采樣速率為10. 24KHz0模塊的時鐘端iCLK輸入為50MHz,計數(shù)器對時鐘輸入計數(shù), 從0計到4883時,即完成一次采樣,并啟動下一次采樣。通過硬件自身計數(shù)實現(xiàn)的10. 24KHz 等間隔采樣要比在處理器系統(tǒng)中使用計數(shù)器啟動每次采用精確。數(shù)據(jù)采集模塊的頂層設(shè)計文件如圖8所示。數(shù)據(jù)采集模塊需要對兩片AD7606-4進(jìn)行驅(qū)動。AC0NVSTA和AC0NVSTB用于啟動第一片AD7606采樣,用于控制電壓采樣的4個通道;BC0NVSTA和BC0NVSTB用于啟動第二片 AD7606采樣,用于控制電流采樣的4個通道。在本系統(tǒng)中,iSwitch用于控制電流電壓的相位偏移,也就是控制AC0NVST和BC0NVST之間的時間差。由于傳感器精度的原因,經(jīng)過傳感器之后的電流電壓信號存在系統(tǒng)性相位誤差,在本模塊中,可以精確控制同一次采樣中電流電壓通道的時間差,校正相位差。當(dāng)iSwitch的輸入為邏輯'0'時,電流電壓通道同時采樣,當(dāng)iSwitch的輸入為邏輯'1',可以設(shè)計電流電壓的采樣通道采樣的時間差,這個參數(shù)由傳感器精度決定。每一輪采樣的流程如圖9所示。在硬件電路的設(shè)計中已經(jīng)知道AD7606的數(shù)據(jù)輸出是三態(tài)的,兩片AD7606公用一條數(shù)據(jù)總線。當(dāng)AD7606完成數(shù)據(jù)轉(zhuǎn)化后,CS_N和RD_N信號配合使用就可以將轉(zhuǎn)換后的數(shù)據(jù)讀取。當(dāng)ACS_N為低電平,RD_N下降沿時,第一片AD7606返回一路轉(zhuǎn)換數(shù)據(jù),4個RD_N下降沿將返回第一片AD7606全部4個通道的采樣數(shù)據(jù);同樣的,當(dāng)BCS_N為低,RD_N下降沿時,第二片AD7606返回轉(zhuǎn)換一路數(shù)據(jù),4個RD_N下降沿將返回第二片AD7606全部4個通道的采樣數(shù)據(jù)。oCLK,ofeite,oDATA是模塊向數(shù)據(jù)緩存模塊輸出的信號和數(shù)據(jù)。在一次數(shù)據(jù)讀取中,在開始輸出數(shù)據(jù)時。ofeite變?yōu)楦唠娖剑?dāng)輸出完8路數(shù)據(jù)后,oWrite變?yōu)榈碗娖?。在oCLK的上升沿,oDATA開始輸出鎖存數(shù)據(jù)。圖10是AD采樣驅(qū)動單元狀態(tài)機流程圖。圖11是數(shù)據(jù)采樣模塊的時序圖。數(shù)據(jù)的緩存模塊主要完成采樣數(shù)據(jù)的存儲,控制,以及輸出,所以該模塊由主要分為3個部分,分別是數(shù)據(jù)的存儲部分、數(shù)據(jù)輸出控制部分、數(shù)據(jù)輸出讀取部分。在數(shù)據(jù)緩存模塊中,有8個長度為IOM的雙口 RAM,用于存儲采樣電路8個通道中的數(shù)據(jù)。數(shù)據(jù)緩存模塊頂層設(shè)計文件如圖12所示。在系統(tǒng)中,數(shù)據(jù)采集模塊的oWrite,oCLK, oDATA,分別與 write_en,data—clk,data—in 相連。數(shù)據(jù)的存儲部分是完成8路數(shù)據(jù)的存儲。數(shù)據(jù)采樣模塊輸出的采樣數(shù)據(jù)按照1到 8路通道依次輸出,系統(tǒng)需要將這8路數(shù)據(jù)分別存儲在對應(yīng)的片內(nèi)RAM中。系統(tǒng)需要為每一片雙口 RAM產(chǎn)生寫使能信號,寫地址。寫地址是由一個10M進(jìn)制計數(shù)器產(chǎn)生,計數(shù)器在每個written信號的下降沿加1。因為written信號表示一幀信號有效,也就是說數(shù)據(jù)采集模塊在輸出8路采樣數(shù)據(jù)時,written信號為邏輯'1',在未輸出數(shù)據(jù)時為邏輯'0'。 因此可以直接將這個計數(shù)器作為每片RAM的寫地址。寫使能信號由一個狀態(tài)機和一個8進(jìn)制計數(shù)器共同產(chǎn)生。在采樣數(shù)據(jù)有效輸出時,也就是written為邏輯'1'時,8進(jìn)制計數(shù)器在每個data_Clk信號的上升沿加'1'。狀態(tài)機根據(jù)8進(jìn)制計數(shù)器的大小,分別為每一片RAM產(chǎn)生寫使能信號。當(dāng)寫完10M幀數(shù)據(jù)后,模塊產(chǎn)生一個done_valid的有效脈沖,表示完成了一輪數(shù)據(jù)的存儲。圖13是寫數(shù)據(jù)的時序圖。數(shù)據(jù)的讀取部分是按照指定順序讀取采樣數(shù)據(jù)。outlay用于控制數(shù)據(jù)輸出方式,當(dāng)OUtjay*' 00'時,只輸出A相的電壓電流,即用于檢測單相諧波;當(dāng)outlay 為1 01',用于檢測三相三線制電路的諧波,當(dāng)outlay*' 10'時,用于檢測三相四線制電路的諧波。當(dāng)datajutjrigger有脈沖輸入時,模塊將開始輸出1個通道的采樣數(shù)據(jù)。 在對單相電路進(jìn)行諧波檢測時,在第一個脈沖輸入data_0ut_trigger,A相的IOM個電壓采樣數(shù)據(jù)將開始從data_0ut輸出,同時data_out_Valid有效;在第二個脈沖輸入data_ out_trigger, A相的10 個電流采樣數(shù)據(jù)將開始從data_out輸出,同時data_out_valid 有效。在對3相3線制電路和3相4線制電路進(jìn)行檢測時,數(shù)據(jù)的輸出方式類推。在數(shù)據(jù)輸出部分主要由兩個狀態(tài)機和一個10M進(jìn)制的計數(shù)器實現(xiàn)。1個狀態(tài)機產(chǎn)生讀使能信號, 1個狀態(tài)機讀取RAM的輸出,并賦值給data_0ut。10 進(jìn)制的計數(shù)器主要用于產(chǎn)生RAM的讀取地址。數(shù)據(jù)類型轉(zhuǎn)化模塊主要將有符號的16位整數(shù)轉(zhuǎn)換成滿足IEEE-7M標(biāo)準(zhǔn)的浮點數(shù)。為了保證計算精度,系統(tǒng)在計算過程的數(shù)據(jù)類型采用單精度的浮點。所有AD7606轉(zhuǎn)后的數(shù)據(jù)都是16位的有符號整數(shù),在采樣數(shù)據(jù)進(jìn)行加窗之前,需要將采樣數(shù)據(jù)轉(zhuǎn)換成浮點數(shù)據(jù)。圖14是數(shù)據(jù)轉(zhuǎn)換模塊的框圖,該模塊是調(diào)用軟件開發(fā)工具Qimrtus II提供的模塊實現(xiàn),使用流水線結(jié)構(gòu),轉(zhuǎn)換過程有6個時鐘的延時。圖15是數(shù)據(jù)類型轉(zhuǎn)換模塊時序仿真。加窗模塊主要實現(xiàn)在采樣數(shù)據(jù)進(jìn)行FFT運算之前,對采樣數(shù)據(jù)加上Blackman窗。 加窗模塊的結(jié)構(gòu)圖如圖16所示。加窗模塊主要包括Blackman窗數(shù)據(jù)產(chǎn)生單元,流水線乘法單元。Blackman窗數(shù)據(jù)產(chǎn)生單元將產(chǎn)生長度為IOM的Blackman窗數(shù)據(jù)。因為Blackman 窗是對稱的,所以Blackman窗的數(shù)據(jù)可以復(fù)用。系統(tǒng)采用復(fù)用方式,只將前513個點的數(shù)據(jù)存儲在片內(nèi)ROM中。片內(nèi)ROM的地址是有一個是由一個計數(shù)器產(chǎn)生。在data_in_valid 有效時,在每個時鐘上升沿,計數(shù)器從0計到512,然后512到1,即完成了運算。流水線乘法單元模塊主要實現(xiàn)對應(yīng)的Blackman窗數(shù)據(jù)與采樣數(shù)據(jù)的乘法運算。 這個單元主要調(diào)用了 Quartus II提供的浮點數(shù)乘法模塊實現(xiàn),具有6級流水結(jié)構(gòu),運算結(jié)果對輸入數(shù)據(jù)有6個時鐘的延時。在本模塊中,確保采樣數(shù)據(jù)和Blackman窗數(shù)據(jù)的同步是關(guān)鍵。FFT運算模塊采用基-4方式,本發(fā)明的FFT模塊采用浮點運算,需要占用較多的邏輯資源。由于FPGA片內(nèi)資源的限制,F(xiàn)FT運算模塊不能采用流水結(jié)構(gòu),蝶形運算單元、復(fù)數(shù)乘法器,以及旋轉(zhuǎn)因子都采用復(fù)用形式。FFT運算模塊的設(shè)計主要包括總體控制單元,蝶形運算單元和地址產(chǎn)生單元??傮w控制單元主要完成總體的時序控制。包括對五級蝶形運算單元,地址產(chǎn)生單元,輸出模塊的是時序控制??傮w控制單元包括一個狀態(tài)機,一個7進(jìn)制計數(shù)器。狀態(tài)機共有4個狀態(tài),配合7進(jìn)制計數(shù)器,控制FFT運算模塊的進(jìn)程。這4個狀態(tài)分別為ctrljdle, stage, stage_idle, stage6。狀態(tài)機的設(shè)計如圖17所示。ctrl_idle是空閑狀態(tài),stage狀態(tài)控制蝶形運算的進(jìn)行,stage_idle等待蝶形運算的完成,stage6狀態(tài)用于控制最后FFT運算結(jié)果的輸出。為了提高系統(tǒng)速度,電路的設(shè)計都是采用流水結(jié)構(gòu),蝶形運算中的乘法器,加法器都是采用流水線結(jié)構(gòu),計算結(jié)果相對計算輸入有時間延時。在本設(shè)計中,蝶形單元加旋轉(zhuǎn)因子一共延時35個時鐘。本系統(tǒng)中Stage_ idle延時50個時鐘,重現(xiàn)進(jìn)入下一個狀態(tài)。當(dāng)IOM個采樣數(shù)據(jù)在FFT中的存儲后完成之后,控制單元狀態(tài)機從空閑的ctrl_ idle狀態(tài)進(jìn)入到控制狀態(tài)stage,同時7進(jìn)制計數(shù)器變?yōu)?1',等待IOM個時鐘長度后, 狀態(tài)機進(jìn)入stagejdle狀態(tài),等待第一級蝶形運算的完成。在stage狀態(tài)中,同時產(chǎn)生1個 1024脈沖寬度的使能控制信號。這個使能信號配合7進(jìn)制計數(shù)器用于產(chǎn)生讀寫時能信號, 并提供給地址產(chǎn)生模塊用于產(chǎn)生讀寫地址,以及控制蝶形單元。在將跳出stagejdle狀態(tài)時,判斷7進(jìn)制計數(shù)器的值,如果7進(jìn)制計數(shù)器不等于5,則重新進(jìn)入stage狀態(tài),開始進(jìn)入下一級蝶形圖運算,且計數(shù)器的值加1。如果計數(shù)器等于5,說明已經(jīng)全部完成了 5級蝶形運算,狀態(tài)機進(jìn)入stage6狀態(tài),開始輸出FFT運算結(jié)果。狀態(tài)機再次進(jìn)行stage狀態(tài)時,運行過程和第一級類似,不在此贅述。狀態(tài)機進(jìn)入Stage6狀態(tài)時,將產(chǎn)生一個IOM個脈沖寬度doutjalid信號,用于控制地址產(chǎn)生單元產(chǎn)生FFT輸出地址,以及讀取使能信號。蝶形運算單元,為了減少系統(tǒng)硬件資源的開銷,本系統(tǒng)的蝶形單元采用流水結(jié)構(gòu)。 蝶形產(chǎn)生單元主要包括3個狀態(tài)機,4個復(fù)數(shù)邏輯運算單元。通過狀態(tài)機控制蝶形單元中的延時。蝶形運算單元的復(fù)數(shù)邏輯運算部分,蝶形運算可以寫成如下形式X (K) = A+BffP+Cff2P+Dff3PX (K+N/4) = A-jBffP-Cff2P+jDff3PX (K+2N/4) = A-jBffP+Cff2P-jDff3PX (K+3N/4) = A+j Bffp-Cff2p-j Dff3p將計算式分割為如下
            χ' (η) = (χ (η)+χ (n+N/2)) + (χ (η+Ν/4)+χ (n+3N/4))χ' (Κ+Ν/2) = (χ (η)+χ (η+Ν/2)) -(χ (η+Ν/4)+χ (η+3Ν/4))χ' (Κ+Ν/4) = (χ (η) -χ (η+Ν/2)) -j* (χ (η+Ν/4) -χ (η+3Ν/4))χ' (Κ+3Ν/4) = (χ (η) -χ (η+Ν/2)) +j* (χ (η+Ν/4) -χ (η+3Ν/4))令a (0) =χ (η) +χ (η+Ν/2) ,a (2) = χ (η) -χ (η+Ν/2),a (1) = χ (η+Ν/4) + (η+3Ν/4), a (3) = χ(η+Ν/4)-(η+3Ν/4)上面的等式可以變?yōu)棣? (η) = b(0) = a(0)+a(l)χ' (n+N/2) = b (2) = a(0)-a(l)x' (n+N/4) = b(l) = a(2)_ja(3)x' (n+3N/4) = b(3) = a(2)+ja(3)通過分析可以看出,蝶形單元最終可以分解為只有a+b,a-b, a+jb, a-jb四種計算。本文設(shè)計了兩種復(fù)數(shù)邏輯運算單元用于實現(xiàn)上述4種運算。這兩種復(fù)數(shù)邏輯運算單元為ALUl禾口 ALU2, ALUl完成a+b禾口 a-jb運算,ALU2完成a-b禾口 a+jb運算。ALUl的結(jié)構(gòu)設(shè)計如圖18所示。OP用于控制ALUl的運算方式,當(dāng)op = 0時,c = a+b,或當(dāng)op = 1時,c = a_jb。 ALUl中的浮點加法器具有8級流水結(jié)構(gòu),從數(shù)據(jù)輸入到計算結(jié)果的輸出有8個時鐘的延時。ALU2的結(jié)構(gòu)設(shè)計如圖19所示,op用于控制ALUl的運算方式,當(dāng)op = 1時,c = a+jb,或當(dāng)op = 0時c = a-b。和ALUl —樣,ALU2中的浮點加法器也是具有8級流水結(jié)構(gòu),從數(shù)據(jù)輸入到計算結(jié)果的輸出有8個時鐘的延時。在蝶形單元內(nèi)部,采用流水線結(jié)構(gòu)。a(0),a(l),W2),a(3)在第一級流水中實現(xiàn), 使用第1個狀態(tài)機控制,b (0),b (1),b (2),b (3)在第二級流水中實現(xiàn),使用第2個狀態(tài)機控制,第3個狀態(tài)機主要讀取第二級狀態(tài)機的運算結(jié)果。表1是狀態(tài)機1的工作流程。狀態(tài)機中一共有4個工作狀態(tài),每個狀態(tài)都需要讀取輸入蝶形單元的四個數(shù)。在第3個狀態(tài)時,分別使用一個ALUl和ALU2,讀取χ (η)和 χ (η+Ν/2)計算a (0)和a (2)。同樣的,在第四個狀態(tài)中,計算出a (1)和a (3)。這個狀態(tài)機主要控制將蝶形單元的4個輸入數(shù)據(jù)送到ALUl和ALU2輸入端口上。
            權(quán)利要求
            1.一種基于FPGA的電力諧波分析儀,其特征在于包括傳感器電路和AD采樣電路,還包括設(shè)在FPGA中的數(shù)據(jù)采集模塊、數(shù)據(jù)緩存模塊、數(shù)據(jù)類型轉(zhuǎn)換模塊、加窗模塊、FFT運算模塊、運算結(jié)果存儲和中斷產(chǎn)生模塊、Nios II軟核處理器模塊、LTM顯示模塊,傳感器電路將電力系統(tǒng)中電流電壓信號轉(zhuǎn)換為_5V、5V之間的電壓信號,并完成信號的調(diào)理;AD采樣電路將模擬的諧波信號轉(zhuǎn)換為16位的數(shù)字信號,提供給FPGA做進(jìn)一步的處理;數(shù)據(jù)采集模塊驅(qū)動AD采樣電路進(jìn)行采樣,并將數(shù)據(jù)讀取到FPGA內(nèi),提供給數(shù)據(jù)緩存模塊;數(shù)據(jù)緩存模塊完成采樣數(shù)據(jù)的存儲、控制及輸出;數(shù)據(jù)類型轉(zhuǎn)換模塊將數(shù)據(jù)類型為有符號整型的采樣數(shù)據(jù)轉(zhuǎn)換為浮點類型數(shù)據(jù),然后將轉(zhuǎn)換后的浮點類型數(shù)據(jù)提供給加窗模塊;加窗模塊將給每一路的采樣數(shù)據(jù)加上1024寬度的Blackman窗,當(dāng)有效信號到達(dá)時,該模塊首先產(chǎn)生 Blackman窗數(shù)據(jù),然后將Blackman窗數(shù)據(jù)與采樣數(shù)據(jù)相乘,并將計算結(jié)果提供給下一級的 FFT運算模塊;FFT運算模塊將加窗后的數(shù)據(jù)進(jìn)行IOM點的浮點FFT運算;運算結(jié)果存儲和中斷產(chǎn)生模塊緩存FFT運算結(jié)果,并通知Nios II軟核處理器模塊讀取運算結(jié)果,該模塊產(chǎn)生緩存的控制信號、數(shù)據(jù)存儲地址,以及時序控制,完成存儲后產(chǎn)生1個時鐘的高電平,作為中斷信號,與Nios II軟核處理器模塊交互;Nios II軟核處理器模塊實現(xiàn)總體的核心控制及各個模塊的協(xié)調(diào),讀取FFT運算結(jié)果,對采樣數(shù)據(jù)進(jìn)行插值運算,計算出諧波的相關(guān)參量,調(diào)用LTM顯示IP核,將諧波參量顯示在LTM液晶屏上。
            2.如權(quán)利要求1所述的基于FPGA的電力諧波分析儀,其特征在于所述傳感器電路包括電壓信號采集電路和電流信號采集電路,電壓信號的采集電路的傳感器選用TV1013,運算放大器選用TI的0PA2227,電壓信號采集電路的輸出端連接采樣電路的模擬輸入端;電流信號采集電路的傳感器選用電流互感器TA1420,放大器選用0PA2227,電流信號采集電路的輸出端連接采樣電路的模擬輸入端。
            3.如權(quán)利要求1所述的基于FPGA的電力諧波分析儀,其特征在于所述AD采樣電路采用兩片AD7606-4并聯(lián),片外參考電壓方式,使用片外的基準(zhǔn)電壓REF5025作為AD7606的參考電壓,采樣電路AD7606與FPGA使用并行通信方式。
            4.如權(quán)利要求1所述的基于FPGA的電力諧波分析儀,其特征在于所述數(shù)據(jù)采集模塊驅(qū)動AD采樣電路進(jìn)行采樣是該數(shù)據(jù)采集模塊驅(qū)動AD7606進(jìn)行8路信號,1024點的同步采樣,然后將數(shù)據(jù)讀取到FPGA內(nèi),提供給下一級的數(shù)據(jù)緩存模塊,模塊通過對輸入時鐘計數(shù), 控制AD7606進(jìn)行采樣速率為10. 24K的等間隔采樣。
            5.如權(quán)利要求1所述的基于FPGA的電力諧波分析儀,其特征在于所述數(shù)據(jù)緩存模塊完成采樣數(shù)據(jù)的存儲、控制及輸出,是從ADC數(shù)據(jù)采樣模塊輸出的采樣數(shù)據(jù)按照1到8路數(shù)據(jù)一次輸出,系統(tǒng)需要將這8路數(shù)據(jù)分別存儲在片內(nèi)RAM中,當(dāng)有輸出觸發(fā)信號時,按設(shè)定順序依次輸出1路采樣數(shù)據(jù)。
            全文摘要
            一種基于FPGA的電力諧波分析儀,屬于電力電子技術(shù)領(lǐng)域,該電力諧波分析儀包括傳感器電路和AD采樣電路以及設(shè)在FPGA中的數(shù)據(jù)采集模塊、數(shù)據(jù)緩存模塊、數(shù)據(jù)類型轉(zhuǎn)換模塊、加窗模塊、FFT運算模塊、運算結(jié)果存儲和中斷產(chǎn)生模塊、NiosⅡ軟核處理器模塊、LTM顯示模塊。本發(fā)明前端的諧波采集和處理使用硬件實現(xiàn),充分發(fā)揮了硬件加速的優(yōu)勢;利用FPGA中的NiosⅡ軟核處理器對處理后的數(shù)據(jù)做進(jìn)一步運算和控制,該過程在FPGA片內(nèi)總線上完成,解決了數(shù)據(jù)的傳輸,充分發(fā)揮了硬件設(shè)計的高速性和NiosⅡ軟核處理器控制的靈活性。
            文檔編號G01R19/25GK102288819SQ201110236979
            公開日2011年12月21日 申請日期2011年8月18日 優(yōu)先權(quán)日2011年8月18日
            發(fā)明者李晶皎, 薛長江, 陸振林 申請人:東北大學(xué)
            網(wǎng)友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品