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

    基于PI 數(shù)據(jù)庫(kù)的燃驅(qū)壓縮機(jī)組性能測(cè)試系統(tǒng)

    2020-04-14 06:26:56
    天然氣工業(yè) 2020年11期
    關(guān)鍵詞:壓氣燃?xì)廨啓C(jī)壓縮機(jī)

    上海交通大學(xué)制冷與低溫工程研究所

    0 引言

    西氣東輸管道工程是目前中國(guó)距離最長(zhǎng)、管徑最大的輸氣干線[1]。在管道天然氣長(zhǎng)距離輸送過(guò)程中,為了使管道中天然氣的壓力保持在一個(gè)合理的范圍內(nèi),需要在管道天然氣輸送途中間隔設(shè)立許多個(gè)壓氣站,相鄰壓氣站之間的距離一般介于120~170 km[2]。壓氣站的關(guān)鍵核心設(shè)備是燃驅(qū)壓縮機(jī)組,目前國(guó)際上大型輸氣項(xiàng)目的壓縮機(jī)組多采用燃?xì)廨啓C(jī)驅(qū)動(dòng)。燃驅(qū)壓縮機(jī)組主要包括燃?xì)廨啓C(jī)裝置、離心壓縮機(jī)組以及相關(guān)的輔助系統(tǒng)。為了改善燃驅(qū)壓縮機(jī)組的性能、降低管道運(yùn)輸?shù)目傒斔统杀?,?duì)壓氣站中燃驅(qū)壓縮機(jī)組的實(shí)際工作情況進(jìn)行實(shí)時(shí)監(jiān)視和對(duì)機(jī)組一段時(shí)間的歷史運(yùn)行趨勢(shì)進(jìn)行分析顯得尤為重要。

    對(duì)此,國(guó)內(nèi)外學(xué)者已經(jīng)開(kāi)展了一系列的研究工作:雷劍宇等[3]開(kāi)發(fā)出了一種離心和軸流式鼓風(fēng)機(jī)以及壓縮機(jī)的熱力學(xué)性能測(cè)試系統(tǒng);Iranian 等[4]將數(shù)據(jù)即時(shí)傳輸?shù)交贚abVIEW 的上位機(jī)上,進(jìn)行數(shù)據(jù)分析和可視化;滿衛(wèi)芳等[5]實(shí)現(xiàn)了對(duì)燃?xì)廨啓C(jī)扭矩和轉(zhuǎn)速信號(hào)的采集傳輸、實(shí)時(shí)顯示等功能;李杰等[6]設(shè)計(jì)了對(duì)800 kW 離心壓縮機(jī)系統(tǒng)進(jìn)行動(dòng)態(tài)壓力測(cè)量的NI 測(cè)試系統(tǒng);郝巖[7]介紹了一種透平式鼓風(fēng)機(jī)以及壓縮機(jī)熱力性能試驗(yàn)的無(wú)線測(cè)試系統(tǒng);Varshini等[8]實(shí)現(xiàn)了采用雙反饋感應(yīng)發(fā)電機(jī)對(duì)風(fēng)力發(fā)電機(jī)的驅(qū)動(dòng);Kassas 等[9]基于LabVIEW 平臺(tái)對(duì)房屋內(nèi)的空調(diào)進(jìn)行功耗測(cè)量和調(diào)節(jié);徐艷雷[10]設(shè)計(jì)了一種基于LabVIEW 和Modbus 總線的燃?xì)廨啓C(jī)故障診斷系統(tǒng);賈廣博[11]采用模塊化的方法進(jìn)行了燃驅(qū)壓縮機(jī)組變工況特性的仿真;鄭濤[12]建立了基于MATLAB/Simulink平臺(tái)的燃驅(qū)壓縮機(jī)組各部件及總體仿真模型。

    國(guó)內(nèi)外的性能測(cè)試系統(tǒng)主要包括兩種:一種用MATLAB 和LabVIEW 等軟件進(jìn)行仿真模擬,為了模擬的準(zhǔn)確性,有些學(xué)者會(huì)利用實(shí)際數(shù)據(jù)進(jìn)行修正計(jì)算模型;另一種則是基于軟件開(kāi)發(fā)出一種應(yīng)用于現(xiàn)場(chǎng)的數(shù)據(jù)采集、分析計(jì)算和可視化的系統(tǒng)。筆者以30 MW國(guó)產(chǎn)燃驅(qū)壓縮機(jī)組為研究對(duì)象,基于MATLAB 平臺(tái)和Excel 軟件,建立了一套可應(yīng)用于遠(yuǎn)端監(jiān)測(cè)、實(shí)時(shí)顯示且可對(duì)機(jī)組歷史運(yùn)行情況進(jìn)行分析的性能測(cè)試系統(tǒng)。

    1 測(cè)試系統(tǒng)設(shè)計(jì)思路

    以西氣東輸工程某壓氣站中的30 MW 級(jí)燃驅(qū)壓縮機(jī)組為測(cè)試對(duì)象,讀取實(shí)時(shí)存儲(chǔ)在PI(Plant Information)數(shù)據(jù)庫(kù)里的數(shù)據(jù)到Excel 表格上,再將其數(shù)據(jù)實(shí)時(shí)傳輸?shù)組ATLAB 上,再基于MATLAB進(jìn)行計(jì)算和可視化,計(jì)算出燃驅(qū)壓縮機(jī)組以及組件的性能參數(shù),達(dá)到實(shí)時(shí)計(jì)算、實(shí)時(shí)顯示和實(shí)時(shí)監(jiān)控的目的。測(cè)試系統(tǒng)設(shè)計(jì)思路原理如圖1 所示。

    圖1 測(cè)試系統(tǒng)設(shè)計(jì)思路原理圖

    PI 數(shù)據(jù)庫(kù)是由美國(guó)OSI 公司開(kāi)發(fā)的一套實(shí)時(shí)(歷史)數(shù)據(jù)庫(kù)系統(tǒng),可在線存儲(chǔ)多個(gè)工藝點(diǎn)的多年歷史數(shù)據(jù),并可對(duì)其進(jìn)行實(shí)時(shí)訪問(wèn)和讀取?;贛ATLAB 開(kāi)發(fā)的測(cè)試程序數(shù)據(jù)處理速度快、效率高、操作方便、界面美觀、功能擴(kuò)展性好[13]。PI_DATALINK 是用于連接PI 數(shù)據(jù)庫(kù)和Excel 的一種軟件[14]。EXCEL_LINK 也是一種可用于Excel 加載的插件,加載完成后,可以直接在MATLAB 中調(diào)用存取Excel 中的數(shù)據(jù)。

    該壓氣站中燃驅(qū)壓縮機(jī)組測(cè)量點(diǎn)的數(shù)據(jù)都存儲(chǔ)在PI 數(shù)據(jù)庫(kù)中,遠(yuǎn)程監(jiān)控中心可直接通過(guò)PI 數(shù)據(jù)庫(kù)查詢機(jī)組的基本參數(shù)。但是想通過(guò)外部接口快速簡(jiǎn)單地實(shí)時(shí)讀取數(shù)據(jù)庫(kù)中的數(shù)據(jù)并不容易。為此,在綜合考察幾種常見(jiàn)讀取PI 數(shù)據(jù)庫(kù)接口方式的基礎(chǔ)上,選取了一種簡(jiǎn)單的實(shí)時(shí)讀取方式。首先通過(guò)插件PI_DATALINK 用Excel 表格實(shí)時(shí)讀取PI 數(shù)據(jù)庫(kù)中的數(shù)據(jù),再通過(guò)插件EXCEL_LINK 用MATLAB 同步讀取Excel 中的數(shù)據(jù),然后在MATLAB 中進(jìn)行計(jì)算和可視化。

    2 測(cè)試系統(tǒng)模塊設(shè)計(jì)與實(shí)現(xiàn)

    筆者設(shè)計(jì)的燃驅(qū)壓縮機(jī)組性能系統(tǒng)程序主要基于兩大軟件,即Excel 和MATLAB。通過(guò)Excel 實(shí)時(shí)讀取PI 數(shù)據(jù)庫(kù)中的數(shù)據(jù),并存儲(chǔ)在表格中,然后利用MATLAB 直接在Excel 中讀取數(shù)據(jù),將讀取的數(shù)據(jù)代入到建立的計(jì)算模型中進(jìn)行計(jì)算并進(jìn)行可視化、存儲(chǔ)以及歷史數(shù)據(jù)的分析。

    2.1 Excel 模塊

    在Excel 模塊中實(shí)現(xiàn)的功能為數(shù)據(jù)庫(kù)數(shù)據(jù)的實(shí)時(shí)讀取以及將數(shù)據(jù)實(shí)時(shí)傳輸?shù)組ATLAB 環(huán)境中。因此,Excel 模塊又可分為數(shù)據(jù)實(shí)時(shí)讀取模塊和數(shù)據(jù)實(shí)時(shí)輸送模塊。

    2.1.1 數(shù)據(jù)實(shí)時(shí)讀取模塊

    通過(guò)在Excel 加載PI_DATALINK 插件后,可以通過(guò)Excel 工具欄中出現(xiàn)的“PI”選項(xiàng)進(jìn)行讀取,但是無(wú)法實(shí)現(xiàn)實(shí)時(shí)讀取。因此采用Excel VBA[15]進(jìn)行編程,通過(guò)調(diào)用函數(shù)的方法達(dá)到實(shí)時(shí)讀取的目標(biāo)。

    所設(shè)計(jì)的讀取模塊的主要功能就是實(shí)時(shí)讀取PI數(shù)據(jù)庫(kù)的數(shù)據(jù),并將一段時(shí)間的數(shù)據(jù)存儲(chǔ)在這張表格中。功能實(shí)現(xiàn)思路是將表格中的一個(gè)單元格定義為當(dāng)前時(shí)間,另一個(gè)單元格定義為與當(dāng)前時(shí)間有一定時(shí)間間隔的時(shí)刻,讓通過(guò)特定調(diào)用函數(shù)從數(shù)據(jù)庫(kù)調(diào)取出來(lái)的數(shù)據(jù)以這兩個(gè)時(shí)間單元格的時(shí)間為起始時(shí)間和終止時(shí)間,按照特定的時(shí)間步長(zhǎng)實(shí)時(shí)存儲(chǔ)。每一列顯示的是一個(gè)參數(shù)在不同的時(shí)間點(diǎn)上的數(shù)據(jù),每一行顯示的是一個(gè)時(shí)刻不同參數(shù)的數(shù)據(jù)。

    2.1.2 數(shù)據(jù)實(shí)時(shí)輸送模塊

    在Excel 加載EXCEL_LINK 插件后,通過(guò)調(diào)用名為MLPutMatrix 的函數(shù),任選一個(gè)空白單元格,輸入命令“=MLPutMatrix(name,mdat)”,其中MLPutMatrix 為Excel 將數(shù)據(jù)以矩陣的形式直接傳輸?shù)組ATLAB 上的內(nèi)部函數(shù),“name”為在MATLAB中該數(shù)據(jù)矩陣的名稱,“mdat”為Excel 傳輸?shù)組ATLAB 的數(shù)據(jù)區(qū)域。通過(guò)該函數(shù),可直接將Excel特定區(qū)域的數(shù)據(jù)以矩陣的形式輸送到MATLAB 中,并且能夠在MATLAB 環(huán)境中與Excel 表格中的數(shù)據(jù)同步變化。

    2.2 MATLAB 模塊

    根據(jù)功能的不同,MATLAB 模塊中又分為天然氣組分含量輸入模塊、機(jī)組性能參數(shù)計(jì)算模塊、機(jī)組性能數(shù)據(jù)存儲(chǔ)模塊、機(jī)組性能可視化模塊以及機(jī)組性能歷史趨勢(shì)分析模塊。

    2.2.1 天然氣組分含量輸入模塊

    計(jì)算燃驅(qū)壓縮機(jī)組的功率、效率等性能參數(shù)需要已知燃?xì)廨啓C(jī)進(jìn)口管道天然氣的組分及其含量,而壓氣站的PI 數(shù)據(jù)庫(kù)里沒(méi)有存儲(chǔ)這些數(shù)據(jù),所以需要建立一個(gè)手動(dòng)輸入模塊。

    2.2.2 機(jī)組性能參數(shù)計(jì)算模塊

    機(jī)組性能參數(shù)計(jì)算模塊根據(jù)組件可分為燃?xì)廨啓C(jī)性能計(jì)算模塊、離心壓縮機(jī)機(jī)組性能計(jì)算模塊、機(jī)組總效率計(jì)算模塊。因?yàn)橛?jì)算離心壓縮機(jī)性能參數(shù)時(shí),進(jìn)口天然氣的壓縮因子是個(gè)重要參數(shù),且實(shí)際氣體的天然氣的壓縮因子的求解較為復(fù)雜,故在此基礎(chǔ)上增加一個(gè)天然氣的壓縮因子的計(jì)算模塊。

    2.2.2.1 燃?xì)廨啓C(jī)性能計(jì)算莫模塊

    燃?xì)廨啓C(jī)的輸出軸功率利用熱平衡的方法求得。假定一個(gè)包括燃?xì)廨啓C(jī)在內(nèi)的合適的“控制體”應(yīng)該便于確定進(jìn)入和離開(kāi)“控制體”的所有熱量和能量的參數(shù),輸出功率是這些參數(shù)中的一項(xiàng),利用熱平衡方程式和其他已知數(shù)據(jù)可計(jì)算得出。假定的燃?xì)廨啓C(jī)熱平衡“控制體”如圖2 所示。

    圖2 燃?xì)廨啓C(jī)熱平衡“控制體”示意圖

    標(biāo)準(zhǔn)條件下壓力為101.325 kPa,溫度為15 ℃??刂企w的熱平衡方程如下:

    式中ma表示進(jìn)入壓氣機(jī)的空氣流量,kg/s;mf表示進(jìn)入控制體的燃料供給量,kg/s;mg表示出口煙氣流量,kg/s;me表示氣封等泄漏和抽氣的流量,kg/s:ha表示進(jìn)入壓氣機(jī)的空氣比焓,kJ/kg;hg表示出口煙氣的比焓,kJ/kg;he表示泄漏空氣的比焓,kJ/kg;ha0表示在標(biāo)準(zhǔn)條件下空氣的比焓,kJ/kg;hf表示進(jìn)口燃?xì)獾谋褥?,kJ/kg;hf0、hg0分別表示在標(biāo)準(zhǔn)條件下燃?xì)狻⒖諝獾谋褥?,kJ/kg;QL表示各種熱損失(約為燃?xì)廨啓C(jī)熱耗量的30%),kW;PS表示輸出軸供給負(fù)載的負(fù)荷,kW;ηtc表示燃燒效率,一般取0.99;LHV表示燃料在標(biāo)準(zhǔn)條件下的低位熱值,kJ/kg。

    出口煙氣流量與進(jìn)口空氣流量、燃料氣流量以及泄漏氣流量的關(guān)系如下:

    壓氣站現(xiàn)場(chǎng)測(cè)量得到的泄漏量的總和約為壓氣機(jī)進(jìn)口空氣量的1%。空氣、燃?xì)夂蜔煔獾谋褥?,可以根?jù)從PI 數(shù)據(jù)庫(kù)讀取的溫度進(jìn)行插值得到。

    根據(jù)式(1)和(2)計(jì)算得出輸出軸功率(PS)時(shí),由于系統(tǒng)熱損失、泄漏量等數(shù)據(jù)和實(shí)際條件有關(guān),并不是一個(gè)恒定值,而且無(wú)法從數(shù)據(jù)庫(kù)中實(shí)時(shí)獲取,為了讓計(jì)算模型的結(jié)果更加接近實(shí)際情況,以一個(gè)時(shí)間段中不同工況下實(shí)際測(cè)量出來(lái)的輸出軸功率為修正標(biāo)準(zhǔn),采用分段插值逼近的方法,通過(guò)給定不同工況下的修正系數(shù)來(lái)修正這些數(shù)據(jù)。表1 為另外的運(yùn)行時(shí)間下不同工況的輸出軸功率。

    表1 不同運(yùn)行工況下的輸出軸功率表

    燃?xì)廨啓C(jī)熱耗量的計(jì)算公式如下:

    式中HI表示單位時(shí)間熱耗量,kJ/s。

    燃?xì)廨啓C(jī)熱效率(對(duì)應(yīng)于凈機(jī)械功率的熱效率)計(jì)算如下:

    式中ηt表示熱效率。

    燃?xì)廨啓C(jī)熱耗率計(jì)算公式如下:

    式中HR表示熱耗率,kJ/kWh。

    2.2.2.2 離心壓縮機(jī)組性能計(jì)算模塊

    壓縮機(jī)實(shí)際熱力過(guò)程中,不可避免地存在能量損耗。因此,實(shí)際過(guò)程均為不可逆過(guò)程,為了便于分析,此處將實(shí)際過(guò)程視為可逆多變過(guò)程計(jì)算其壓縮功和多變效率。

    離心壓縮機(jī)的多變系數(shù)計(jì)算如下:

    式中m表示多變指數(shù);T1、T2分別表示壓縮機(jī)入口、出口溫度,K;p1、p2分別表示壓縮機(jī)入口、出口壓力,Pa。

    壓縮機(jī)的多變壓縮功計(jì)算如下:

    式中wpol表示壓縮機(jī)多變壓縮功,J/kg;Z 表示實(shí)際氣體壓縮因子;Rg表示氣體常數(shù),J/(kg·K),Rg=8 314.41/M。

    多變壓縮功率計(jì)算如下:

    式中Ppol表示多變壓縮功率,kW;ρ表示標(biāo)準(zhǔn)條件下,進(jìn)入壓縮機(jī)的天然氣的密度,kg/m3;qv表示進(jìn)入壓縮機(jī)的氣體流量,m3/h。

    壓縮機(jī)的多變效率計(jì)算如下:

    式中ηpol表示多變效率;κ表示絕熱指數(shù),和氣體組分有關(guān)。這里取值范圍為1.25~1.32。

    2.2.2.3 機(jī)組總效率的計(jì)算模塊

    燃驅(qū)壓縮機(jī)組的總效率計(jì)算如下:

    式中ηtot表示機(jī)組總效率。

    2.2.2.4 天然氣壓縮因子的計(jì)算模塊

    天然氣壓縮因子的計(jì)算采用AGA8-92DC 方程[16],其應(yīng)用范圍主要為管道運(yùn)輸氣體,適合的壓力范圍為0~12 MPa,溫度范圍為263~338 K[17]。其計(jì)算公式如下:

    式中B表示第二維利系數(shù),m3/kmol;ρm表示摩爾密度,kmol/m3;ρr表示對(duì)比密度;bn、cn、kn表示常數(shù),其值可由文獻(xiàn)[17]直接查得;Cn*表示與溫度和組分有關(guān)的系數(shù)。ρm、ρr、B和Cn*也可由本文參考文獻(xiàn)[17]相應(yīng)的公式計(jì)算得到。

    2.2.3 機(jī)組性能數(shù)據(jù)存儲(chǔ)模塊

    數(shù)據(jù)存儲(chǔ)模塊的功能是把計(jì)算出來(lái)的燃驅(qū)壓縮機(jī)組的各項(xiàng)性能參數(shù)按照時(shí)間排序,以表格的形式存儲(chǔ)在本地文件中。表格每一行為同一時(shí)刻燃驅(qū)壓縮機(jī)組的各項(xiàng)性能參數(shù)的計(jì)算值,每一列為同一指標(biāo)在不同時(shí)刻的計(jì)算值,時(shí)間間隔為1 s,表格的名稱直接以當(dāng)前時(shí)間命名,具體形式為“年_月_日_時(shí)_分_秒”,便于查找和讀取。

    2.2.4 機(jī)組性能可視化模塊

    性能可視化模塊的功能是實(shí)時(shí)顯示燃驅(qū)壓縮機(jī)基本參數(shù)的數(shù)值以及一段時(shí)間內(nèi)各項(xiàng)性能指標(biāo)的數(shù)值以及趨勢(shì)圖,其主界面如圖3 所示。圖3 上半部分實(shí)時(shí)顯示燃驅(qū)壓縮機(jī)組的流量等現(xiàn)場(chǎng)測(cè)量的基本數(shù)值,一直隨時(shí)間變動(dòng),為當(dāng)前時(shí)間的準(zhǔn)確數(shù)值,下半部分以折線圖的形式顯示燃?xì)廨啓C(jī)、壓縮機(jī)和機(jī)組的一段時(shí)間內(nèi)性能參數(shù)隨著動(dòng)力渦輪轉(zhuǎn)速的變化趨勢(shì)圖。因此,從顯示界面中可直觀看出燃驅(qū)壓縮機(jī)組的運(yùn)行狀態(tài),從而達(dá)到實(shí)時(shí)監(jiān)控的目的。

    圖3 燃驅(qū)壓縮機(jī)組性能指標(biāo)顯示主界面圖

    2.2.5 機(jī)組性能歷史趨勢(shì)分析模塊

    歷史趨勢(shì)分析模塊的功能是將一段時(shí)間內(nèi)燃驅(qū)壓縮機(jī)組性能指標(biāo)參數(shù)的計(jì)算數(shù)值以趨勢(shì)圖的形式展現(xiàn)出來(lái)。歷史趨勢(shì)分析模塊的主要界面如圖4 所示。通過(guò)選定需要分析的時(shí)間段、機(jī)組名稱和具體性能指標(biāo),該項(xiàng)指標(biāo)會(huì)以折線圖的形式,以時(shí)間為橫坐標(biāo)顯示出其該段時(shí)間內(nèi)的趨勢(shì)分布,以一種非常直觀的形式分析這段時(shí)間內(nèi)燃驅(qū)壓縮機(jī)組運(yùn)行的實(shí)際情況,便于機(jī)組的改進(jìn)和優(yōu)化。

    圖4 燃驅(qū)壓縮機(jī)組性能歷史趨勢(shì)分析主界面圖

    3 結(jié)論

    1)通過(guò)Excel 作為中間接口軟件,用一種較為簡(jiǎn)單的方法實(shí)現(xiàn)了PI 數(shù)據(jù)庫(kù)的數(shù)據(jù)讀取,并通過(guò)Excel VBA 實(shí)現(xiàn)了對(duì)數(shù)據(jù)庫(kù)實(shí)時(shí)讀取的功能;通過(guò)Excel 內(nèi)部函數(shù)實(shí)現(xiàn)了向MATLAB 實(shí)時(shí)傳遞數(shù)據(jù)的功能。

    2)數(shù)據(jù)的實(shí)時(shí)傳輸保存以及實(shí)時(shí)顯示皆正常,可視化監(jiān)控界面直觀,人機(jī)交互性好,軟件操作簡(jiǎn)單方便,可遠(yuǎn)程監(jiān)控燃驅(qū)壓縮機(jī)組的實(shí)時(shí)運(yùn)行狀態(tài)。

    3)設(shè)計(jì)的燃驅(qū)壓縮機(jī)組性能測(cè)試程序已經(jīng)應(yīng)用于西氣東輸管道工程中,并根據(jù)現(xiàn)場(chǎng)數(shù)據(jù)對(duì)性能計(jì)算模型進(jìn)行了修正,計(jì)算出來(lái)的燃驅(qū)壓縮機(jī)組的性能指標(biāo)數(shù)值誤差范圍滿足實(shí)際壓氣站的要求。

    4)建立了機(jī)組性能歷史數(shù)據(jù)分析趨勢(shì)模塊,為管道運(yùn)輸長(zhǎng)距離運(yùn)輸過(guò)程中燃驅(qū)壓縮機(jī)組的優(yōu)化和改進(jìn)提供參考。

    猜你喜歡
    壓氣燃?xì)廨啓C(jī)壓縮機(jī)
    小棚蝦“壓氣頭”現(xiàn)象分析及探究
    長(zhǎng)輸天然氣管道壓氣站運(yùn)行分析
    渦輪增壓器壓氣葉輪爆裂轉(zhuǎn)速數(shù)值分析與試驗(yàn)研究
    BOG壓縮機(jī)在小型LNG船舶上的應(yīng)用
    淺析特大型往復(fù)壓縮機(jī)組的安裝
    壓縮機(jī)組智能控制與節(jié)能增效
    長(zhǎng)輸天然氣壓氣站進(jìn)出站管道應(yīng)力分析
    《燃?xì)廨啓C(jī)技術(shù)》2014年索引
    SGT5-4000F(4)燃?xì)廨啓C(jī)夏季最大負(fù)荷研究及應(yīng)用
    輕型燃?xì)廨啓C(jī)LM6000PC與重型燃?xì)廨啓C(jī)PG6581B研究與對(duì)比分析
    国产精品免费视频内射| 日韩欧美国产一区二区入口| 丰满的人妻完整版| 色播亚洲综合网| 黄色片一级片一级黄色片| 日韩大尺度精品在线看网址| 亚洲成人久久性| 亚洲国产看品久久| 亚洲av电影在线进入| 欧美色视频一区免费| 99久久综合精品五月天人人| 亚洲专区字幕在线| 亚洲国产精品sss在线观看| 午夜福利免费观看在线| 久久午夜亚洲精品久久| 精品第一国产精品| 女同久久另类99精品国产91| 宅男免费午夜| 国产成人精品无人区| 欧美3d第一页| 亚洲国产精品合色在线| 99re在线观看精品视频| 色精品久久人妻99蜜桃| 成年免费大片在线观看| 亚洲自偷自拍图片 自拍| 久久久久亚洲av毛片大全| 亚洲中文av在线| 最近在线观看免费完整版| 黑人欧美特级aaaaaa片| 悠悠久久av| 国产精品,欧美在线| 最近最新免费中文字幕在线| 亚洲欧洲精品一区二区精品久久久| 又黄又爽又免费观看的视频| 丰满的人妻完整版| 十八禁人妻一区二区| 成人18禁在线播放| 夜夜看夜夜爽夜夜摸| 日韩欧美一区二区三区在线观看| 久久久久九九精品影院| 精品国产乱子伦一区二区三区| 亚洲中文字幕日韩| 亚洲熟妇熟女久久| 亚洲人成伊人成综合网2020| 熟妇人妻久久中文字幕3abv| 男人舔女人下体高潮全视频| 久久久国产欧美日韩av| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| 久久久久精品国产欧美久久久| 国产黄片美女视频| 丰满人妻一区二区三区视频av | 一级毛片高清免费大全| 男女午夜视频在线观看| 亚洲成人精品中文字幕电影| 亚洲七黄色美女视频| 免费看日本二区| 国产伦在线观看视频一区| 国产三级在线视频| 国产成年人精品一区二区| 色综合站精品国产| x7x7x7水蜜桃| 丰满人妻熟妇乱又伦精品不卡| 天天一区二区日本电影三级| 99久久精品热视频| 在线观看免费日韩欧美大片| 夜夜躁狠狠躁天天躁| 一边摸一边做爽爽视频免费| 波多野结衣高清无吗| 婷婷精品国产亚洲av在线| av天堂在线播放| 色噜噜av男人的天堂激情| 色综合亚洲欧美另类图片| 国产主播在线观看一区二区| 国产真人三级小视频在线观看| 国产精品久久久久久精品电影| 国产精品自产拍在线观看55亚洲| 又爽又黄无遮挡网站| 亚洲欧美日韩东京热| 美女免费视频网站| 十八禁人妻一区二区| 精品久久久久久久久久免费视频| 国产爱豆传媒在线观看 | 一本综合久久免费| 婷婷六月久久综合丁香| 亚洲成av人片免费观看| 亚洲成av人片在线播放无| 动漫黄色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日本一区二区免费在线视频| 国内精品一区二区在线观看| 免费电影在线观看免费观看| 免费在线观看成人毛片| 国产精品1区2区在线观看.| 国产真人三级小视频在线观看| 99久久无色码亚洲精品果冻| 国内精品久久久久精免费| 免费在线观看影片大全网站| 丝袜美腿诱惑在线| 黄色丝袜av网址大全| 熟女电影av网| 亚洲中文日韩欧美视频| 国产成+人综合+亚洲专区| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机午夜福利在线观看视频| 亚洲国产精品久久男人天堂| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av日韩精品久久久久久密| 精品一区二区三区四区五区乱码| 国内精品一区二区在线观看| 精品乱码久久久久久99久播| 很黄的视频免费| 国产精品亚洲美女久久久| 国产精品久久久av美女十八| 一边摸一边抽搐一进一小说| 成人欧美大片| 18禁黄网站禁片免费观看直播| 久久人妻福利社区极品人妻图片| av超薄肉色丝袜交足视频| 最近最新免费中文字幕在线| 午夜激情福利司机影院| 精品久久蜜臀av无| 日本黄色视频三级网站网址| 90打野战视频偷拍视频| 美女扒开内裤让男人捅视频| 不卡一级毛片| 俄罗斯特黄特色一大片| 午夜久久久久精精品| 手机成人av网站| 国产成人系列免费观看| 久久性视频一级片| 一本综合久久免费| 少妇被粗大的猛进出69影院| 狂野欧美白嫩少妇大欣赏| 国产片内射在线| 成熟少妇高潮喷水视频| 久99久视频精品免费| 亚洲avbb在线观看| 久久人妻av系列| 中文字幕av在线有码专区| 手机成人av网站| 老熟妇仑乱视频hdxx| 国产高清激情床上av| 色综合亚洲欧美另类图片| 露出奶头的视频| 毛片女人毛片| 色av中文字幕| 国产亚洲精品综合一区在线观看 | 国产一区在线观看成人免费| 亚洲男人天堂网一区| 亚洲中文av在线| 亚洲一卡2卡3卡4卡5卡精品中文| www.熟女人妻精品国产| 国产伦在线观看视频一区| 亚洲熟女毛片儿| 亚洲国产精品sss在线观看| 亚洲精品中文字幕一二三四区| 伦理电影免费视频| 老司机靠b影院| 在线国产一区二区在线| 亚洲精品av麻豆狂野| 欧美成人午夜精品| 精品久久久久久久毛片微露脸| 少妇的丰满在线观看| 国产精品免费视频内射| 亚洲性夜色夜夜综合| 亚洲av电影不卡..在线观看| 18禁观看日本| 亚洲精华国产精华精| 亚洲在线自拍视频| 18禁观看日本| 日韩高清综合在线| 在线观看美女被高潮喷水网站 | 色av中文字幕| 搡老熟女国产l中国老女人| 亚洲精品中文字幕在线视频| 午夜福利18| 亚洲精品色激情综合| 黄片大片在线免费观看| 老司机午夜十八禁免费视频| 伦理电影免费视频| 天天一区二区日本电影三级| 久久欧美精品欧美久久欧美| 黄色视频不卡| 免费看十八禁软件| 亚洲精品中文字幕一二三四区| a级毛片a级免费在线| 久久久久九九精品影院| 老司机午夜十八禁免费视频| 日日爽夜夜爽网站| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站| 18美女黄网站色大片免费观看| 亚洲人与动物交配视频| 中文字幕最新亚洲高清| 婷婷精品国产亚洲av| 日本 av在线| 九九热线精品视视频播放| 国产欧美日韩一区二区精品| 欧美成人免费av一区二区三区| 国产精品久久电影中文字幕| 国产成人av教育| 国产精品野战在线观看| 岛国在线观看网站| 国产三级在线视频| 国产精品九九99| 哪里可以看免费的av片| 亚洲国产欧美一区二区综合| 国产成人欧美在线观看| 免费一级毛片在线播放高清视频| 成人特级黄色片久久久久久久| 久久精品国产亚洲av香蕉五月| or卡值多少钱| 香蕉国产在线看| 亚洲第一电影网av| 精品福利观看| 日韩 欧美 亚洲 中文字幕| cao死你这个sao货| 真人做人爱边吃奶动态| 丰满人妻一区二区三区视频av | 久久精品国产亚洲av香蕉五月| 免费在线观看黄色视频的| 国产精品一区二区三区四区免费观看 | 精品国产超薄肉色丝袜足j| 啦啦啦免费观看视频1| 亚洲狠狠婷婷综合久久图片| 久99久视频精品免费| 国产精品野战在线观看| 久久久久久人人人人人| 久久99热这里只有精品18| 又大又爽又粗| 欧美绝顶高潮抽搐喷水| 久久精品综合一区二区三区| 两个人看的免费小视频| 在线视频色国产色| 热99re8久久精品国产| av免费在线观看网站| 午夜福利免费观看在线| 99精品在免费线老司机午夜| 久久天堂一区二区三区四区| 亚洲欧美精品综合久久99| 999精品在线视频| 成人午夜高清在线视频| 国产亚洲精品av在线| 在线免费观看的www视频| 国产v大片淫在线免费观看| 国产一区二区三区视频了| cao死你这个sao货| 免费在线观看视频国产中文字幕亚洲| 国产精品日韩av在线免费观看| 国产成人啪精品午夜网站| 亚洲av成人不卡在线观看播放网| 亚洲精品久久国产高清桃花| 国产精品一区二区免费欧美| 在线观看免费视频日本深夜| 午夜福利18| 一个人免费在线观看的高清视频| 欧美色视频一区免费| 日日爽夜夜爽网站| 一边摸一边做爽爽视频免费| 他把我摸到了高潮在线观看| 国模一区二区三区四区视频 | 久久久久久久久免费视频了| 最好的美女福利视频网| 在线观看美女被高潮喷水网站 | 亚洲欧美一区二区三区黑人| e午夜精品久久久久久久| 最新美女视频免费是黄的| 国产午夜福利久久久久久| av福利片在线观看| 亚洲天堂国产精品一区在线| 日韩国内少妇激情av| 亚洲精品在线观看二区| 草草在线视频免费看| 国产亚洲精品第一综合不卡| 麻豆一二三区av精品| 女人高潮潮喷娇喘18禁视频| 变态另类丝袜制服| 免费在线观看成人毛片| 欧美高清成人免费视频www| www.自偷自拍.com| 成在线人永久免费视频| 又爽又黄无遮挡网站| 精品福利观看| 天堂动漫精品| 国产成人一区二区三区免费视频网站| 妹子高潮喷水视频| 午夜免费激情av| 手机成人av网站| 久久久水蜜桃国产精品网| 亚洲成av人片免费观看| 亚洲人成网站高清观看| 国产不卡一卡二| 亚洲免费av在线视频| 每晚都被弄得嗷嗷叫到高潮| 91av网站免费观看| 91大片在线观看| 少妇粗大呻吟视频| 国产精品av久久久久免费| 国内精品久久久久精免费| 日韩欧美一区二区三区在线观看| 久久伊人香网站| 日韩欧美 国产精品| 精品乱码久久久久久99久播| 99精品久久久久人妻精品| 看黄色毛片网站| 丁香欧美五月| 最近最新免费中文字幕在线| 黄色丝袜av网址大全| 亚洲成人精品中文字幕电影| 久久久久久亚洲精品国产蜜桃av| 日韩欧美 国产精品| 欧美在线一区亚洲| 成人永久免费在线观看视频| 我要搜黄色片| 欧美绝顶高潮抽搐喷水| 日本a在线网址| 国产成人av激情在线播放| 999精品在线视频| 国产精品九九99| www.精华液| 人人妻,人人澡人人爽秒播| 免费搜索国产男女视频| 国产av一区在线观看免费| 国产一级毛片七仙女欲春2| 香蕉丝袜av| 人人妻人人看人人澡| 国产一区二区在线av高清观看| 亚洲av成人精品一区久久| 免费av毛片视频| 少妇的丰满在线观看| 动漫黄色视频在线观看| 国产一区二区在线av高清观看| 亚洲欧美一区二区三区黑人| 亚洲一区高清亚洲精品| 国内精品一区二区在线观看| 亚洲国产欧洲综合997久久,| 两个人看的免费小视频| 国产私拍福利视频在线观看| 欧美大码av| 日本 av在线| 色av中文字幕| 草草在线视频免费看| 成人亚洲精品av一区二区| 午夜影院日韩av| e午夜精品久久久久久久| 波多野结衣高清无吗| 黄色片一级片一级黄色片| 国产爱豆传媒在线观看 | 国语自产精品视频在线第100页| 成人午夜高清在线视频| 国语自产精品视频在线第100页| 国产成人精品久久二区二区91| 欧美绝顶高潮抽搐喷水| 一个人免费在线观看电影 | 亚洲自拍偷在线| 日本撒尿小便嘘嘘汇集6| 天天躁夜夜躁狠狠躁躁| 亚洲欧美激情综合另类| 国产av在哪里看| 很黄的视频免费| 波多野结衣高清无吗| 久久久精品国产亚洲av高清涩受| 人人妻人人澡欧美一区二区| 可以免费在线观看a视频的电影网站| www.熟女人妻精品国产| 在线视频色国产色| 最新在线观看一区二区三区| 啦啦啦免费观看视频1| 婷婷丁香在线五月| 男女做爰动态图高潮gif福利片| www国产在线视频色| 美女大奶头视频| 很黄的视频免费| 夜夜看夜夜爽夜夜摸| 90打野战视频偷拍视频| 亚洲国产精品久久男人天堂| 久久天堂一区二区三区四区| 观看免费一级毛片| 国产av不卡久久| 看免费av毛片| 琪琪午夜伦伦电影理论片6080| 麻豆一二三区av精品| 亚洲熟妇熟女久久| 国产精品电影一区二区三区| 国产av在哪里看| 久久亚洲真实| 亚洲欧美日韩高清在线视频| 18禁黄网站禁片免费观看直播| aaaaa片日本免费| 亚洲精品在线观看二区| 可以免费在线观看a视频的电影网站| 国产在线观看jvid| 国内精品久久久久精免费| 欧美午夜高清在线| 国产精品永久免费网站| 国产私拍福利视频在线观看| 免费观看精品视频网站| 久久久久久大精品| 老司机深夜福利视频在线观看| 我要搜黄色片| 国产黄片美女视频| 国产精华一区二区三区| 999久久久国产精品视频| 非洲黑人性xxxx精品又粗又长| av天堂在线播放| 女人被狂操c到高潮| 巨乳人妻的诱惑在线观看| 老司机福利观看| 在线观看一区二区三区| 最近最新中文字幕大全电影3| 男女床上黄色一级片免费看| 中文字幕av在线有码专区| 久久中文字幕人妻熟女| 一本综合久久免费| 母亲3免费完整高清在线观看| av免费在线观看网站| 高清毛片免费观看视频网站| 免费在线观看影片大全网站| 久久久久性生活片| 亚洲乱码一区二区免费版| 久久久国产成人精品二区| 99re在线观看精品视频| 国产又黄又爽又无遮挡在线| 久久精品影院6| 精品久久久久久久人妻蜜臀av| 90打野战视频偷拍视频| 亚洲欧美日韩东京热| a级毛片在线看网站| 在线观看免费视频日本深夜| 国产真实乱freesex| 亚洲中文字幕一区二区三区有码在线看 | 此物有八面人人有两片| 久久精品国产亚洲av高清一级| 国产亚洲精品综合一区在线观看 | 国产激情偷乱视频一区二区| 国产一区二区在线av高清观看| 99久久99久久久精品蜜桃| 一夜夜www| 黄色成人免费大全| 香蕉国产在线看| 欧美+亚洲+日韩+国产| 成年免费大片在线观看| 两个人免费观看高清视频| 九色国产91popny在线| АⅤ资源中文在线天堂| 亚洲一区二区三区不卡视频| 人成视频在线观看免费观看| 妹子高潮喷水视频| 日本撒尿小便嘘嘘汇集6| 美女午夜性视频免费| 法律面前人人平等表现在哪些方面| 国产精品免费一区二区三区在线| 50天的宝宝边吃奶边哭怎么回事| 男女下面进入的视频免费午夜| 99久久国产精品久久久| 一二三四在线观看免费中文在| 国产午夜精品久久久久久| 色尼玛亚洲综合影院| 国产免费av片在线观看野外av| 成年免费大片在线观看| 一二三四在线观看免费中文在| 国产高清有码在线观看视频 | 国产亚洲精品久久久久久毛片| 国产精品久久视频播放| 亚洲一区中文字幕在线| 久久中文字幕一级| 中文字幕av在线有码专区| 欧美午夜高清在线| 日本 欧美在线| 欧美日韩瑟瑟在线播放| 51午夜福利影视在线观看| 久久精品亚洲精品国产色婷小说| 一夜夜www| 欧美日韩精品网址| 亚洲五月天丁香| 熟女电影av网| 亚洲色图 男人天堂 中文字幕| 一进一出抽搐gif免费好疼| 久99久视频精品免费| 日韩精品青青久久久久久| 亚洲免费av在线视频| 巨乳人妻的诱惑在线观看| 国产成人av教育| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| 在线a可以看的网站| 欧美成人性av电影在线观看| 91在线观看av| 欧美一区二区国产精品久久精品 | 午夜免费成人在线视频| 欧洲精品卡2卡3卡4卡5卡区| 老汉色av国产亚洲站长工具| 国产成+人综合+亚洲专区| 久久精品成人免费网站| 97碰自拍视频| 成人18禁在线播放| 欧美日韩中文字幕国产精品一区二区三区| 成人国产综合亚洲| 久9热在线精品视频| 精品人妻1区二区| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三| 99国产精品99久久久久| 老鸭窝网址在线观看| 老汉色∧v一级毛片| 国内精品久久久久久久电影| 久久99热这里只有精品18| 国产激情偷乱视频一区二区| 国产精品久久久久久人妻精品电影| 人妻夜夜爽99麻豆av| 国产区一区二久久| 亚洲成人国产一区在线观看| 欧美乱码精品一区二区三区| 一边摸一边抽搐一进一小说| 亚洲成人久久性| 婷婷丁香在线五月| 亚洲精品美女久久av网站| 熟妇人妻久久中文字幕3abv| 给我免费播放毛片高清在线观看| 女人高潮潮喷娇喘18禁视频| 午夜久久久久精精品| 国产一区二区三区在线臀色熟女| 日本成人三级电影网站| 久久久久性生活片| 1024香蕉在线观看| 欧美日本视频| 男人的好看免费观看在线视频 | √禁漫天堂资源中文www| 国产黄a三级三级三级人| 久久久久国产精品人妻aⅴ院| 亚洲成av人片免费观看| 一级a爱片免费观看的视频| 91老司机精品| 欧美成人性av电影在线观看| 久久久久久人人人人人| 欧美中文综合在线视频| videosex国产| 麻豆成人av在线观看| aaaaa片日本免费| 亚洲成人国产一区在线观看| 国产成+人综合+亚洲专区| 给我免费播放毛片高清在线观看| 国产精品久久久av美女十八| 中文字幕久久专区| 岛国在线免费视频观看| 每晚都被弄得嗷嗷叫到高潮| 美女扒开内裤让男人捅视频| 久久久国产精品麻豆| 在线a可以看的网站| 99久久精品国产亚洲精品| 丰满的人妻完整版| 女人爽到高潮嗷嗷叫在线视频| 国产精华一区二区三区| 97人妻精品一区二区三区麻豆| 国内精品久久久久精免费| av视频在线观看入口| 两个人的视频大全免费| www.精华液| or卡值多少钱| 欧美三级亚洲精品| 精品日产1卡2卡| 国产亚洲欧美在线一区二区| 一级黄色大片毛片| 在线看三级毛片| 深夜精品福利| 不卡av一区二区三区| 1024手机看黄色片| 欧美黑人精品巨大| 在线观看舔阴道视频| 亚洲成av人片免费观看| 国产又色又爽无遮挡免费看| 国产精品野战在线观看| 一边摸一边抽搐一进一小说| 久久国产乱子伦精品免费另类| 国产成人av激情在线播放| 久久久久久久精品吃奶| 亚洲av成人不卡在线观看播放网| 最新美女视频免费是黄的| 18禁美女被吸乳视频| 午夜福利高清视频| 久久伊人香网站| 天天添夜夜摸| 亚洲一区二区三区色噜噜| 亚洲男人的天堂狠狠| 亚洲狠狠婷婷综合久久图片| 亚洲人成网站在线播放欧美日韩| 日本一区二区免费在线视频| 波多野结衣高清无吗| 国产精品永久免费网站| 精品国产乱码久久久久久男人| 亚洲成人久久性| 国产av一区在线观看免费| 免费看美女性在线毛片视频| 成人av一区二区三区在线看| www.熟女人妻精品国产| 制服丝袜大香蕉在线| 99久久久亚洲精品蜜臀av| www.熟女人妻精品国产| 性色av乱码一区二区三区2| 日本免费一区二区三区高清不卡| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美免费精品| 国内少妇人妻偷人精品xxx网站 | 精品国产超薄肉色丝袜足j| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av电影在线进入| 国产一区二区三区视频了| 亚洲国产精品久久男人天堂| 久久中文字幕一级| ponron亚洲| 中文在线观看免费www的网站 | 老汉色∧v一级毛片| 丰满人妻熟妇乱又伦精品不卡| 国产伦人伦偷精品视频|