含水率對(duì)油水乳化液傳輸頻率特性影響分析

陳 彬， 劉 閣

(重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心 重慶,400067)

10.16450/j.cnki.issn.1004-6801.2017.05.002

含水率對(duì)油水乳化液傳輸頻率特性影響分析

陳 彬， 劉 閣

(重慶工商大學(xué)廢油資源化技術(shù)與裝備教育部工程研究中心 重慶,400067)

為了掌握油水乳化液中含水率對(duì)其動(dòng)態(tài)特性的影響，建立含水率與乳化液傳輸?shù)念l率特性之間的非線性關(guān)系。首先,推導(dǎo)了油水乳化液在管體中壓力波傳播速度公式，分析了含水率對(duì)乳化液傳輸速度的影響規(guī)律;然后,利用分布參數(shù)理論以及傳遞矩陣法建立了油水乳化液的傳輸頻率模型，表明在PVC透明鋼絲軟管中考慮流固耦合作用時(shí)含水率增加，其傳播速度呈降低趨勢(shì)；反之，鋼管內(nèi)乳化液的傳輸速度隨著含水率的增加而呈遞增趨勢(shì)；含水率對(duì)系統(tǒng)的諧振頻率影響也較為明顯，乳化液中含水率越大，系統(tǒng)的同階諧振頻率越低，且隨著含水率的增加，諧振頻率的階次越大，同階次的諧振頻率之間的差值也相應(yīng)變大。實(shí)驗(yàn)驗(yàn)證所建模型是有效的，為有效檢測(cè)油中含水率以及保障系統(tǒng)安全運(yùn)行提供一定的理論基礎(chǔ)。

含水率； 油水乳化液； 傳播速度； 頻率特性

引 言

油液在運(yùn)行過程中會(huì)因?yàn)榉N種原因混入水分，使油液理化性能發(fā)生變化[1]，極易生成穩(wěn)定的乳化液、加速油的氧化過程和降低油潤(rùn)滑性，仍然會(huì)影響整個(gè)系統(tǒng)的正常工作，且油水乳化液的非恒定流的動(dòng)力學(xué)特性與純油液有很多不同，是僅次于固體顆粒物的一種污染。采取合理方法分析含水率不同的油水乳化液的非恒定流的動(dòng)力學(xué)特性，及時(shí)采取合理的治理措施一直是工業(yè)用油研究熱點(diǎn)之一。

目前，對(duì)油水乳化液的研究主要是液滴粒徑及其分布、相界面形狀及其變化、相含率和相界面濃度等方面。Kang 等[2]研究了方形管道內(nèi)氣液兩相流流型和多孔介質(zhì)內(nèi)的侵滲現(xiàn)象。Suga 等[3]研究了壁面粗糙度和滲透性對(duì)湍流的影響，利用粒子圖像測(cè)速法(particle image velocimetry,簡(jiǎn)稱PIV)獲得了湍流流場(chǎng)的特性。何利民等[4-5]采用等動(dòng)量取樣與顯微照相相結(jié)合的方法研究了水平管中原油-水的液滴粒徑和壓降梯度與混合流量、含水率等參數(shù)及表面活性劑的關(guān)系。徐廣麗等[6]對(duì)管道中油水系統(tǒng)的特性進(jìn)行了研究，建立了水相厚度梯度的計(jì)算模型能很好地預(yù)測(cè)相界面分布。許偉偉等[7]研究了 U 形管路在內(nèi)流作用下的振動(dòng)特性，認(rèn)為流體速度增加，管道位移增加，管道最大應(yīng)力增加，管道基頻減小。

油液作為管路系統(tǒng)的傳輸介質(zhì)，其輸出的動(dòng)力學(xué)特征的研究方法有瞬態(tài)法和頻域法兩種。Sun等[8]建立了空氣彈簧連接管元模型，用瞬態(tài)法分析了模型的動(dòng)態(tài)特征，驗(yàn)證了所建模型的有效性。文獻(xiàn)[9]利用頻域法對(duì)光纖光學(xué)傳感器使用流動(dòng)引起的管道振動(dòng)進(jìn)行了實(shí)驗(yàn)研究，獲取了管道的振動(dòng)信號(hào)的振動(dòng)頻率特性。文獻(xiàn)[10-11]采用頻域法對(duì)可變節(jié)流閥-蓄能器子系統(tǒng)建立了數(shù)學(xué)模型，得出在最佳吸收效果時(shí)可變節(jié)流閥的開度與電機(jī)轉(zhuǎn)速的關(guān)系。柳貢民等[12]應(yīng)用傳遞矩陣方法研究了蒸汽參數(shù)對(duì)于蒸汽管路系統(tǒng)固有特性的影響，發(fā)現(xiàn)蒸汽壓力對(duì)于管路系統(tǒng)固有頻率的影響程度主要與管截面參數(shù)有關(guān)，而流速的影響相對(duì)不大。賀尚紅等[13]為解決液壓系統(tǒng)中由于壓力脈動(dòng)而引起的振動(dòng)和噪聲問題，基于管路動(dòng)態(tài)特性，建立薄板振動(dòng)式液壓脈動(dòng)衰減器的傳遞矩陣模型，對(duì)壓力脈動(dòng)的衰減特性進(jìn)行仿真。母東杰等[14]針對(duì)伺服閥控液壓管路系統(tǒng)內(nèi)部流體振蕩問題，采用一維流體瞬變理論分析了動(dòng)態(tài)過程中油液壓縮性對(duì)油液動(dòng)量的改變，建立了閥控液壓管路系統(tǒng)動(dòng)態(tài)數(shù)學(xué)模型。前期的這些研究都取得了一定的成果，但是很少涉及到油中水分含量對(duì)系統(tǒng)特性的影響，因而本課題利用分布參數(shù)理論建立油水乳化液的傳輸頻率模型，深入分析含水率對(duì)油水乳化液傳輸?shù)膭?dòng)態(tài)特性的影響關(guān)系，為有效檢測(cè)油中含水率提供科學(xué)依據(jù)。

1 油水乳化液的壓力波傳遞速度

油水乳化液傳輸?shù)膭?dòng)態(tài)特性與油水乳化液的壓縮性和質(zhì)量密切相關(guān)，這二者的相互關(guān)系又與壓力波速相關(guān)，因而首先分析油水乳化液中含水率與其傳輸?shù)膲毫Σㄋ俅笮〉年P(guān)系。

取油水乳化液的一段端面為A的管體如圖1所示，0-0斷面為乳化液靜止時(shí)所在的位置，當(dāng)乳化液在壓力P作用下以速度Vf向右運(yùn)動(dòng)dt時(shí)間后到達(dá)1-1斷面。假定管體是彈性體，又由于乳化液可壓縮的特性，則乳化液在dt時(shí)間內(nèi)由0-0斷面移動(dòng)到1-1斷面的過程中，1-1斷面的右側(cè)乳化液受到壓縮，也引起管體膨脹，使管體截面積由A增大了dA，乳化液的壓力增加了dP，從而產(chǎn)生了壓力波以速度cf向右傳播，在dt時(shí)間內(nèi)到達(dá)2-2斷面。

圖1 油水乳化液的一段斷面為A的管體Fig.1 A section of pipe body A in oil-water emulsion

定義油水乳化液中水分含量為cw，ρo，ρw分別為油液和水分的密度。根據(jù)連續(xù)性方程可得1-1斷面和2-2斷面之間的乳化液質(zhì)量守恒方程

ρwcwAVfdt+ρo(1-cw)AVfdt=

(dρwA+ρwdA)cw(cf-Vf)dt+

(dρoA+ρodA)(1-cw)(cf-Vf)dt?

ρwcw+ρo(1-cw)= [(dρwA+ρwdA)cw+

(1)

油水乳化液在壓力P作用下，從0-0斷面到2-2斷面內(nèi)的乳化液在dt時(shí)間內(nèi)從靜止變?yōu)樗俣萔f的運(yùn)動(dòng)狀態(tài)，則有

[ρwcw+ρo(1-cw)]cfVf=dp

(2)

(3)

定義Cw為油水乳化液的物性系數(shù)

可見，當(dāng)油液中含水率cw取不同值時(shí)，物性系數(shù)Cw隨之線性變化。如cw=0(或1)時(shí)，即流體為純油液(或純水)在彈性管道中壓力波水擊的波速表達(dá)式與Kortweg推導(dǎo)的公式一致；當(dāng)cw取值為0～1之間時(shí)，即含水率對(duì)油水乳化液的壓力波波速的影響較大。

2 油水乳化液的傳輸頻率特性模型

為了對(duì)油水乳化液進(jìn)行傳輸頻率特性分析，從乳化液的運(yùn)動(dòng)方程、連續(xù)性方程通過零初始條件的Laplace變換，由時(shí)域變換到s域，得到油水乳化液傳輸?shù)念l域描述。研究對(duì)象為油水真空分離裝置的進(jìn)油管路內(nèi)非恒定流乳化液，進(jìn)油管路系統(tǒng)示意圖如圖2所示。圖2中:①,⑤為聚氯乙烯(polyrinyl chloride,簡(jiǎn)稱PVC)透明鋼絲軟管;②,③,④為鋼管作為PVC透明鋼絲軟管的連接管件;⑥,⑧為鋼管分別與⑦流量傳感器相連接，⑨為真空室。油水乳化液在真空室內(nèi)真空壓力的作用下，從系統(tǒng)左端的油箱依次經(jīng)過各段管件進(jìn)入真空室，由于乳化液的流速遠(yuǎn)小于壓力波的傳播速度，且各段管件半徑小于油液流動(dòng)的波長(zhǎng)，因而假定管內(nèi)乳化液為無旋運(yùn)動(dòng)并且不受重力影響，為軸對(duì)稱流動(dòng)，壓力沿管橫斷面相等，沒有切向速度分量，不考慮乳化液和管壁間的熱傳導(dǎo)，即管壁是絕熱的，管內(nèi)油液與外界無熱交換，溫度在徑向均勻分布。

圖2 油水乳化液管路系統(tǒng)示意圖Fig.2 A diagram of oil-water emulsion pipe system

2.1 管件的頻率特性模型

(4)

2.2 建立油水乳化液管路系統(tǒng)的頻率特性模型

由于油水乳化液管路系統(tǒng)的動(dòng)力是真空室與油箱之間的壓差，因而將系統(tǒng)左端的油箱看作負(fù)載端，且是開端管路，系統(tǒng)右端的真空室⑨看作是系統(tǒng)的源元件，設(shè)Pz,Vf z分別為真空室⑨容腔的壓力、流速。則根據(jù)流體傳輸耦合理論，P(l(i+1)1,s)=P(li2,s)，Vf(l(i+1)1,s)=Vf(li2,s)，可得油水乳化液管路系統(tǒng)的頻率特性模型[16]。將s=jω代入油水乳化液管路系統(tǒng)的頻率特性模型，則油水乳化液管路系統(tǒng)的源阻抗頻率方程為

(5)

其中

3 油水乳化液傳輸?shù)膭?dòng)態(tài)特性分析

以油水真空分離裝置的進(jìn)油管路內(nèi)非恒定流乳化液為例，系統(tǒng)的結(jié)構(gòu)示意圖如圖2所示。系統(tǒng)的相關(guān)參數(shù)如下：①，⑤PVC透明鋼絲軟管的彈性模量為Kp=1.5×106Pa ，二者的管內(nèi)徑和壁厚一致，分別為D1=D5=0.026 m，δ1=δ5=0.003 m，二者的管長(zhǎng)分別為l1=4 m，l5=3.1 m；②,③,④,⑥,⑦,⑧鋼管的彈性模量Kc=1.9×1011Pa，其管路內(nèi)徑分別為D2=D4=0.02 m，D3=0.15 m，D6=D8=0.028 m，D7=0.006 m，各段管長(zhǎng)分別為l2=0.2 m，l3=0.75 m，l4=0.15 m，l6=0.4 m，l7=0.2 m，l8=2.05 m；油液密度ρo=870 kg/m3, 油液的體積彈性模量Ko=1.2×109Pa，運(yùn)動(dòng)粘度υo=4.6×10-5m2/s，水的密度ρw=1 000 kg/m3, 水的體積彈性模量Kw=2.2×109Pa，水的運(yùn)動(dòng)粘度υw=1.01×10-6m2/s。

3.1 水分含量對(duì)其壓力波傳播速度的影響

將相關(guān)參數(shù)代入式(3)中，并考慮流固(油水乳化液與管體)耦合作用與非流固耦合作用下，乳化液中含水率對(duì)其傳播速度的影響如圖3所示。

圖3 乳化液中含水率對(duì)其傳播速度的影響Fig.3 Effect of the rate of water content of emulsion on velocity

圖3(a)表示油水乳化液與PVC透明鋼絲軟管相互耦合作用下壓力波速隨含水率的變化規(guī)律。圖3(b)表示鋼管與油水乳化液在耦合作用下(虛線)和非耦合作用下(實(shí)線)壓力波速隨含水率的變化情況。由圖3可以看出，隨著乳化液中含水率的增大，PVC透明鋼絲軟管內(nèi)乳化液的傳播速度逐漸降低，當(dāng)cw=0時(shí)，即純油狀態(tài)下，其傳播速度為11.52 m/s，則當(dāng)cw=1時(shí)，即純水狀態(tài)下，其傳播速度為10.75 m/s，二者相差0.77 m/s。由于PVC透明鋼絲軟管的彈性模量遠(yuǎn)小于油和水的彈性模量，因而可將式(3)變換為

cf=

(6)

由于油的密度小于水的密度，所以考慮流固耦合作用的油水乳化液的傳播速度隨著含水率的增加而呈下降趨勢(shì)，如表1所示。這與不考慮流固耦合所得到的乳化液壓力波速(如圖3(b))相比，差距很大，因而對(duì)PVC透明鋼絲軟管內(nèi)乳化液的傳播速度的計(jì)算不能忽略流固耦合作用的影響。

表1乳化液中含水率與壓力波速的關(guān)系

Tab.1Therelationofwatercontentandpressurewavevelocityinemulsion

壓力波速/(m·s-1)含水率/%020406080100PVC透明鋼絲軟管(耦合)11.5211.3511.1911.0310.8810.75鋼管(耦合)113911801226127813361404非耦合作用117412211273133314021483

由圖3(b)可以看出，對(duì)于鋼管內(nèi)油水乳化液的傳播速度考慮耦合作用與非耦合作用下，都隨含水率的增加而增加，這是由于鋼管的彈性模量遠(yuǎn)大于油和水的彈性模量，則式(3)可變換為

cf=

(7)

當(dāng)cw=0時(shí)，即純油狀態(tài)下，其傳播速度考慮耦合作用時(shí)為1 139 m/s，非耦合作用為1 174 m/s，則當(dāng)cw=1時(shí)，即純水狀態(tài)下，其傳播速度考慮耦合作用時(shí)為1 483 m/s，非耦合作用為1 404 m/s。可見含水率對(duì)鋼管內(nèi)乳化液的傳播速度的影響很大，考慮耦合作用下二者相差265 m/s，不考慮耦合作用下二者相差309 m/s；而對(duì)于考慮耦合作用與不考慮耦合作用時(shí)，純油狀態(tài)下二者相差35 m/s，純水狀態(tài)下二者相差79 m/s。對(duì)于鋼管內(nèi)油水乳化液的速度計(jì)算，考慮含水率的影響是非常必要的。

3.2 水分含量對(duì)其傳輸頻率特性的影響

將油水真空分離裝置的進(jìn)油管路的真空室作為源端，油箱為開端管路，設(shè)置油水乳化液的含水率分別為0%，20%，40%，60%，80%，100%，將其余相關(guān)參數(shù)代入式(5)中，利用Matlab編寫求解程序，從而繪出系統(tǒng)源阻抗頻率特性曲線如圖4所示。

由圖4可知，油水乳化液管路系統(tǒng)的同階諧振頻率隨著含水率的增加而呈遞減趨勢(shì)，且階次越大，這種差異性就越明顯。從表2可以看到，在含水率為40%時(shí)所建模型獲取了油水乳化液管路系統(tǒng)的基頻為0.98 Hz，在含水率為0%，20%，60%，80%時(shí)獲取系統(tǒng)的頻率接近系統(tǒng)基頻3倍頻為2.93 Hz，100%時(shí)獲取系統(tǒng)的頻率接近系統(tǒng)基頻2倍頻為1.95 Hz；不同含水率的乳化液在前6階諧振頻率中有相互交叉，不易分辨的特點(diǎn)，而高于6階后其諧振頻率明顯相差1 Hz，這是由于本項(xiàng)研究對(duì)象的管路系統(tǒng)總長(zhǎng)為10.85 m，而兩段PVC透明鋼絲軟管共長(zhǎng)7.1 m，鋼管共長(zhǎng)3.75 m，油水乳化液在本系統(tǒng)中的傳播速度主要受到PVC透明鋼絲軟管的耦合作用影響，而純油與純水在該管路中的傳播速度差別不大，但含水率越大其傳播速度越小，因而從總體上看系統(tǒng)的諧振頻率會(huì)隨含水率增加而減低。

表2不同含水率對(duì)應(yīng)的各階諧振頻率

Tab.2Theresonancefrequencyofdifferentwatercontentcorresponding

各階頻率/Hz含水率/%02040608010012.932.930.982.932.931.95211.7211.7210.74322.4622.4622.4622.46447.8545.90570.3161.5257.6159.57679.1071.4666.40797.6689.848117.19115.23114.26110.3591.799278.32274.41180.66146.48144.5310218.7511373.05346.68238.2812443.36

3.3 實(shí)驗(yàn)驗(yàn)證與分析

分別配置不同的油液含水濃度的46#汽輪機(jī)油樣，攪拌30 min后，根據(jù)不同的實(shí)驗(yàn)條件(如表3所示)進(jìn)行實(shí)驗(yàn)，開啟真空泵進(jìn)入一個(gè)工作循環(huán)。該實(shí)驗(yàn)系統(tǒng)的測(cè)試部分由進(jìn)液管的流量傳感器、溫度傳感器、真空濾油機(jī)內(nèi)真空傳感器以及配套的多路傳感器，帶數(shù)據(jù)采集卡的PXI測(cè)試機(jī)(America，NI)等一系列硬件及與測(cè)試功能相匹配的軟件共同構(gòu)成。其中，在圖2中PT為高精度壓力傳感器，將采集的實(shí)時(shí)乳化液的傳輸動(dòng)態(tài)信號(hào)輸出給NI-PCI 6221采集卡的一體機(jī)，經(jīng)過數(shù)據(jù)分析模塊(數(shù)據(jù)濾波、Hilbert變換分析和數(shù)據(jù)計(jì)算等)獲取了不同含水率的乳化液的壓力波動(dòng)態(tài)信號(hào)如圖5所示。

表3 實(shí)驗(yàn)參數(shù)取值范圍

含水油液的諧振信號(hào)是一種典型的調(diào)制信號(hào)，而含水油液的諧振信號(hào)的載波頻率(如壓力傳感器動(dòng)態(tài)響應(yīng)及測(cè)試系統(tǒng)固有頻率)通常不易獲取，因而采用希爾伯特(Hilbert)變換方法進(jìn)行包絡(luò)解調(diào)。Hilbert變換主要由2個(gè)部分組成：經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition，簡(jiǎn)稱EMD)和Hilbert譜分析。EMD是一種自適應(yīng)的、高效的數(shù)據(jù)分解方法，是以局部時(shí)間尺度為基礎(chǔ)，將任何復(fù)雜的數(shù)據(jù)集被分解為個(gè)數(shù)有限的幾個(gè)固有模函數(shù)(intrinsic mode functions，簡(jiǎn)稱IMF)的線性疊加。在此基礎(chǔ)上，對(duì)各個(gè)IMF進(jìn)行Hilbert變換，取模后得到的就為含水油液的調(diào)制信號(hào)，就可以有效地將壓力傳感器動(dòng)態(tài)響應(yīng)及測(cè)試系統(tǒng)固有頻率即載波頻率去除。

圖5(a)，(c)是含水率分別為20%,40%的乳化液經(jīng)過Hilbert包絡(luò)處理后的信號(hào)，圖5(b)，(d)是對(duì)含水率分別為20%，40%乳化液的柏拉圖經(jīng)過傅里葉變換后的頻譜，含水率為20%，40%乳化液的各階諧振頻率實(shí)測(cè)值與理論值對(duì)比如表4所示。

圖5 不同含水率的乳化液的壓力波動(dòng)態(tài)信號(hào)Fig.5 The pressure wave dynamic signal of different water content of emulsion

Tab.4Thecontrastofresonancefrequencyofdifferentwatercontent

各階頻率/Hz含水率/%20(理論)20(實(shí)測(cè))40(理論)40(實(shí)測(cè))12.931.950.980.97211.7212.6710.747.80322.4624.3622.4621.44447.8549.7149.70561.5264.3257.6155.55679.1077.9771.4670.17797.66100.3889.8487.818115.23115.00114.26116.219274.41180.66184.5710218.75221.6811346.68

可見與20%,40%含水率的油水乳化液的各階諧振頻率理論值和實(shí)測(cè)值很接近，因而所建立的油水乳化液傳輸頻率模型能夠反映出含水率對(duì)其動(dòng)態(tài)特性的影響，為進(jìn)一步根據(jù)管體的頻率特性獲取乳化液中含水率的大小提供了一種新的途徑。

4 結(jié) 論

1) 推導(dǎo)了油水乳化液在管體中壓力波傳播速度公式，分析了含水率對(duì)乳化液傳輸速度的影響規(guī)律。表明在PVC透明鋼絲軟管中的乳化液有必要考慮流固耦合作用，含水率增加，其傳播速度呈降低趨勢(shì)；反之，含水率對(duì)于鋼管內(nèi)乳化液的傳輸速度的影響較考慮耦合作用與不考慮耦合作用的影響大，且隨著含水率的增加，其傳播速度呈遞增趨勢(shì)。

2) 以油水真空分離裝置的進(jìn)油管路內(nèi)非恒定流乳化液為例，考慮了PVC透明鋼絲軟管的耦合作用，利用分布參數(shù)原理以及傳遞矩陣法建立了乳化液傳輸?shù)念l率特性方程，獲得了乳化液中含水率與系統(tǒng)動(dòng)態(tài)特性之間的非線性關(guān)系。

3) 結(jié)合實(shí)際工程，設(shè)置了乳化液系統(tǒng)的相關(guān)物理參數(shù)，對(duì)含水率與乳化液傳輸系統(tǒng)的動(dòng)態(tài)特性之間的關(guān)系進(jìn)行了仿真，獲得了含水率越大，系統(tǒng)的同階諧振頻率越低，且隨著含水率增加，諧振頻率的階次越大，同階的諧振頻率之間的差值也相應(yīng)變大。

4) 為了驗(yàn)證所建模型的有效性，對(duì)油水乳化液的管路系統(tǒng)進(jìn)行了驗(yàn)證，通過采集含水率為20%，40%的乳化液的壓力波實(shí)測(cè)信號(hào)，經(jīng)過濾波/Hilbert包絡(luò)分解等處理，得到的諧振頻率與仿真的數(shù)據(jù)較為吻合，為進(jìn)一步根據(jù)管體的頻率特性獲取乳化液中含水率的大小提供了一種新的途徑。

[1] Alinezhad K, Hosseini M, Movagarnejad K,et al. Experimental and modeling approach to study separation of water in crude oil emulsion under non-uniform electrical field[J].Korean Journal of Chemical Engineering, 2010, 27(1): 198-205.

[2] Kang J H, Lee K J, Nam J H,et al. Visualization of invasion percolation drainage process in porous media using density-matched immiscible fluids and refractive index matched solid structures[J].Journal of Power Sources, 2010, 195(9):2608-2612.

[3] Suga K Y, Matsumura A,Shitaka Y,et al. Effects of wall permeability on turbulence[J].International Journal of Heat and Fluid Flow, 2010, 31(6):974-984.

[4] 呂宇玲，何利民，丁慎圓，等.原油-水分散流中液滴粒徑及壓降規(guī)律研究[J].工程熱物理學(xué)報(bào)，2013，34(3)：472-475.

Lv Yuling,He Limin,Ding Shenyuan,et al.Investigation of droplet size and pressure drop in crude oil-water dispersed flows [J]. Journal of Engineering Thermophysics, 2013,34(3):472-475.(in Chinese)

[5] 呂宇玲，何利民，程浩，等.油水兩相分散流液滴粒徑預(yù)測(cè)模型[J].化工學(xué)報(bào)，2012，63(7)：2059-2063.

Lü Yuling,He Limin,Cheng Hao,et al. Oil and water two-phase flow droplet size dispersion prediction model[J]. Journal of Chemical Industry and Engineering, 2012,63(7):2059-2063.(in Chinese)

[6] 徐廣麗，張國(guó)忠，Brauner N，等.管道油流攜水系統(tǒng)的界面分布[J].中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版，2011，35(1)：124-129.

Xu Guangli,Zhang Guozhong, Brauner N,et al. Pipeline carrying oil flow water distribution system interface [J]. Journal of China University of Petroleum：Natural Science Edition, 2011, 35(1):124-129.(in Chinese)

[7] Xu Weiwei, Wu Bo, Wu Turn, et al, Unsteady fluid and the U-shaped pipe coupling vibration characteristics [J]. Journal of Chemical Engineering, 2012, 26(5):770-774.

[8] Sun Liqing, Li Zhongxing, Shen Xufeng, et al. Simulation and test study on dynamic characteristic of air spring with auxiliary chamber[J]. Applied Mechanics & Materials, 2013,341-342(341-342):391-394.

[9] Tang Xiaochun, Zhong Wei, Zhuang Dachang, et al. Study of optical fiber fluid flow monitoring system using pipe vibration frequency characteristic [J]. Applied Mechanics and Materials, 2013, 295(2):2054-2057.

[10] 廖健，何琳，祁曉野，等.一種變頻率流量脈動(dòng)吸收方法[J].艦船科學(xué)技術(shù)，2013，35(2)：39-42.

Liao Jian, He Lin , Qi Xiaoye, et al. A variable flow pulsation frequency absorption method[J].Ship Science and Technology,2013,35(2):39-42.(in Chinese)

[11] 蔣祺， 張方，姜金輝，等.動(dòng)態(tài)載荷位置識(shí)別當(dāng)量載荷誤差判別法[J]. 振動(dòng)、測(cè)試與診斷，2015，35(4)：620-625.

Jiang Qi, Zhang Fang, Jiang Jinhui,et al.Identification of equivalent load error for dynamic load position identification[J]. Journal of Vibration ,Measurement & Diagnosis ,2015,35(4):620-625.(in Chinese)

[12] Liu Gongmin, Chen Hao, Li Shuaijun. Steam piping system parameters on the dynamic characteristics of research[J]. Journal of Solid Mechanics,2012,33(2):168-175.

[13] 賀尚紅，王雪芝，何志勇，等.薄板振動(dòng)式液壓脈動(dòng)衰減器濾波特性[J].機(jī)械工程學(xué)報(bào)，2013，49(4)：148-153.

He Shanghong,Wang Xuezhi, He Zhiyong, et al. Vibrating plate type hydraulic damper filtering characteristics[J].Journal of Mechanical Engineering,2013,49(4):148-153.(in Chinese)

[14] 母東杰，李長(zhǎng)春，延皓，等.基于特征線理論的閥控液壓管路瞬變過渡流數(shù)值分析[J].兵工學(xué)報(bào)，2012，33(12)：1455-1460.

Mu Dongjie, Li Changchun , Yan Hao,et al. SUN meng feature-based line theory transitional flow valve hydraulic line transients numerical analysis[J].Acta Armamentarii ,2012,33(12):1455-1460 .(in Chinese)

[15] Chen Bin, Yi Menglin. A mathematics model investigation of ideal fluids in a pipeline of moving ending[J].Journal of Harbin Institute of Technology,2007,14(1):259-262.

[16] 陳彬，劉閣，張賢明.基于動(dòng)邊界的閥控缸系統(tǒng)水擊振動(dòng)傳輸特性分析[J].機(jī)械科學(xué)與技術(shù)，2010，29(1)：123-128.

Chen Bin,Liu Ge,Zhang Xianming. Based on dynamic boundary valve cylinder system vibration transmission characteristics of water hammer [J]. Mechanical Science and Technology for Aerospace Engineering, 2010,29(1):123-128.(in Chinese)

國(guó)家自然科學(xué)基金資助項(xiàng)目(51375516);重慶基礎(chǔ)與前沿研究資助項(xiàng)目(cstc2016jcyjA0185)

2015-04-27；

2015-10-15

TH13

陳彬,男，1972年9月生，博士、教授。主要研究方向?yàn)橛鸵何廴究刂萍夹g(shù)。曾發(fā)表《基于HHT包絡(luò)譜的油中水分含量檢測(cè)研究》(《儀器儀表學(xué)報(bào)》2014年第35卷第2期)等論文。E-mail：hustchb@163.com