張 媛 王碩亮 張賢松

(1. 中國石化石油勘探開發(fā)研究院 北京 100083； 2. 中國地質(zhì)大學(xué)(北京)能源學(xué)院 北京 100083； 3. 中海油研究總院 北京 100028)

微圓管實驗壓力誤差校正方法探討*

張 媛1,2王碩亮2張賢松3

(1. 中國石化石油勘探開發(fā)研究院 北京 100083； 2. 中國地質(zhì)大學(xué)(北京)能源學(xué)院 北京 100083； 3. 中海油研究總院 北京 100028)

目前微圓管實驗的流量計量管中氣液兩相毛管力采用靜態(tài)毛管力公式計算，沒有考慮流體流動速度差異造成的毛管力動態(tài)變化機理。本文提出了采用數(shù)值方法計算動態(tài)毛管力來校正微圓管實驗壓力誤差的新方法。該方法針對微圓管實驗裝置的幾何特征，建立了數(shù)值模擬模型，利用Level-set方法計算兩相界面處的流體流動，得到了不同流動速度條件下的動態(tài)毛管力。對比表明，隨著流速的增加，動態(tài)毛管力數(shù)值增加的幅度降低，動態(tài)毛管力數(shù)值與靜態(tài)毛管力數(shù)值差別增大。利用動態(tài)毛管力校正了微圓管實驗結(jié)果，校正后的實驗數(shù)據(jù)與校正前的規(guī)律一致，但是數(shù)值存在明顯偏差。為得到更準確的結(jié)果，建議利用本文方法校正微圓管實驗的驅(qū)替壓差。

微圓管；實驗壓力誤差；動態(tài)毛管力；壓差校正

由于低、特低滲透油氣田儲層的孔道半徑處于微米甚至納米級別，導(dǎo)致經(jīng)典的滲流模型無法描述致密儲層中流體的滲流特征[1-2]，因此微圓管實驗(是研究微觀流體滲流規(guī)律的基礎(chǔ)實驗)得到了廣泛應(yīng)用[3-4]。微圓管實驗測試目標是微圓管兩端的壓差和流速，由于物質(zhì)守恒原理，流量測量的數(shù)值相對準確，而壓力測量存在系統(tǒng)誤差。目前大部分學(xué)者認為微圓管兩端的壓力差為壓力計顯示數(shù)值與大氣壓之間的差值[5]。王斐[6]指出，壓力計與大氣壓之間的壓力差值并不能表征微圓管兩端的壓差，中間存在系統(tǒng)誤差，即由于管壁粗糙形成的沿程壓力損失、管徑突然擴大縮小形成的局部壓力損失、流速計量管中存在兩相(液相、氣相)而形成的毛管力壓力損失；并對各系統(tǒng)誤差進行了計算，發(fā)現(xiàn)沿程壓力損失和局部壓力損失造成的誤差可以忽略，而毛管力壓力損失造成的壓力誤差約有4%(管徑為320 μm)，不可忽略，并且隨著計量管管徑的縮小，計量誤差會增加；進而提出在毛管力計算過程中，采用經(jīng)典靜態(tài)毛管力計算公式進行校正。但整個微圓管實驗過程中驅(qū)替流體沿著微圓管流動，流速計量管中的氣液界面沿著流動方向發(fā)生位移，管中的毛管力應(yīng)該為動態(tài)毛管力，因此用靜態(tài)毛管力公式的計算結(jié)果進行實驗壓力誤差校正必然會影響實驗結(jié)果的準確性。筆者提出采用數(shù)值方法計算動態(tài)毛管力來對實驗壓力誤差進行校正。

1 微圓管實驗壓力誤差計算新方法

動態(tài)毛管力在水動力學(xué)領(lǐng)域研究較多[7-9]，但測試得到的數(shù)據(jù)和建立的數(shù)學(xué)模型[10]無法直接應(yīng)用于微圓管實驗中的壓力誤差校正。要計算微圓管中的動態(tài)毛管力，必須要選擇一種合適的微孔道計算流體力學(xué)方法。目前研究兩相流動的計算流體力學(xué)方法可以分為3類：宏觀方法、微觀方法[11]和介觀方法[12]。本文針對微圓管中氣液兩相驅(qū)替開展研究，由于液體分子半徑自由程較小，流量計量管中的孔道半徑較大(一般為300 μm)，因此，流量計量管內(nèi)部的液體流動符合達西定律，可用宏觀方法開展計算。Level-set方法是最為常用的研究兩相流動的計算方法[13]，也是現(xiàn)今在利用網(wǎng)格數(shù)值模擬方法解決自由表面流動問題中較為流行的方法[14-16]，該方法在解決多尺度氣液兩相流動過程中的自由界面追蹤問題已較為成熟[17]。因此，本文選用Level-set方法計算兩相界面處的流體流動。

1.1 Level-set方法原理

圖1 計算區(qū)域物質(zhì)的界面

(1)

用尋找φ零等值線的方法確定兩相界面的位置Г(t)。

(2)

(3)

(4)

(5)

式(5)即為Level-set方程。目前求解Level-set方程的方法主要是基于有限差分法的高分辨率的差分格式，這些方法[18-19]具有高精度的優(yōu)點，而且已非常成熟。

1.2 基于Level-set方法的數(shù)學(xué)模型建立

利用Level-set方法可以計算兩相驅(qū)替過程中接觸界面的變化情況，分別統(tǒng)計兩相接觸面兩側(cè)的壓力，即可得到不同條件下的動態(tài)毛管力。采用文獻[20]中微觀實驗裝置尺寸，根據(jù)物理模擬中流量計量管的尺度，建立了數(shù)值模擬模型，模型參數(shù)取值見表1，模型網(wǎng)格剖分結(jié)果見圖2。

表1 數(shù)值模擬模型參數(shù)賦值表

圖2 數(shù)值模擬模型的網(wǎng)格剖分結(jié)果

數(shù)值模型入口端為恒定壓力邊界，壓力為0.15 MPa。模型的出口端與大氣相通，壓力恒定為0.1 MPa。根據(jù)物理模型中流量計量管的兩相分布狀態(tài)，數(shù)值模擬模型中流量計量管的前部為液體，流量計量管中部存在一個氣泡。

2 微圓管實驗驅(qū)替壓差校正

利用Level-set有限元方法，選取文獻[20]中不同微圓管半徑的實驗方案(流量計量管尺寸相同，不同方案流速不同)進行計算，得到了不同驅(qū)替速度條件下的動態(tài)毛管力(圖3)。

從圖3可以看出，動態(tài)毛管力的數(shù)值均高于靜態(tài)毛管力，這與水動力學(xué)領(lǐng)域(氣水動態(tài)毛管力)的實驗結(jié)果規(guī)律性一致[21-22]。隨著流速的增加，動態(tài)毛管力數(shù)值增加的幅度降低。利用動態(tài)毛管力的計算結(jié)果，對實驗數(shù)據(jù)進行校正，將實驗測試得到的驅(qū)替壓差減去對應(yīng)方案的動態(tài)毛管力，得到了微圓管兩端的真實壓差，并且計算了是否利用動態(tài)毛管力進行實驗數(shù)據(jù)校正的誤差，具體結(jié)果見表2。

圖3 利用Level-set計算動態(tài)毛管力結(jié)果

表2 動態(tài)毛管力計算結(jié)果與實驗誤差分析

注：誤差*=[1-(實驗測試壓差-動態(tài)毛管力)/實驗測試壓差]×100%。 從表2可以看出，由于流量計量管的尺寸不變，動態(tài)毛管力的數(shù)值僅與實驗流體的流動速度有關(guān)，實驗流體的流動速度越大，對應(yīng)的壓差越大，動態(tài)毛管力數(shù)值也越大。隨著實驗流體流速的增加，動態(tài)毛管力的變化速度小于微圓管兩端壓差的變化速度，因此，實驗流體流速越快，誤差越小。但是，在流速較低時，誤差最大可以達到85%，有必要利用動態(tài)毛管力對實驗數(shù)據(jù)進行校正。

經(jīng)過動態(tài)毛管力校正后的實驗數(shù)據(jù)如圖4所示?？梢钥闯觯Ｕ蟮膶嶒灁?shù)據(jù)與校正前的實驗數(shù)據(jù)規(guī)律上保持一致，但是數(shù)值存在明顯偏差，在進行實驗數(shù)據(jù)分析和滲流模型建立時，若不進行動態(tài)毛管力校正，會得出錯誤的結(jié)果。

圖4 動態(tài)毛管力校正前后結(jié)果對比

3 結(jié)論

提出了利用動態(tài)毛管力進行微圓管實驗壓力誤差校正的思路方法，根據(jù)文獻[20]中實驗物理模型實際尺寸和流體性質(zhì)，建立了動態(tài)毛管力計算模型，利用Level-set計算方法，計算得到了不同流動速度對應(yīng)的動態(tài)毛管力。計算結(jié)果表明，動態(tài)毛管力的數(shù)值均高于靜態(tài)毛管力，隨著流速的增加，動態(tài)毛管力數(shù)值增加的幅度降低，動態(tài)毛管力數(shù)值與靜態(tài)毛管力數(shù)值差別增大。利用動態(tài)毛管力對微圓管實驗的驅(qū)替壓差進行了校正，校正后的實驗數(shù)據(jù)與校正前的實驗數(shù)據(jù)規(guī)律上保持一致，但是數(shù)值存在明顯偏差，建議應(yīng)采用本文提出的方法對微圓管實驗的驅(qū)替壓差進行校正。

[1] BIKWEMU R,WOLFE A J,XING X,et al.Facilitated translocation of polypeptides through a single nanopore[J].Journal of Physics Condensed Matter,2010,22(45):2889-2898.

[2] JAVADPOUR F,FISHER D,UNSWORTH M.Nanoscale gas flow in shale gas sediments[J].Journal of Canadian Petroleum Technology,2007,46(10):55-61.

[3] 王小鋒,朱維耀,鄧慶軍,等.考慮固液范德華力作用的微圓管流動數(shù)學(xué)模型[J].東北石油大學(xué)學(xué)報,2013,37(5):85-89.

Wang Xiaofeng,Zhu Weiyao,Deng Qingjun,et al.Micro circular pipe flow mathematical model considering the effect of Van der Waals force[J].Journal of Northeast Petroleum University,2013,37(5):85-89.

[4] 趙清民,岳湘安,王斐,等.高壓氣體在微圓管中的流動特性[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2012,36(4):129-133.

Zhao Qingmin,Yue Xiang’an,Wang Fei,et al.Flow characteristics of high pressure real gas in microtubes [J].Journal of China University of Petroleum(Edition of Natural Science),2012,36(4):129-133.

[5] 朱維耀,田英愛,于明旭,等.微圓管中流體的微觀流動機制[J].科技導(dǎo)報,2014,32(27):23-27.

Zhu Weiyao,Tian Ying’ai,Yu Mingxu,et al.Mechanism of microscopic fluid flow in microtubes[J].Science and Technology Review,2014,32(27):23-27.

[6] 王斐,岳湘安,徐紹良,等.潤濕性對水在微管和巖芯中流動特性的影響[J].科學(xué)通報,2009,54(7):972-977.

Wang Fei,Yue Xiang’an,Xu Shaoliang,et al.Influence of wettability on flow characteristics of water through microtubes and cores[J].Chinese Science Bulletin,2009,54(7):972-977.

[7] WEITZ D A,STOKES J P,BALL R C,et al.Dynamic capillary pressure in porous media:origin of the viscous-fingering length scale[J].Physical Review Letters,1988,59(26):2967-2970.

[8] HASSANIZADEH S M,CELIA M A,DAHLE H K.Dynamic effect in the capillary pressure-saturation relationship and its impacts on unsaturated flow[J].Vadose Zone Journal， 2002,1(1):38-57.

[9] BENTSEN R G.Effect of hydrodynamic forces on capillary pressure and relative permeability[J].Transport in Porous Media， 1994,17(2):121-132.

[10] CAMPS-ROACH G,O'CARROLL D M,NEWSON T A,et al.Experimental investigation of dynamic effects in capillary pressure:Grain size dependency and upscaling[J].Water Resources Research，2010,46(8):863-864.

[11] 任美鵬,李相方,劉書杰,等.新型深水鉆井井噴失控海底搶險裝置概念設(shè)計及方案研究[J].中國海上油氣,2014,26(2):66-71.

Ren Meipeng,Li Xiangfang,Liu Shujie,et al.Research on the conceptual design of new seabed rescue equipment for uncontrolled blowout in deep water drilling [J].China Offshore Oil and Gas,2014,26(2):66-71.

[12] AIDUN C K,CLAUSEN J R.Lattice-Boltzmann method for complex flows[J].Annual Review of Fluid Mechanics，2009,42(1):439-472.

[13] 張超,張智,曾春珉,等.潿洲11-4油田含CO2氣井油管柱腐蝕分析[J].中國海上油氣,2015,27(4):122-125.

Zhang Chao,Zhang Zhi,Zeng Chunmin,et al.Analysis on tubing corrosion for gas wells with CO2in WZ 11-4 oilfield[J].China Offshore Oil and Gas,2015,27(4):122-125.

[14] 徐喬威,胡志強,謝彬,等.LNG運輸船旁靠FLNG卸載作業(yè)時的水動力性能試驗[J].中國海上油氣,2015,27(2):112-119.

Xu Qiaowei,Hu Zhiqiang,Xie Bin,et al.Hydrodynamics of an FLNG system in side-by-side offloading operation with a LNG carrier[J].China Offshore Oil and Gas,2015,27(2):112-119.

[15] 王成官,王嘉松,田中旭,等.不同預(yù)緊力時隔水管渦激振動特性三維數(shù)值模擬研究[J].中國海上油氣,2011,23(6):415-419.

Wang Chengguan,Wang Jiasong,Tian Zhongxu,et al.Three-dimensional numerical simulation of marine riser VIV characteristics with different pre-tightening force[J].China Offshore Oil and Gas,2011,23(6):415-419.

[16] 陳欣,孫旭,李東芳,等.海上平臺直升機甲板受環(huán)境影響的安全分析評估[J].中國海上油氣,2012,24(1):60-64.

Chen Xin,Sun Xu,Li Dongfang,et al.Safety assessment of the environmental impact on helideck of offshore platform[J].China Offshore Oil and Gas,2012,24(1):60-64.

[17] SUSSMAN M,FATEMI E,SMEREKA P,et al.An improved level set method for incompressible two-phase flows[J].Computers and Fluids,1998,27(5):663-680.

[18] LUO J,LUO Z,CHEN L,et al.A semi-implicit level set method for structural shape and topology optimization[J].Journal of Computational Physics,2008,227(11):5561-5581.

[19] KEES C E,AKKERMAN I,FARTHING M W,et al.A conservative level set method suitable for variable-order approximations and unstructured meshes[J].Journal of Computational Physics,2011,230(12):4536-4558.

[20] 張雪齡.考慮液—固界面作用的微可壓縮流體的滲流理論研究[D].北京：北京科技大學(xué),2015.

Zhang Xueling.Percolation theory research of weakly compressible fluid flow considering wall-liquid interaction[D].Beijing:University of Science and Technology Beijing,2015.

[21] ABIDOYE L K,DAS D B.Scale dependent dynamic capillary pressure effect for two-phase flow in porous media[J].Advances in Water Resources,2014,74(1):212-230.

[22] CUESTA C,DUIJN C J V,HULSHOF J.Infiltration in porous media with dynamic capillary pressure:travelling waves[J].European Journal of Applied Mathematics,2000,11(4):397-398.

(編輯：楊 濱)

Discussion on method for correcting the pressure error in micro tube experiment

Zhang Yuan1, 2Wang Shuoliang2Zhang Xiansong3

(1.SinopecPetroleumExploration&ProductionResearchInstitute,Beijing100083,China;2.SchoolofEnergyResources,ChinaUniversityofGeosciences,Beijing100083,China;3.CNOOCResearchInstitute,Beijing100028,China)

The gas-liquid two-phase capillary force is calculated by static capillary force formula in micro tube experiment, without considering the dynamic capillary force mechanism caused by the difference of flow velocity. This paper presents the numerical calculation method of dynamic capillary force to correct micro tube experiment pressure error. The numerical simulation model is established based on the geometric characteristics of the micro tube，the Level-set method is used to calculate the fluid flow at the two phase interface，and the dynamic capillary force under different flow velocity is obtained. Results show that with the increase of flow velocity, the increasing degree of the dynamic capillary force reduce, the difference between dynamic capillary force and static capillary force increases. The results of micro tube experiment are corrected using dynamic capillary force，the corrected experimental data are in agreement with the law before the correction, but the numerical value has obvious deviation. In order to get more accurate results, it is recommended to correct the displacement pressure of micro tube experiment using this method.

micro tube; experimental pressure error; dynamic capillary force; pressure difference correction

*國家自然科學(xué)基金“疏松砂巖儲層深部調(diào)剖堵劑分布模式研究(編號：51504223)”、基本科研優(yōu)秀教師資助基金“竄流通道識別與描述方法研究(編號：53200859060)”部分研究成果。

張媛，女，工程師，2008年畢業(yè)于中國石油大學(xué)(北京)，并獲油氣田開發(fā)工程碩士學(xué)位，目前主要從事油藏工程方面的工作。地址：北京市海淀區(qū)學(xué)院路31號(郵編：100083)。E-mail：zhangyuan-101@163.com。

王碩亮，男，講師， 2011年博士畢業(yè)于中國石油大學(xué)(北京)，主要從事油氣田開發(fā)方面的研究。地址：北京市海淀區(qū)學(xué)院路29號(郵編：100083)。 E-mail:wangshuoliang@cugb.edu.cn。

1673-1506(2016)02-0094-05

10.11935/j.issn.1673-1506.2016.02.012

TE 312

A

2015-11-16 改回日期：2015-12-20

張媛,王碩亮,張賢松.微圓管實驗壓力誤差校正方法探討[J].中國海上油氣，2016,28(2)：94-98.

Zhang Yuan,Wang Shuoliang,Zhang Xiansong.Discussion on method for correcting the pressure error in micro tube experiment[J].China Offshore Oil and Gas,2016,28(2):94-98.