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

    沼氣制油橇裝裝置的燃爆風(fēng)險分析與優(yōu)化設(shè)計

    2019-11-12 06:18:06鄭曉云陳國明付建民張新琪孔得朋
    天然氣工業(yè) 2019年10期
    關(guān)鍵詞:制油隔板沼氣

    鄭曉云 陳國明 付建民 張 帥 楊 燦 張新琪 孔得朋

    中國石油大學(xué)(華東)海洋油氣裝備與安全技術(shù)研究中心

    0 引言

    沼氣作為一種可再生的生物質(zhì)能源,具有很廣闊的應(yīng)用前景,以沼氣為原料生產(chǎn)合成輕質(zhì)燃油和重質(zhì)燃油,可以有效降低目前我國對于化石燃料的依賴程度,優(yōu)化目前的能源結(jié)構(gòu)[1]。我國對于沼氣制油技術(shù)的研究還處于起步階段,尚未大規(guī)模投入運營。目前來看,國內(nèi)可利用的氣源多為大型養(yǎng)殖場的沼氣池。這類氣源地點分散、氣量不穩(wěn)定,同時還受養(yǎng)殖場經(jīng)營狀況的影響。氣源的以上特點決定了如果采用固定式規(guī)?;朴脱b置不僅建設(shè)成本較高,而且不方便移動。因此,基于節(jié)約建設(shè)成本、方便移動等方面的考慮,國內(nèi)外認(rèn)可的建設(shè)方案均為小型模塊式橇裝化裝置。但是對于該種方案的運行風(fēng)險和安全設(shè)施設(shè)計,尤其是對沼氣制油橇裝裝置泄漏后果的研究目前仍存在著大量空白。沼氣制油裝置的工藝模塊主要有沼氣凈化處理、壓縮與費托合成、儲罐共3個模塊。裝置的主要危險來自于可燃物質(zhì)的泄漏,泄漏的沼氣容易在裝置區(qū)域聚集,形成甲烷高濃度燃爆區(qū)域。

    沼氣制油橇裝化裝置的物料泄漏和火災(zāi)爆炸等風(fēng)險是安全設(shè)計的主要考慮因素,目前相關(guān)針對性研究較少。計算流體動力學(xué)(Computational Fluid Dynamics,以下簡稱CFD)仿真模擬軟件在泄漏以及爆炸后果研究中有著廣泛的應(yīng)用,經(jīng)常被用于計算工業(yè)爆炸產(chǎn)生的超壓。該方法被認(rèn)為比現(xiàn)象學(xué)模型更為準(zhǔn)確[2]。因此,筆者選擇CFD方法對沼氣制油工藝裝置進(jìn)行氣體爆炸風(fēng)險分析。爆炸風(fēng)險的分析多以穩(wěn)態(tài)泄漏速率為基礎(chǔ)對泄漏量進(jìn)行計算,再基于CFD方法建立現(xiàn)場模型進(jìn)行氣體擴(kuò)散和爆炸模擬。但是,基于穩(wěn)態(tài)泄漏模擬得到的爆炸后果過于保守,Gupta和Chan[3]研究發(fā)現(xiàn)全過程動態(tài)泄漏得出的等效可燃?xì)庠企w積Q9值準(zhǔn)確性更高。Нansen等[4]通過等效云的方法,分析了Q5、Q8、Q9等效云方法的優(yōu)缺性并在改變通風(fēng)、點火位置、氣云位置和氣云體積的條件下驗證了Q9等效方法在爆炸測試和理論計算方面的保守性。Shi等[5]研究了BRANN和LMANN兩種方法在估算擁擠海上平臺易燃云尺寸的適用性中,根據(jù)Spouge[6]提出的公式計算瞬態(tài)泄漏速率,但該泄漏速率并未與實際動態(tài)工藝結(jié)合,在爆炸模擬及預(yù)測中造成一定的保守估計。2016年北京市海淀區(qū)紅聯(lián)南村小區(qū)管道破裂導(dǎo)致天然氣爆炸事故,Ji等[7]根據(jù)經(jīng)驗公式計算了泄漏參數(shù),再利用CFD軟件對事故現(xiàn)場進(jìn)行了重建。該研究忽視了由工藝因素影響的泄漏速率與理論計算值的差別。為解決這一問題,將Aspen НYSYS軟件(以下簡稱НYSYS)引入后果分析中,利用其功能模塊模擬化工工藝流程,并指定泄漏點實時模擬工藝過程中發(fā)生的泄漏。作為較成熟的化工流程模擬軟件,НYSYS已廣泛應(yīng)用于化工模擬與優(yōu)化,Sunny等[8]利用НYSYS模擬了再氣化液化天然氣(LNG)的合成氣生產(chǎn)單元,預(yù)測出機(jī)組運行最優(yōu)參數(shù)與合成氣氣體的組成。在動態(tài)泄漏擴(kuò)散數(shù)據(jù)方面,付建民等[9]構(gòu)建了基于НYSYS的泄漏等效模型,就ESD閥關(guān)斷順序和延遲對管道流量、泄漏速率在不同泄漏孔徑下的影響進(jìn)行了研究。但是,上述研究均以相對簡單的天然氣管道或生產(chǎn)裝置為對象,未考慮過程控制系統(tǒng)的干預(yù)對泄漏過程的影響,導(dǎo)致在后果分析上產(chǎn)生較大的偏差。

    筆者以某設(shè)計院設(shè)計的沼氣制油裝置為研究對象,將НYSYS軟件用于沼氣制油過程的泄漏模擬,考慮過程控制對泄漏速率的影響,通過準(zhǔn)確預(yù)測不同工藝安全控制條件下的泄漏量,結(jié)合泄漏點所在區(qū)域的場景得到泄漏影響范圍和火災(zāi)爆炸后果,以期為后續(xù)風(fēng)險分析和本質(zhì)安全研究優(yōu)化提供依據(jù)。

    1 氣體泄漏擴(kuò)散爆炸數(shù)學(xué)模型

    1.1 p-R狀態(tài)方程

    沼氣制油工藝泄漏模擬選擇p-R狀態(tài)方程(Рeng-Robinson)[10]建模,p-R方程式是НYSYS常用的物性方法,可以對單相和多相物流進(jìn)行精確計算,其表達(dá)式為:

    式中p表示壓力,Рa;R表示摩爾氣體常數(shù);T表示開氏溫度,K;V表示摩爾體積,mol/L;a、b均表示特征參數(shù);B表示常數(shù)。

    1.2 氣相擴(kuò)散守恒方程

    沼氣制油裝置泄漏氣體擴(kuò)散和爆炸過程由FLACS軟件計算,F(xiàn)LACS通過有限體積法在笛卡爾網(wǎng)格下求解納維—斯托克斯方程(N-S方程),泄漏過程符合連續(xù)性方程和質(zhì)量、能量、動量守恒方程,流體控制方程為:

    式中φ表示求解變量(質(zhì)量、動量和能量等變量);ρ表示氣體密度,kg/m3;μi表示i方向的速度矢量;xj表示j方向的積分;Γφ表示擴(kuò)散系數(shù);Sφ表示源項。

    1.3 燃燒爆炸方程

    沼氣制油裝置泄漏后產(chǎn)生的可燃?xì)庠圃陟o止的條件下被點燃時,燃燒初期為層流狀態(tài),火焰表面光滑,火焰?zhèn)鞑ビ蔁醾鬟f和分子擴(kuò)散決定。在之后的過渡階段,火焰表面起皺,燃燒速度增加變?yōu)闇?zhǔn)湍流狀態(tài),燃燒后期火焰狀態(tài)變?yōu)橥牧?。層流燃燒速度取決于燃料類型、燃料—空氣混合比和壓力,對于每種燃料,不同當(dāng)量比條件下的層流燃燒速度是確定的,可燃?xì)怏w—空氣混合物的層流燃燒速度用體積加權(quán)平均法計算。壓力條件下層流燃燒速率為:

    式中SL表示層流燃燒速率,kg/s;SL0表示初始層流燃燒速率,kg/s;p表示設(shè)定壓力,Рa;p0表示初始壓力,Рa;γp是燃料決定的經(jīng)驗參數(shù)。

    在準(zhǔn)湍流條件下的燃燒速率:

    式中SQL表示準(zhǔn)湍流燃燒速率,kg/s;χ表示燃料決定的常數(shù);R表示火焰半徑,m。

    湍流燃燒速率(ST)的一般表達(dá)式為:

    其中Karlovitz拉伸因子(K)可表示為:

    式中u'表示湍流強(qiáng)度;Re表示湍流雷諾數(shù)。

    最終可得到:

    式中l(wèi)I表示湍流尺度,m-1;μ表示流體黏度,g/(cm·s)。

    FLACS采用的燃燒速率公式為:

    2 甲烷動態(tài)泄漏擴(kuò)散規(guī)律數(shù)值仿真模型

    2.1 沼氣壓縮工藝模塊泄漏模型

    在沼氣制油裝置的幾個主要工藝模塊中,選擇有代表性的沼氣壓縮工藝模塊做動態(tài)泄漏分析。在НYSYS中,將管線的泄漏位置處設(shè)置閥門,閥門出口處的壓力設(shè)置為大氣壓,打開閥門即可進(jìn)行泄漏模擬。利用НYSYS對沼氣壓縮模塊進(jìn)行過程控制系統(tǒng)干預(yù)下物料動態(tài)泄漏強(qiáng)度模擬,其模型如圖1所示。考慮到沼氣資源工藝比較復(fù)雜,為了方便泄漏強(qiáng)度模擬仿真,在不違背實際工藝原則基礎(chǔ)上對裝置進(jìn)行合理簡化。忽略設(shè)備內(nèi)部和工藝管道之中壓力水頭損失和環(huán)境的影響,并假設(shè)沼氣制油工藝所有裝置與環(huán)境熱交換可以忽略。

    2.2 動態(tài)泄漏強(qiáng)度計算

    設(shè)備發(fā)生泄漏后的計算時間域取300 s,動態(tài)泄漏過程分為2個階段:①第一階段,前60 s在前端ESD未關(guān)斷進(jìn)料時,僅在過程控制系統(tǒng)影響下泄漏;②第二階段,后240 s ESD系統(tǒng)檢測到泄漏自動關(guān)斷進(jìn)料,物料在重力作用下泄漏。根據(jù)筆者所取的泄漏孔徑標(biāo)準(zhǔn),沼氣壓縮模塊發(fā)生3種孔徑泄漏的泄漏速率動態(tài)變化曲線如圖2所示,不同孔徑泄漏前60 s和后240 s泄漏速率下降如圖3所示。

    圖1 沼氣壓縮模塊НYSYS模型圖

    圖2 不同泄漏孔徑下沼氣壓縮模塊泄漏速率圖

    圖3 泄漏速率下降率圖

    對于小孔泄漏(5 mm),在前端泄漏開始一分鐘內(nèi)未關(guān)閉進(jìn)料,泄漏速率在60 s時比0秒時下降0.44%,泄漏60~240 s內(nèi)泄漏速率下降3.92%。因此5 mm小孔泄漏對裝置內(nèi)物流壓力影響較小,未進(jìn)行前端ESD進(jìn)料關(guān)斷時裝置仍可保壓運行。對于20 mm和50 mm泄漏孔徑,前端進(jìn)料關(guān)斷前的60 s內(nèi)泄漏速率分別下降5.55%和21.06%,不能保壓運行,當(dāng)進(jìn)行前端關(guān)斷后泄漏速率急劇下降,在關(guān)斷后240 s泄漏時間內(nèi),泄漏速率分別下降46.03%和84.29%。因此對于中孔和大孔泄漏,后300 s的泄漏時間里裝置內(nèi)物料會大量外泄。

    3 橇裝模塊安全間距優(yōu)化

    3.1 沼氣制油裝置物理模型與網(wǎng)格劃分

    參照工藝資料,在FLACS中建立3個小型化橇裝模型(圖4),分別代表沼氣制油裝置的3個工藝模塊,即沼氣凈化處理模塊(左)、壓縮與費托合成模塊(中)、儲罐模塊(右),每個橇裝模塊(以下簡稱橇塊)尺寸為12 m×5 m×5 m。對泄漏擴(kuò)散場景,采用0.4 m均勻網(wǎng)格,與泄漏點射流方向垂直的網(wǎng)格進(jìn)行0.2 m局部加密,為了能讓FLACS正確識別并計算,通過Smooth功能將疏密網(wǎng)格之間有限體積元的厚度變化控制在40%以內(nèi)[11]。

    圖4 大當(dāng)量可燃?xì)庠颇P蛨D

    3.2 安全間距和點火位置對蒸氣云爆炸超壓的影響

    安全間距作為一個開放空間,是一項廣泛應(yīng)用的安全設(shè)計,其原理可以有效防止湍流的產(chǎn)生,切斷增強(qiáng)反應(yīng)物混合、火焰加速和超壓氣體爆炸的正反饋機(jī)制。為了對比不同的安全間距對爆炸超壓的減緩作用,相鄰橇塊之間安全間距分別按照3 m、6 m、9 m進(jìn)行計算(圖4)。圖4中紅色區(qū)域為可燃?xì)庠莆恢?,在不同安全間距條件下研究橇塊內(nèi)充滿大量可燃?xì)怏w(充滿3個橇塊)和橇塊充滿少量可燃?xì)怏w(充滿1個橇塊)兩種情景,可燃?xì)怏w充滿3個橇塊的情景定義為大當(dāng)量可燃?xì)庠颇P停?種不同安全間距場景采用相同體積的可燃?xì)庠?,泄漏點都位于圖4中右側(cè)儲罐橇塊,以便于進(jìn)行對比。少量場景假設(shè)可燃?xì)庠浦怀錆M左邊的沼氣凈化處理橇塊。

    在FLACS中可以通過定義非活性泄壓板來監(jiān)測給定區(qū)域面積上的爆炸壓力載荷,大體積的防爆墻、甲板、設(shè)備外殼等都可直接通過泄壓板獲得爆炸壓力載荷監(jiān)測值,而對于小型設(shè)備,則需要通過定義多個泄壓板,對監(jiān)測數(shù)據(jù)進(jìn)行積分后獲得。文中通過定義非活性泄壓板獲得爆炸過程中沖擊波在平面內(nèi)的平均壓力載荷,可作為裝置承壓載荷設(shè)計的標(biāo)準(zhǔn)。監(jiān)測壓力的非活性泄壓板設(shè)置在每2個橇塊安全間距的中心位置,以監(jiān)測爆炸沖擊波對相鄰橇塊的影響及其在安全間距中的衰減過程。

    合理的安全間距可降低爆炸后火焰?zhèn)鞑ピ斐傻耐牧鲾_動,但FLACS在計算大安全間距時對爆炸沖擊波傳播計算失真,需要通過Data-Dump(數(shù)據(jù)轉(zhuǎn)儲)方法進(jìn)行修正,對于大型裝置的小安全間距(安全間距與裝置尺寸之比介于0.125~1),F(xiàn)LACS的計算值則與實驗值吻合良好。本研究中的橇塊安全間距為3 m、6 m、9 m,單橇塊裝置最大尺寸為 12 m,屬于小安全間距范圍,可直接通過FLACS進(jìn)行爆炸仿真模擬,其結(jié)果具有較高的可信度。點火位置對可燃?xì)庠频谋ê蠊簿哂泻艽笥绊?,考慮到這一因素,筆者研究了點火源位于橇塊內(nèi)和橇塊間這2種點火方式,得到點火源Y/Z截面上爆炸時間內(nèi)的最大爆炸超壓(pmax)如圖5、6所示。

    在上述2種點火方式中,對比3種不同安全間距和點火位置,在可燃?xì)庠飘?dāng)量體積(Q9)不變的前提下,可以看出:①隨著安全間距的增大,爆炸超壓區(qū)域減小,核心裝置區(qū)域內(nèi)的爆炸超壓明顯降低;②點火位置位于橇塊內(nèi)的pmax和影響面積均小于橇塊間,但差別較小。上述研究結(jié)果表明:對于小型橇裝化沼氣制油裝置,增大橇塊間的安全間距可以明顯降低爆炸火焰?zhèn)鞑ピ斐傻耐牧鲾_動,可以降低爆炸超壓造成的破壞,縮小影響區(qū)域。全區(qū)域可燃?xì)庠票▓鼍胺腔钚孕箟喊灞O(jiān)測到的爆炸超壓平均壓強(qiáng)如表1所示。

    圖5 點火源位于不同安全間距橇塊內(nèi)部時最大爆炸超壓圖

    圖6 點火源位于不同安全間距橇塊間最大爆炸超壓圖

    點火源位于橇塊內(nèi)部時,在安全間距相同的情況下泄壓監(jiān)測板的平均超壓值隨著與點火源距離的增加,超壓值有先增加后降低的趨勢,而橇塊間點火方式則未出現(xiàn)明顯下降。這是由于隨著火焰?zhèn)鞑ゾ嚯x的增加,湍流擾動越來越強(qiáng)烈,加劇了火焰燃燒速度[12-13],而隨著與點火源距離的增加,可燃?xì)庠仆獠康目扇細(xì)怏w被爆炸沖散,使氣體濃度降低,故泄壓監(jiān)測板監(jiān)測到的平均超壓會出現(xiàn)下降趨勢。橇塊與橇塊間的擁塞度極低,若點火源設(shè)置在橇塊之間,燃燒初級階段不足以造成大的湍流擾動[14-15]。因此該點火方式僅在小型化裝置中能夠產(chǎn)生湍流擾動,在筆者研究的監(jiān)測區(qū)域內(nèi)未出現(xiàn)泄壓監(jiān)測板與點火源距離增加爆炸超壓數(shù)值明顯下降的情況。對分布在單一橇塊內(nèi)的少量可燃?xì)庠?,研究結(jié)果表明,安全間距大小基本不影響爆炸超壓的擴(kuò)散范圍。

    表1 全區(qū)域可燃?xì)庠票▓鼍胺腔钚孕箟喊灞O(jiān)測的平均超壓值表 Рa

    4 保護(hù)隔板排布方式優(yōu)化

    4.1 隔板排布方式

    對于橇裝化沼氣制油裝置,在四周添加隔板可降低設(shè)備外部環(huán)境造成的腐蝕,在生產(chǎn)和吊裝時起到保護(hù)設(shè)備的作用,但增加保護(hù)隔板會增大橇塊的擁塞度,形成半密閉空間,影響泄漏氣體擴(kuò)散,對爆炸后果產(chǎn)生影響[16-17]。筆者對保護(hù)隔板進(jìn)行優(yōu)化設(shè)計,通過改變橇塊兩側(cè)及兩端隔板的數(shù)量,計算以下6組不同保護(hù)隔板類型:3個橇塊均無保護(hù)隔板(圖7-a)、3個橇塊兩側(cè)都有保護(hù)隔板(圖7-b)、壓縮與費托合成模塊的兩側(cè)、沼氣凈化處理模塊的左側(cè)和儲罐模塊的右側(cè)設(shè)置隔板(圖7-c)、只在沼氣凈化模塊和儲罐模塊的一側(cè)設(shè)置隔板(圖7-d、e)、在每個橇塊的兩端布置隔板(圖7-f)。采用3 m橇塊間距,每組保護(hù)隔板進(jìn)行6組不同泄漏方向研究(±X、±Y、±Z),共得到108組不同泄漏場景,通過FLACS泄漏擴(kuò)散模塊(Release and Dispersion)對其進(jìn)行擴(kuò)散仿真模擬[18],以期得到最優(yōu)的保護(hù)隔板排布方式,將動態(tài)泄漏模擬得到的數(shù)據(jù)作為FLACS泄漏模塊的輸入數(shù)據(jù)。

    4.2 不同保護(hù)隔板排布下爆炸后果

    圖7 6組不同保護(hù)隔板類型示意圖

    對于保護(hù)隔板設(shè)計對爆炸后果的影響,將大當(dāng)量可燃?xì)庠坪驼託庵朴凸に嚄l件下動態(tài)泄漏產(chǎn)生可燃?xì)庠?種情況進(jìn)行對比[20]。在實際工藝條件下凈化處理模塊和壓縮與費托合成模塊發(fā)生5 mm孔徑泄漏都不能產(chǎn)生有效可燃蒸氣云,不具有爆炸風(fēng)險。凈化處理模塊中20 mm孔徑和50 mm孔徑凈化沼氣裝置泄漏在不同保護(hù)隔板類型、不同泄漏方向下的最大等效可燃?xì)庠企w積(Q9)FLACS計算值如圖8所示。沼氣壓縮與費托合成模塊20 mm孔徑泄漏(圖7-b場景兩側(cè)加隔板) 產(chǎn)生的最大等效可燃?xì)庠企w積為5.7 m3,其余保護(hù)隔板排列方式下Q9<1 m3,其爆炸風(fēng)險很小。沼氣壓縮模塊50 mm孔徑泄漏最大爆炸超壓pmax與可燃?xì)庠瞥錆M橇塊的大當(dāng)量氣云爆炸超壓對比如表2所示。

    圖7中(a~f)所示的保護(hù)隔板排布類型,在氣云爆炸中心高度水平最大超壓值對應(yīng)的截面如圖9所示,在大當(dāng)量氣云爆炸條件下,無保護(hù)隔板(圖9-a)時產(chǎn)生的爆炸超壓最小,雙側(cè)增加保護(hù)隔板(圖9-b)會大幅度提高裝置的封閉性,造成最高17.55 kРa的爆炸超壓,遠(yuǎn)大于其他場景(圖9-c~f)。在實際工況泄漏條件下,根據(jù)FLACS動態(tài)泄漏擴(kuò)散結(jié)果,將等效氣云體積值導(dǎo)入爆炸模塊,得到的爆炸超壓為5.784 kРa,也遠(yuǎn)大于其余保護(hù)隔板排布方式。

    由此可知,增加保護(hù)隔板可以明顯控制沖擊波在該方向的傳播。雙側(cè)加保護(hù)隔板爆炸的pmax較之不加隔板保護(hù)措施顯著提升,并且隔板可以顯著降低火焰在該方向的傳播。但是,單側(cè)增加保護(hù)隔板可以在不明顯增加爆炸后果基礎(chǔ)上起到相應(yīng)的保護(hù)作用。

    5 結(jié)論與建議

    1)沼氣制油實際工藝條件下沼氣壓縮模塊小孔孔徑(5 mm)泄漏時前60 s的泄漏速率基本不變,可保壓運行,前端進(jìn)料閥門關(guān)斷后240 s泄漏時間內(nèi)泄漏速率降低3.92%,流量壓力儀表無法發(fā)現(xiàn)泄漏特征。當(dāng)泄漏孔徑為20 mm和50 mm時該值分別為46.03%和84.29%,說明中孔和大孔泄漏會導(dǎo)致裝置內(nèi)物料大量外泄,流量壓力等過程監(jiān)控儀表會有顯著響應(yīng)。設(shè)計中應(yīng)對易發(fā)小孔泄漏部位重點輔以泄漏監(jiān)控。

    圖8 不同隔板排布類型下兩種孔徑泄漏時的Q9值圖

    表2 不同保護(hù)隔板排布類型下最大爆炸超壓圖 kРa

    2)對于沼氣制油橇裝化裝置,適當(dāng)增大橇塊安全間距會大大降低氣體爆炸超壓強(qiáng)度,其中爆炸遠(yuǎn)端壓力監(jiān)測板的爆炸超壓最大值在安全間距為3 m時為 721.1 Рa,而安全間距為 9 m 時降低為 83.68 Рa ;其他變量一致時,點火源位于橇塊內(nèi)部產(chǎn)生的爆炸超壓最大值明顯大于兩橇塊間。

    3)橇塊兩側(cè)加保護(hù)隔板(圖9-b場景),在裝置實際工況下50 mm大孔動態(tài)泄漏導(dǎo)致的爆炸超壓最大值為5.784 kРa,當(dāng)研究可燃?xì)怏w充滿整個橇塊理想狀況時該值為17.55 kРa,兩種爆炸場景下產(chǎn)生的超壓值都遠(yuǎn)高于其余排布方式。

    圖9 不同保護(hù)隔板排布類型下爆炸超壓最大值二維截圖

    4)建議在橇塊單側(cè)增加保護(hù)隔板,在起到泄壓作用同時一旦發(fā)生爆炸不會產(chǎn)生過高爆炸超壓,有利于降低爆炸后果;橇塊間安全間距建議設(shè)置為6 m或更大,在降低爆炸后果前提下盡可能節(jié)省裝置的空間占位。

    猜你喜歡
    制油隔板沼氣
    基于有限元分析的隔板工藝支撐技術(shù)研究
    壓力容器隔板的一種設(shè)計方法
    第四章 化糞土為力量——沼氣能
    第四章 化糞土為力量——沼氣能
    橫隔板參數(shù)對裝配式箱梁寬橋受力的影響
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    《中國沼氣》征稿簡則
    中國沼氣(2017年4期)2017-11-08 01:55:08
    水泡糞工藝制沼氣及沼氣的凈化提純
    上海煤氣(2016年1期)2016-05-09 07:12:35
    國內(nèi)首套百萬噸級煤間接制油項目試車成功
    國內(nèi)首套百萬噸級煤間接制油項目試車成功
    80岁老熟妇乱子伦牲交| 国产男人的电影天堂91| 在线观看免费日韩欧美大片 | 亚洲色图综合在线观看| 校园人妻丝袜中文字幕| 美女福利国产在线| 高清毛片免费看| 亚洲久久久国产精品| 久久99一区二区三区| 日本-黄色视频高清免费观看| 国产精品久久久久久久电影| 交换朋友夫妻互换小说| 国产欧美另类精品又又久久亚洲欧美| 2021少妇久久久久久久久久久| 免费黄网站久久成人精品| 中文字幕免费在线视频6| 老熟女久久久| 欧美日韩综合久久久久久| 亚洲天堂av无毛| 久久久精品94久久精品| 亚洲色图综合在线观看| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看 | 亚洲一级一片aⅴ在线观看| 色视频在线一区二区三区| 午夜激情福利司机影院| 久久久国产欧美日韩av| 精品人妻在线不人妻| 女性生殖器流出的白浆| 91aial.com中文字幕在线观看| 五月伊人婷婷丁香| 精品久久蜜臀av无| 日日摸夜夜添夜夜爱| 免费大片18禁| 免费看av在线观看网站| 久久久精品免费免费高清| 91aial.com中文字幕在线观看| 久久99一区二区三区| 少妇人妻 视频| 人人澡人人妻人| 亚洲国产色片| 午夜影院在线不卡| 精品熟女少妇av免费看| 久久久久久久久大av| 国产综合精华液| 日韩 亚洲 欧美在线| 极品少妇高潮喷水抽搐| 人人妻人人澡人人爽人人夜夜| 国产一区二区三区av在线| 狂野欧美激情性xxxx在线观看| 亚洲欧美中文字幕日韩二区| 一区二区日韩欧美中文字幕 | 亚洲av欧美aⅴ国产| 日日撸夜夜添| 夜夜骑夜夜射夜夜干| 国产爽快片一区二区三区| 欧美激情 高清一区二区三区| 天天操日日干夜夜撸| 99久久中文字幕三级久久日本| 日日爽夜夜爽网站| 高清毛片免费看| av有码第一页| 精品久久久久久久久亚洲| 国产片内射在线| 曰老女人黄片| 欧美变态另类bdsm刘玥| 中文欧美无线码| 欧美少妇被猛烈插入视频| 欧美日韩视频高清一区二区三区二| 久久影院123| 精品一区在线观看国产| 午夜久久久在线观看| 日本vs欧美在线观看视频| 最黄视频免费看| 亚洲精品日韩av片在线观看| 久久久午夜欧美精品| 91久久精品电影网| 日产精品乱码卡一卡2卡三| av播播在线观看一区| 少妇人妻 视频| 国产精品免费大片| 简卡轻食公司| 久久久精品免费免费高清| 男女啪啪激烈高潮av片| 丰满迷人的少妇在线观看| 午夜激情福利司机影院| 18+在线观看网站| 一级,二级,三级黄色视频| 97在线视频观看| 最近2019中文字幕mv第一页| 欧美激情 高清一区二区三区| 免费人妻精品一区二区三区视频| 激情五月婷婷亚洲| 麻豆精品久久久久久蜜桃| 日韩欧美一区视频在线观看| 2022亚洲国产成人精品| 国产高清国产精品国产三级| 嘟嘟电影网在线观看| 在线精品无人区一区二区三| 亚洲天堂av无毛| av国产精品久久久久影院| 久久精品久久久久久久性| 国产成人一区二区在线| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| av在线app专区| 一级毛片黄色毛片免费观看视频| 亚洲一级一片aⅴ在线观看| 久久 成人 亚洲| 亚洲不卡免费看| 成人免费观看视频高清| 免费不卡的大黄色大毛片视频在线观看| 最后的刺客免费高清国语| 乱人伦中国视频| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 欧美少妇被猛烈插入视频| 有码 亚洲区| 国产色爽女视频免费观看| 久久精品国产亚洲网站| 久久久国产欧美日韩av| 欧美人与善性xxx| 午夜福利网站1000一区二区三区| 国产熟女午夜一区二区三区 | av不卡在线播放| 韩国高清视频一区二区三区| av线在线观看网站| 亚洲精品一区蜜桃| 亚洲精品国产av蜜桃| 香蕉精品网在线| 免费看av在线观看网站| 婷婷成人精品国产| 在线观看人妻少妇| 美女中出高潮动态图| 蜜桃在线观看..| 丝袜喷水一区| 在线精品无人区一区二区三| 中文字幕人妻丝袜制服| 日韩伦理黄色片| 亚洲精品av麻豆狂野| 久久久久久久大尺度免费视频| 日韩,欧美,国产一区二区三区| 国产精品无大码| 国产免费视频播放在线视频| 夜夜爽夜夜爽视频| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产色婷婷电影| 狂野欧美白嫩少妇大欣赏| 国产又色又爽无遮挡免| 欧美人与性动交α欧美精品济南到 | 夫妻午夜视频| 丰满饥渴人妻一区二区三| 国产成人午夜福利电影在线观看| 2022亚洲国产成人精品| 黄色怎么调成土黄色| 超碰97精品在线观看| 91精品三级在线观看| 欧美人与性动交α欧美精品济南到 | 97超碰精品成人国产| 国产黄色免费在线视频| 人人妻人人添人人爽欧美一区卜| 中文字幕制服av| 久久女婷五月综合色啪小说| 99热网站在线观看| 狂野欧美激情性bbbbbb| 国产欧美亚洲国产| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区国产| 人妻一区二区av| 美女大奶头黄色视频| 亚洲欧洲国产日韩| 中文精品一卡2卡3卡4更新| 国产男女内射视频| 免费av中文字幕在线| 国国产精品蜜臀av免费| 久久国产亚洲av麻豆专区| 性色av一级| 免费黄频网站在线观看国产| 99re6热这里在线精品视频| 亚洲精品aⅴ在线观看| 天天躁夜夜躁狠狠久久av| 国产精品一国产av| 伊人亚洲综合成人网| 我的老师免费观看完整版| 国产精品成人在线| 又粗又硬又长又爽又黄的视频| 美女福利国产在线| 国产 精品1| 久久精品国产a三级三级三级| 日韩在线高清观看一区二区三区| 欧美最新免费一区二区三区| 99国产综合亚洲精品| 成人手机av| 日韩三级伦理在线观看| 一区二区三区四区激情视频| 国产老妇伦熟女老妇高清| 日日啪夜夜爽| 久久久午夜欧美精品| 中文字幕av电影在线播放| tube8黄色片| 男女啪啪激烈高潮av片| 一区二区三区免费毛片| 亚洲av成人精品一区久久| 国产极品天堂在线| 国产高清有码在线观看视频| 久久国产精品男人的天堂亚洲 | 日本猛色少妇xxxxx猛交久久| 久久久久国产精品人妻一区二区| 波野结衣二区三区在线| 色视频在线一区二区三区| 久久久久久久久久久丰满| 久久99一区二区三区| 人人澡人人妻人| 亚洲精品一区蜜桃| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 中文字幕久久专区| av福利片在线| 人妻少妇偷人精品九色| 国产成人精品在线电影| 男女高潮啪啪啪动态图| 国产亚洲av片在线观看秒播厂| 亚洲精品aⅴ在线观看| 黄片播放在线免费| 一级爰片在线观看| 超色免费av| 国产免费视频播放在线视频| 在线播放无遮挡| 亚洲伊人久久精品综合| 丝袜脚勾引网站| 伦精品一区二区三区| 国产一区二区在线观看av| 久久久久精品性色| 国产男女超爽视频在线观看| 在线观看国产h片| 亚洲av在线观看美女高潮| 美女大奶头黄色视频| 啦啦啦中文免费视频观看日本| av黄色大香蕉| 男女边摸边吃奶| 免费观看性生交大片5| 乱码一卡2卡4卡精品| 日产精品乱码卡一卡2卡三| 精品国产国语对白av| 精品一品国产午夜福利视频| 一区二区av电影网| 免费看av在线观看网站| 欧美精品人与动牲交sv欧美| 亚洲国产最新在线播放| 国产精品人妻久久久久久| 久久久久精品久久久久真实原创| 人妻系列 视频| 久久久久久久久久人人人人人人| 亚洲欧洲日产国产| 中文字幕亚洲精品专区| 男女边摸边吃奶| av卡一久久| 一区二区日韩欧美中文字幕 | 国产欧美另类精品又又久久亚洲欧美| 欧美 日韩 精品 国产| 亚洲国产精品一区二区三区在线| 亚洲欧美日韩卡通动漫| 97超碰精品成人国产| 国产精品一国产av| 欧美日韩av久久| 成人免费观看视频高清| 久久人人爽av亚洲精品天堂| 久久狼人影院| 国产成人免费观看mmmm| 麻豆精品久久久久久蜜桃| 美女脱内裤让男人舔精品视频| 中国国产av一级| 在线观看人妻少妇| 久久久久精品久久久久真实原创| 最近2019中文字幕mv第一页| 久久久久久久久久人人人人人人| 最近中文字幕高清免费大全6| 国产一区有黄有色的免费视频| 欧美精品一区二区免费开放| 亚洲精品国产色婷婷电影| 王馨瑶露胸无遮挡在线观看| 欧美人与性动交α欧美精品济南到 | 99热6这里只有精品| 老熟女久久久| 99热6这里只有精品| 精品人妻在线不人妻| 高清午夜精品一区二区三区| 精品久久久久久电影网| 国产精品免费大片| 大片免费播放器 马上看| 99久久人妻综合| 国产白丝娇喘喷水9色精品| 欧美精品一区二区免费开放| 青春草亚洲视频在线观看| 99热全是精品| 亚洲精品一二三| 亚洲人成77777在线视频| 精品久久久久久久久亚洲| 黄色怎么调成土黄色| 香蕉精品网在线| 女性生殖器流出的白浆| 99热网站在线观看| videossex国产| 涩涩av久久男人的天堂| 最近中文字幕2019免费版| 久久久久久伊人网av| av电影中文网址| 欧美97在线视频| 2022亚洲国产成人精品| 欧美 亚洲 国产 日韩一| www.av在线官网国产| 中文字幕最新亚洲高清| 国产精品人妻久久久久久| 国产精品人妻久久久影院| 在线亚洲精品国产二区图片欧美 | 国产精品久久久久久久久免| 国产精品成人在线| 久久精品国产亚洲网站| 天天躁夜夜躁狠狠久久av| 中文字幕精品免费在线观看视频 | 99热网站在线观看| 久久精品国产a三级三级三级| 欧美日韩在线观看h| 曰老女人黄片| 国产成人精品无人区| 午夜福利,免费看| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 亚洲精品一区蜜桃| 精品国产一区二区久久| 如何舔出高潮| 69精品国产乱码久久久| 亚洲精品456在线播放app| 亚洲av欧美aⅴ国产| 一本大道久久a久久精品| videosex国产| 美女福利国产在线| 久久精品国产亚洲av涩爱| a级毛片免费高清观看在线播放| 国产精品麻豆人妻色哟哟久久| 最黄视频免费看| 午夜激情福利司机影院| 永久免费av网站大全| 日韩免费高清中文字幕av| 精品人妻熟女毛片av久久网站| 免费久久久久久久精品成人欧美视频 | 久久午夜综合久久蜜桃| 日韩精品免费视频一区二区三区 | 午夜久久久在线观看| 老女人水多毛片| 少妇被粗大猛烈的视频| 如何舔出高潮| 日本色播在线视频| 久久99热6这里只有精品| 国产免费一级a男人的天堂| av国产久精品久网站免费入址| 久久这里有精品视频免费| 亚洲一区二区三区欧美精品| 欧美bdsm另类| 蜜桃在线观看..| 国产精品麻豆人妻色哟哟久久| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| 国产成人aa在线观看| 一个人免费看片子| 高清黄色对白视频在线免费看| 18禁裸乳无遮挡动漫免费视频| 一本一本综合久久| 99久久精品一区二区三区| 这个男人来自地球电影免费观看 | 丰满迷人的少妇在线观看| 纯流量卡能插随身wifi吗| 欧美xxⅹ黑人| 丰满乱子伦码专区| 国产精品久久久久久精品古装| 五月开心婷婷网| 五月伊人婷婷丁香| 99九九在线精品视频| 三上悠亚av全集在线观看| 国产午夜精品一二区理论片| 精品酒店卫生间| 亚洲av男天堂| 在线观看免费高清a一片| 成人国语在线视频| 国产精品国产三级国产av玫瑰| 亚洲国产精品成人久久小说| 色5月婷婷丁香| 欧美日韩成人在线一区二区| 久久久久精品久久久久真实原创| 高清欧美精品videossex| 国产av一区二区精品久久| 这个男人来自地球电影免费观看 | 午夜免费鲁丝| 日日啪夜夜爽| 日韩av在线免费看完整版不卡| 国产高清有码在线观看视频| 午夜福利视频在线观看免费| av网站免费在线观看视频| 观看美女的网站| 国产精品麻豆人妻色哟哟久久| 成人手机av| 国产视频内射| 在线观看免费日韩欧美大片 | 精品久久久久久久久亚洲| 国产精品久久久久久久久免| av福利片在线| 老熟女久久久| 青春草视频在线免费观看| 亚州av有码| 国产高清三级在线| 99热网站在线观看| 又黄又爽又刺激的免费视频.| 精品一区二区免费观看| 丁香六月天网| 亚洲一级一片aⅴ在线观看| videos熟女内射| 少妇人妻久久综合中文| 成年av动漫网址| 成人国产麻豆网| 国产欧美日韩一区二区三区在线 | 中文字幕人妻熟人妻熟丝袜美| 少妇人妻 视频| 高清毛片免费看| 精品国产一区二区久久| 国产亚洲精品第一综合不卡 | 久久国内精品自在自线图片| 两个人的视频大全免费| 狂野欧美激情性bbbbbb| av在线观看视频网站免费| 国产视频内射| 观看av在线不卡| 18禁观看日本| 亚洲怡红院男人天堂| 女的被弄到高潮叫床怎么办| 少妇被粗大猛烈的视频| 乱人伦中国视频| 99热网站在线观看| 亚洲第一av免费看| 男女边吃奶边做爰视频| a级片在线免费高清观看视频| 久久午夜综合久久蜜桃| 国产成人一区二区在线| 国产免费视频播放在线视频| 最近中文字幕2019免费版| 成人无遮挡网站| 最近中文字幕2019免费版| 亚洲欧美清纯卡通| 久久久久久久久久久久大奶| 麻豆乱淫一区二区| 日本av免费视频播放| 如何舔出高潮| 日日撸夜夜添| 一级毛片aaaaaa免费看小| 91精品国产国语对白视频| 人人妻人人爽人人添夜夜欢视频| 黄片播放在线免费| 国产免费现黄频在线看| 99久国产av精品国产电影| 纵有疾风起免费观看全集完整版| 亚洲精品第二区| 国产极品天堂在线| 飞空精品影院首页| 国产成人91sexporn| 国产高清三级在线| 久久久国产一区二区| 国产精品国产av在线观看| 国产精品女同一区二区软件| 大香蕉久久网| 国产精品久久久久成人av| 午夜福利网站1000一区二区三区| 秋霞在线观看毛片| 欧美日韩精品成人综合77777| 蜜桃久久精品国产亚洲av| 七月丁香在线播放| 日韩精品免费视频一区二区三区 | av卡一久久| 一级毛片黄色毛片免费观看视频| 99久国产av精品国产电影| 高清午夜精品一区二区三区| 欧美精品亚洲一区二区| 欧美日韩视频精品一区| 高清毛片免费看| 亚洲一级一片aⅴ在线观看| 成人综合一区亚洲| 夫妻午夜视频| 在线精品无人区一区二区三| 大片免费播放器 马上看| 免费观看av网站的网址| 亚洲欧美日韩另类电影网站| 日本午夜av视频| 97在线人人人人妻| 精品久久久久久久久av| 日韩成人伦理影院| 人人妻人人添人人爽欧美一区卜| 日韩熟女老妇一区二区性免费视频| av电影中文网址| av网站免费在线观看视频| 麻豆精品久久久久久蜜桃| 日韩中字成人| 日本91视频免费播放| 精品午夜福利在线看| 国产不卡av网站在线观看| 国精品久久久久久国模美| 成人亚洲欧美一区二区av| 老司机亚洲免费影院| 汤姆久久久久久久影院中文字幕| 久久精品国产a三级三级三级| 精品久久久久久久久av| 欧美激情极品国产一区二区三区 | 在线观看一区二区三区激情| 国产精品女同一区二区软件| 在线天堂最新版资源| 51国产日韩欧美| 最新中文字幕久久久久| 国产色婷婷99| 91在线精品国自产拍蜜月| 亚洲第一区二区三区不卡| 丰满饥渴人妻一区二区三| 久久久久久久久大av| 国产精品国产av在线观看| 国产精品久久久久久av不卡| a级毛片在线看网站| 18禁观看日本| 多毛熟女@视频| av有码第一页| 欧美精品亚洲一区二区| 亚洲美女视频黄频| 国产午夜精品一二区理论片| 大片免费播放器 马上看| 日韩一本色道免费dvd| 欧美日韩亚洲高清精品| 2018国产大陆天天弄谢| 免费看光身美女| 日韩亚洲欧美综合| 午夜av观看不卡| 在线观看一区二区三区激情| 黄色怎么调成土黄色| 亚洲第一区二区三区不卡| 国产无遮挡羞羞视频在线观看| 日日撸夜夜添| 久久av网站| 美女cb高潮喷水在线观看| 中文欧美无线码| 久久久久久久亚洲中文字幕| 人人澡人人妻人| 日本-黄色视频高清免费观看| 免费看av在线观看网站| 亚洲国产欧美在线一区| 看免费成人av毛片| 国产片特级美女逼逼视频| 3wmmmm亚洲av在线观看| 久久这里有精品视频免费| 美女福利国产在线| 国产黄色免费在线视频| 涩涩av久久男人的天堂| 亚洲情色 制服丝袜| 99国产综合亚洲精品| 丰满饥渴人妻一区二区三| 一本色道久久久久久精品综合| 欧美激情 高清一区二区三区| 亚洲欧美清纯卡通| 日本色播在线视频| 在线观看免费视频网站a站| 91精品伊人久久大香线蕉| 亚洲精品自拍成人| 亚洲精品一区蜜桃| 一级二级三级毛片免费看| 欧美三级亚洲精品| 晚上一个人看的免费电影| 亚洲国产av新网站| 免费看av在线观看网站| 色婷婷久久久亚洲欧美| 亚洲精品日韩在线中文字幕| 亚洲人成网站在线观看播放| 一本色道久久久久久精品综合| 在线观看免费视频网站a站| 久久久国产一区二区| videosex国产| 一本大道久久a久久精品| 日本色播在线视频| 免费看av在线观看网站| 最新中文字幕久久久久| 免费看不卡的av| 黑人欧美特级aaaaaa片| 成人毛片60女人毛片免费| 高清毛片免费看| 久久久久视频综合| 免费看光身美女| 好男人视频免费观看在线| av播播在线观看一区| 午夜福利视频在线观看免费| 又粗又硬又长又爽又黄的视频| 国产精品.久久久| 日韩欧美精品免费久久| 男人操女人黄网站| 亚洲国产成人一精品久久久| 亚洲欧洲精品一区二区精品久久久 | 黑人猛操日本美女一级片| 国产精品一区二区三区四区免费观看| 午夜老司机福利剧场| 亚洲天堂av无毛| 国产爽快片一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 国产亚洲一区二区精品| 日韩人妻高清精品专区| 亚洲图色成人| 日韩制服骚丝袜av| 亚洲久久久国产精品| 黄色视频在线播放观看不卡| 街头女战士在线观看网站| 一级毛片aaaaaa免费看小| 国产精品一区二区在线不卡| 亚洲av日韩在线播放| 亚洲精品一区蜜桃| 中文精品一卡2卡3卡4更新|