高頻波動(dòng)對(duì)多孔彈性介質(zhì)中稠油黏溫參數(shù)影響規(guī)律研究

蒲春生，徐加祥，劉 璽，鄭黎明，劉 靜

(1.中國(guó)石油大學(xué)(華東) 石油工程學(xué)院，山東 青島 266580；2.中國(guó)石油大學(xué)(華東) 重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室，山東 青島 266580；3.延長(zhǎng)油田股份有限公司 西區(qū)采油廠，陜西 延安 717500)

?

高頻波動(dòng)對(duì)多孔彈性介質(zhì)中稠油黏溫參數(shù)影響規(guī)律研究

蒲春生1,2，徐加祥1,2，劉 璽3，鄭黎明1,2，劉 靜1,2

(1.中國(guó)石油大學(xué)(華東) 石油工程學(xué)院，山東 青島 266580；2.中國(guó)石油大學(xué)(華東) 重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室，山東 青島 266580；3.延長(zhǎng)油田股份有限公司 西區(qū)采油廠，陜西 延安 717500)

為了定量描述高頻波動(dòng)激勵(lì)井熱力學(xué)降黏解堵的機(jī)理，基于高頻波在多孔彈性介質(zhì)中能量傳播、損耗和轉(zhuǎn)化的動(dòng)力學(xué)特征，建立了高頻波動(dòng)激勵(lì)下稠油油藏近井帶儲(chǔ)層中稠油溫度及黏度變化的穩(wěn)態(tài)和非穩(wěn)態(tài)耦合動(dòng)力學(xué)模型。該模型由能量衰減和轉(zhuǎn)化方程、非均勻熱源圓筒壁導(dǎo)熱微分方程和溫黏狀態(tài)方程組成，采用耦合迭代方法進(jìn)行數(shù)值計(jì)算，揭示了高頻波對(duì)多孔彈性介質(zhì)中稠油溫度和黏度影響的主控因素及規(guī)律。結(jié)果表明：頻率為20 kHz的高頻波動(dòng)在r<0.6 m范圍內(nèi)作用效果明顯，頻率越高作用效果越好，但處理范圍越??；與波動(dòng)頻率相比，振幅對(duì)其作用效果更顯著；波動(dòng)作用由非穩(wěn)態(tài)向穩(wěn)態(tài)過(guò)渡所用時(shí)間約400 s。該模型的建立與求解，為稠油油藏近井地帶高頻波動(dòng)解堵增產(chǎn)增注技術(shù)提供了重要的理論依據(jù)。

稠油開(kāi)采；高頻波動(dòng)；降黏解堵

蒲春生，徐加祥，劉璽，等.高頻波動(dòng)對(duì)多孔彈性介質(zhì)中稠油黏溫參數(shù)影響規(guī)律研究[J].西安石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016，31(6)：54-59,107.

PU Chunsheng,XU Jiaxiang,LIU Xi,et al.Study on influence of high frequency vibration on viscosity and temperature parameters of heavy oil in porous elastic medium[J].Journal of Xi'an Shiyou University (Natural Science Edition),2016,31(6):54-59,107.

引 言

近幾十年來(lái)，高頻波動(dòng)采油技術(shù)在油田開(kāi)發(fā)過(guò)程中取得了比較好的礦場(chǎng)應(yīng)用效果，成為一種比較典型的物理法提高采收率技術(shù)。利用高頻波動(dòng)處理生產(chǎn)井和注水井的近井地帶，憑借其力學(xué)效應(yīng)、熱學(xué)效應(yīng)和化學(xué)效應(yīng)，可以解除地層堵塞，在油層產(chǎn)生微裂縫，改變流體的物性及流動(dòng)狀態(tài)[1]，從而改善近井地帶的滲流狀況，達(dá)到增產(chǎn)增注的目的。

稠油由于其黏度高、凝固點(diǎn)高的特點(diǎn)，礦場(chǎng)普遍采用注熱流體加熱的方式開(kāi)發(fā)。但是在此過(guò)程中，大量高溫流體改變了油藏內(nèi)部孔隙條件，特別是在近井地區(qū)引起瀝青質(zhì)沉積，堵塞孔喉，造成注氣壓力升高[2-3]。高頻波動(dòng)作為一種物理降黏解堵方法，通過(guò)波傳播過(guò)程中能量的衰減及向熱能的轉(zhuǎn)化，給稠油加熱升溫，降低其黏度，達(dá)到近井地帶降黏解堵的效果。由于飽和流體多孔彈性介質(zhì)的復(fù)雜性，前人的研究普遍局限于室內(nèi)實(shí)驗(yàn)，對(duì)其熱力學(xué)降黏機(jī)理大多停留在定性解釋的層面[4-7]，對(duì)高頻波動(dòng)熱力學(xué)降黏規(guī)律及降黏模型還未形成成熟的理論。

本文從高頻波動(dòng)能量由機(jī)械能向熱能轉(zhuǎn)化的過(guò)程出發(fā)，將波動(dòng)能量衰減和轉(zhuǎn)化方程、非均勻熱源圓筒壁導(dǎo)熱微分方程和溫黏狀態(tài)方程組合，建立了定量表征高頻波參數(shù)對(duì)多孔彈性介質(zhì)中稠油黏度和溫度影響的耦合動(dòng)力學(xué)模型。采用耦合迭代方法進(jìn)行數(shù)值計(jì)算，揭示高頻波對(duì)多孔彈性介質(zhì)中稠油溫度和黏度影響的主控因素及作用規(guī)律。

1 高頻波動(dòng)作用下多孔彈性介質(zhì)能量衰減與轉(zhuǎn)化動(dòng)力學(xué)模型

在無(wú)限大彈性介質(zhì)中前進(jìn)的波稱為行波，其在傳播過(guò)程中輸送能量和動(dòng)量，可用波動(dòng)方程進(jìn)行描述。假設(shè)在無(wú)限均勻彈性介質(zhì)中有一個(gè)無(wú)限大剛性物體沿法線方向振動(dòng)，這時(shí)所產(chǎn)生的是柱面波動(dòng)，波動(dòng)方程為

(1)

式中：ζ為某質(zhì)點(diǎn)某時(shí)刻位移，m；A為某質(zhì)點(diǎn)最大振幅，m；ω為角頻率，rad/s；τ為時(shí)間，s；r為質(zhì)點(diǎn)距振源的距離，m；C0為波動(dòng)傳播速度，m/s。

高頻波動(dòng)在彈性介質(zhì)中傳播時(shí)，會(huì)帶動(dòng)彈性介質(zhì)中的質(zhì)點(diǎn)振動(dòng)，導(dǎo)致彈性介質(zhì)中壓強(qiáng)出現(xiàn)交替變化。波動(dòng)場(chǎng)中某一點(diǎn)在某一瞬時(shí)所產(chǎn)生的壓強(qiáng)p1與沒(méi)有波場(chǎng)時(shí)同一點(diǎn)的靜態(tài)壓強(qiáng)p0之差稱為該點(diǎn)的聲壓[8]，用p表示，即p=p1-p0，

(2)

式中：p為聲壓，MPa；ρ0為振動(dòng)質(zhì)點(diǎn)密度，g/cm3。

設(shè)在波動(dòng)場(chǎng)中取一塊足夠小的單位體積，在未加載波動(dòng)前體積為V0，壓強(qiáng)為p0，密度為ρ0，在加載波動(dòng)后，由于波動(dòng)的擾動(dòng)使該單位體積得到動(dòng)能

(3)

式中：V0為體積，cm3；v為質(zhì)點(diǎn)振動(dòng)速度，m/s。

同時(shí)在該擾動(dòng)下，該單位體積壓強(qiáng)由p0升高到p0+p，于是該單位體積具有的位能

(4)

單位體積內(nèi)總的高頻波動(dòng)能量為動(dòng)能和位能之和，即

(5)

式中：ΔE為單位體積機(jī)械能，10-3J；ΔEk為單位體積動(dòng)能，10-3J；ΔEp為單位體積勢(shì)能，10-3J。

單位體積中高頻波動(dòng)能量稱為波動(dòng)能量密度，表示為

(6)

式中：ε為波動(dòng)能量密度，10-3J/cm3。

其中振動(dòng)速度

由式(6)得到高頻波動(dòng)下波動(dòng)能量密度[9]

(7)

波的衰減主要表現(xiàn)為其振幅的衰減，衰減形式表現(xiàn)為A=A0e-αr。

其中高頻波動(dòng)與低頻相比有較高的衰減系數(shù)，其衰減系數(shù)[10-11]

exp[-(0.001 622 8.79ρKfμ2)ρKpμ-1]。

(8)

式中：α為衰減系數(shù)； K為地層滲透率，μm2； μ為流體黏度，Pa·s； f為振動(dòng)頻率， kHz。

那么，波動(dòng)能量密度

(9)

能量密度在一個(gè)周期內(nèi)的平均值稱為平均能量密度，熱能的平均能量密度

(10)

熱能的平均能量密度即為導(dǎo)熱微分方程中的內(nèi)熱源，表示為

(11)

2 非均勻熱源徑向?qū)釀?dòng)力學(xué)模型

通過(guò)上一節(jié)的推導(dǎo)，可以得到高頻波動(dòng)在傳播過(guò)程中不斷衰減，其所具有的機(jī)械能向熱能轉(zhuǎn)化。為了計(jì)算該熱能的有效值，推導(dǎo)得到熱能的平均能量密度，即單位時(shí)間內(nèi)熱能的平均值，也就是導(dǎo)熱微分方程的熱源項(xiàng)，該熱源項(xiàng)中包含了波動(dòng)的振幅、頻率、衰減系數(shù)等相關(guān)參數(shù)。同時(shí)高頻波動(dòng)能量在徑向上的衰減隨振源距離的不同而不同，即導(dǎo)熱微分方程中內(nèi)熱源隨距離的不同而變化，所以利用穩(wěn)態(tài)非均勻熱源的圓筒壁導(dǎo)熱微分方程進(jìn)行推導(dǎo)，并作如下假設(shè)：

(1)振動(dòng)彈性介質(zhì)是無(wú)限大各向同性均勻連續(xù)的；

(2)彈性介質(zhì)的熱導(dǎo)率、密度等參數(shù)已知；

(3)彈性介質(zhì)內(nèi)熱源在徑向上隨振源距離不同而變化，在周向上各點(diǎn)內(nèi)熱源相等，不存在能量傳遞；

(4)高頻波動(dòng)僅對(duì)彈性介質(zhì)溫度和流體黏度有所影響，對(duì)其他參數(shù)的影響忽略。

采用非均勻熱源圓筒壁導(dǎo)熱微分方程[12]將熱能平均密度作為熱源項(xiàng)引入，設(shè)定邊界條件為無(wú)窮遠(yuǎn)處彈性介質(zhì)溫度不受影響且為原始溫度，同時(shí)無(wú)窮遠(yuǎn)處導(dǎo)熱熱能流量為零，即

(12)

式中：λ為導(dǎo)熱系數(shù)，kW/(K·m)；t為波動(dòng)作用下彈性介質(zhì)溫度，℃；t0為彈性介質(zhì)初始溫度，℃。

上述方程描述了在高頻波動(dòng)作用一段時(shí)間后，多孔彈性介質(zhì)中稠油溫度隨作用距離變化的穩(wěn)態(tài)方程?，F(xiàn)利用非均勻熱源圓筒壁非穩(wěn)態(tài)導(dǎo)熱微分方程[13-15]結(jié)合邊界條件和初始條件建立模型(假設(shè)條件相同)，描述高頻波動(dòng)作用到達(dá)到穩(wěn)態(tài)的過(guò)程，即

(13)

式中：τ為振動(dòng)時(shí)間，s。

3 稠油溫黏狀態(tài)方程

上一節(jié)中得到了高頻波動(dòng)作用下稠油溫度變化的偏微分方程及求解條件，從中可以看出溫度是波動(dòng)衰減系數(shù)的函數(shù)，而高頻波動(dòng)衰減系數(shù)又是黏度的函數(shù)，同時(shí)黏度對(duì)溫度比較敏感，因此，需要利用黏度預(yù)測(cè)模型對(duì)溫度變化后流體黏度進(jìn)行預(yù)測(cè)。油氣工業(yè)中常用的黏度預(yù)測(cè)模型有LBC黏度模型、CS黏度模型、LLS黏度模型、PT 黏度模型和PR黏度模型。經(jīng)過(guò)前人的研究發(fā)現(xiàn)PR黏度模型在預(yù)測(cè)黏度時(shí)總平均誤差最小。本文采用PR模型來(lái)預(yù)測(cè)波動(dòng)作用下儲(chǔ)層稠油黏度的變化，并進(jìn)行耦合迭代運(yùn)算。

已被廣泛應(yīng)用的PR模型的形式如下[16-18]：

(14)

Q1～Q3已經(jīng)普遍化為偏心因子的關(guān)聯(lián)式。

對(duì)于Ec<0.3，有

(15)

對(duì)于Ec≥0.3，有

(16)

式中：T為溫度，K；Ec為偏心因子；Z為壓縮因子；pc為臨界壓力，MPa；Tc為臨界溫度，K；pr為對(duì)比壓力，MPa；Tr為對(duì)比溫度，K；Mw為分子量，g/mol。

可用解析法求出μ的3個(gè)根，在對(duì)應(yīng)溫度、壓力下，流體黏度取值的判據(jù)如下：

(1)在汽、液兩相區(qū)，飽和汽、液相的黏度分別大于b′的最小實(shí)根和最大實(shí)根；

(2)在單一液相區(qū)和超臨界區(qū)(T>Tc)，流體黏度取方程的最大實(shí)根；

(3)在壓力低于流體所處溫度下的飽和蒸汽壓的氣相區(qū)時(shí)，流體黏度大于b′的最小實(shí)根。

應(yīng)用于混合物的黏度計(jì)算時(shí)，PR黏度方程可以表示為

(17)

4 稠油黏溫耦合模型的求解

基于上述所建模型，采用耦合迭代的方法進(jìn)行求解，對(duì)于距離振源ri處的某點(diǎn)，不同頻率和振幅下其黏度和溫度變化計(jì)算流程如下：

(1)查詢并計(jì)算模型求解所必須的基礎(chǔ)物性參數(shù)，如油藏溫度、壓力，稠油臨界溫度、壓力等；

(2)根據(jù)PR黏度預(yù)測(cè)模型，即式(14)計(jì)算未加載波動(dòng)作用下稠油在油藏條件下的黏度，式(14)中的必要參數(shù)在該小節(jié)中已列出，同時(shí)根據(jù)判別條件選取適合的解；

(3)根據(jù)衰減系數(shù)方程，即式(8)計(jì)算該黏度下高頻波動(dòng)的衰減系數(shù)；

(4)由所得到的衰減系數(shù)以及波動(dòng)的振幅、頻率等參數(shù)，利用非均勻熱源徑向?qū)釀?dòng)力學(xué)方程，即式(12)、式(13)可以分別計(jì)算ri處的稠油溫度；

(5)利用PR黏度預(yù)測(cè)模型，即式(14)預(yù)測(cè)該溫度下稠油的黏度；

(6)將步驟(5)中得到的稠油黏度與之前的稠油黏度比較并計(jì)算誤差，若在誤差允許范圍內(nèi)則停止計(jì)算，若不滿足誤差要求則重復(fù)步驟(3)—(6)；

(7)由此得到高頻波動(dòng)作用下距離振源ri處稠油黏度及溫度的變化，繼續(xù)向油藏深處推進(jìn)計(jì)算距離振源ri+1處稠油黏度及溫度的變化。

5 結(jié)果與討論

基于上述迭代計(jì)算流程給出計(jì)算所需要的基本數(shù)據(jù)[19]如表1所示，在基本數(shù)據(jù)的基礎(chǔ)上利用MATLAB編程運(yùn)算并繪制圖形，討論高頻波動(dòng)在多孔彈性介質(zhì)中對(duì)稠油黏溫參數(shù)的影響規(guī)律。

5.1 波動(dòng)頻率對(duì)稠油黏溫參數(shù)影響規(guī)律

不同波動(dòng)頻率(10 kHz、20 kHz、30 kHz)作用下稠油溫度及黏度隨作用距離的變化見(jiàn)圖1，這里保持振幅為1×10-5m。

表1 波動(dòng)作用下稠油黏溫模擬基本數(shù)據(jù)表Tab.1 Basic data for viscosity-temperature simulation of heavy oil under vibration

由圖1可以看出，高頻波動(dòng)對(duì)稠油溫度和黏度的影響近似呈指數(shù)變化，距振源由近至遠(yuǎn)溫度呈指數(shù)衰減，流體黏度呈指數(shù)升高，有效作用范圍在0.6 m左右，與實(shí)驗(yàn)中超聲波處理油層半徑范圍近似，這是由于波動(dòng)較高的頻率導(dǎo)致其衰減系數(shù)較大，波動(dòng)能量在近井地帶就得到較大損耗，無(wú)法作用于較遠(yuǎn)距離的油層。

圖1 不同波動(dòng)頻率下稠油溫度及黏度隨作用距離變化曲線Fig.1 Variation of temperature and viscosity of heavy oil with radial distance under different vibration frequencies

同時(shí)由圖1可以看出，隨著波動(dòng)頻率的提高，近井地帶溫度升高，流體黏度降低，但是溫度的增加量和黏度的減小量在降低。同時(shí)隨著頻率的提高，遠(yuǎn)井地帶在頻率較高的波動(dòng)作用下，油層溫度較低，稠油黏度較高。20 kHz曲線與10 kHz曲線的交匯點(diǎn)(即2個(gè)頻率的等效作用點(diǎn))在0.18 m左右，30 kHz曲線與20 kHz曲線的交匯點(diǎn)在0.10 m左右，這是由于波動(dòng)頻率越高，衰減系數(shù)越大，其能量集中損耗在近井地帶，對(duì)于遠(yuǎn)井地帶的作用反而不如頻率較低的波動(dòng)，同時(shí)也印證了高頻波動(dòng)有較好的近井處理效果，而低頻波動(dòng)有較大處理范圍的結(jié)論。

5.2 波動(dòng)振幅對(duì)稠油黏溫參數(shù)影響規(guī)律

不同振幅(0.5×10-5m、1×10-5m、1.5×10-5m)的波動(dòng)作用下稠油溫度和黏度隨作用距離的變化見(jiàn)圖2，這里保持波動(dòng)頻率為20 kHz。

由圖2可以看出，不同振幅下，波動(dòng)對(duì)稠油溫度及黏度的影響依然表現(xiàn)為指數(shù)形式，同時(shí)相對(duì)于波動(dòng)頻率，振幅對(duì)波動(dòng)作用效果有比較大的影響(在頻率和振幅提高幅度相同的條件下比較)，隨著振幅的提高，近井地帶多孔彈性介質(zhì)溫度加速升高而流體黏度加速降低，這是由于波動(dòng)的能量與振幅的平方成正比，振幅的改變會(huì)導(dǎo)致波動(dòng)能量的大幅變化。另一個(gè)方面，隨著振幅的提高，高頻波動(dòng)的處理范圍增大，振幅為5 μm時(shí)其作用范圍約為0.4 m，15 μm時(shí)其作用范圍約為0.7 m。原因是在衰減因子相同的情況下，振幅提高，波動(dòng)能量提高，可以傳播更遠(yuǎn)的距離。通過(guò)計(jì)算可以得到，3種振幅由高到低，平均降黏率分別為31.06%、16.96%、6.46%，平均溫度升高率分別為58.45%、29.94%、9.69%。

圖2 不同振幅下稠油溫度及黏度隨作用距離變化曲線Fig.2 Variation of temperature and viscosity of heavy oil with radial distance under different vibration amplitudes

5.3 高頻波動(dòng)下稠油溫度和黏度隨作用時(shí)間、作用距離的變化規(guī)律

在振幅為10 μm、頻率為20 kHz的高頻波動(dòng)作用下，稠油溫度和黏度隨作用時(shí)間和作用距離的變化如圖3所示。

由圖3可以看出，近井地帶隨著作用時(shí)間的增長(zhǎng)，稠油溫度不斷升高，但是升高的速度逐漸變緩，相對(duì)應(yīng)的，稠油黏度不斷降低，但降低的速度逐漸變緩。分析原因是由于隨著高頻波動(dòng)的作用，近井地區(qū)的溫度升高，單位距離的溫差增大，導(dǎo)致徑向熱量傳遞加快，最終在某一溫差下波動(dòng)提供的熱量與徑向散失的熱量平衡，近井地帶稠油的溫度和黏度便趨向于平穩(wěn)。由高頻波動(dòng)開(kāi)始作用至稠油溫度和黏度趨于穩(wěn)定，所用時(shí)間大約是400 s，最終在穩(wěn)態(tài)下稠油溫度和黏度變化規(guī)律(作用效果、變化趨勢(shì)、有效作用距離等)與5.1中的計(jì)算結(jié)果相符。

圖3 高頻波動(dòng)下稠油溫度和黏度隨作用時(shí)間、作用距離變化關(guān)系Fig.3 Variation of temperature and viscosity of heavy oil with time and radial distance under high frequency vibration

6 結(jié) 論

(1)高頻波在多孔彈性介質(zhì)中的指數(shù)式衰減導(dǎo)致近井地帶溫度及流體黏度也呈指數(shù)形式變化，有效作用距離在0.6 m左右，與以往的實(shí)驗(yàn)結(jié)果近似。同時(shí)30 kHz與20 kHz波動(dòng)的等效作用距離為0.1 m，而20 kHz與10 kHz為0.18 m，印證了高頻波動(dòng)有較好的近井處理效果，而低頻波動(dòng)有較大的處理范圍。

(2)相對(duì)于頻率，振幅對(duì)近井地帶溫度和黏度的影響更加明顯。在波動(dòng)頻率為20 kHz、振幅15 μm下，近井地帶平均降黏率為31.06%，平均溫度升高率為58.45%。由高頻波動(dòng)開(kāi)始作用至稠油溫度和黏度趨于穩(wěn)定，所用時(shí)間大約是400 s。

(3)根據(jù)本文結(jié)果，波動(dòng)近井解堵不宜采用過(guò)高的振動(dòng)頻率，否則波場(chǎng)在近井地區(qū)衰減過(guò)快，無(wú)法實(shí)現(xiàn)良好的解堵效果。同時(shí)提高振動(dòng)功率也就是振幅可以明顯改善作用距離和效果。若要處理油藏內(nèi)部宜采用作用距離更遠(yuǎn)的低頻振動(dòng)。

[1] 胡博仲.波場(chǎng)采油[M].北京:石油工業(yè)出版社,1996:31-36.

[2] GUO Xiao,DU Zhimin,FU Yu,et al.Effect of formation damage to production performance in heavy oil reservoir via steam-injection after water-flooding[C].SPE 97869,2005.

[3] 張人雄,向陽(yáng),李曉梅.超稠油油藏?zé)岵蛇^(guò)程中瀝青質(zhì)損害室內(nèi)模擬[J].石油勘探與開(kāi)發(fā),2001,28(1):85-86. ZHANG Renxiong,XIANG Yang,LI Xiaomei.Asphaltene damage simulation in super heavy oil reservoir during thermal recovery process[J].Petroleum Exploration and Development,2001,28(1):85-86.

[4] 朱靜,李傳憲,楊飛,等.稠油降黏新技術(shù)的研究進(jìn)展[J].西安石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,27(1):64-70. ZHU Jing,LI Chuanxian,YANG Fei,et al.Research progress in new viscosity-reducing techniques for heavy oil[J].Journal of Xi’an Shiyou University(Natural Science Edition),2012,27(1):64-70.

[5] 谷瀟雨,蒲春生,王蓓,等.超聲波解除巖心鉆井液堵塞實(shí)驗(yàn)研究[J].西安石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,29(1):76-79. GU Xiaoyu,PU Chunsheng,WANG Bei,et al.Experimental study on removing drilling fluid damage of core by ultrasonic[J].Journal of Xi’an Shiyou University(Natural Science Edition),2014,29(1):76-79.

[6] 蒲春生,石道涵,趙樹(shù)山,等.大功率超聲波近井處理無(wú)機(jī)垢堵塞技術(shù)[J].石油勘探與開(kāi)發(fā),2011,38(2):243-248. PU Chunsheng,SHI Daohan,ZHAO Shushan,et al.Technology of removing near wellbore inorganic scale damage by high power ultrasonic treatment[J].Petroleum Exploration and Development,2011,38(2):243-248.

[7] 傅俊萍,李錄平,劉澤利,等.超聲波除垢與強(qiáng)化傳熱實(shí)驗(yàn)研究[J].熱能動(dòng)力工程,2006,21(4):355-357. FU Junping,LI Luping,LIU Zeli,et al.Experimental study on ultrasonic scale removal and heat transfer enhancement[J].Journal of Engineering for Thermal Energy and Power,2006,21(4):355-357.

[8] 王仲茂.振動(dòng)采油[M].北京:石油工業(yè)出版社,2000.

[9] 杜功煥,朱哲民,龔秀芬.聲學(xué)基礎(chǔ)[M].南京：南京大學(xué)出版社,2001:190-192.

[10] BIOT M A.Mechanics of deformation and acoustic propagation in porous media[J].Journal of Applied Physics,1962,33(4):1482-1498.

[11] 張全興.超聲波非均勻介質(zhì)傳播衰減特性研究[D].沈陽(yáng):沈陽(yáng)工業(yè)大學(xué),2015.

[12] 黃曉齊,龔德鴻,李環(huán).非均勻內(nèi)熱源圓筒壁穩(wěn)態(tài)導(dǎo)熱微分方程的熱流分流和疊加[J].貴州大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,30(3):9-12. HUANG Xiaoqi,GONG Dehong,LI Huan.Heat flow distributing and stacking of non-uniform heat source cylinder steady heat-conduction differential equations[J].Journal of Guizhou University(Natural Science Edition),2001,30(3):9-12.

[13] 戴鍋生.傳熱學(xué)[M].2版.北京:高等教育出版社,1999.

[14] 王永巖,秦楠,蘇傳奇.無(wú)限大平板非穩(wěn)態(tài)導(dǎo)熱過(guò)程的數(shù)字特征[J].青島科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,34(5):511-515. WANG Yongyan,QIN Nan,SU Chuanqi.Research on numerical characteristics of unsteady heat conduction of an infinite plat[J].Journal of Qingdao University of Science and Technology(Natural Science Edition),2013,34(5):511-515.[15] 顧祥紅.無(wú)限大平板非穩(wěn)態(tài)導(dǎo)熱理論分析[J].沈陽(yáng)工業(yè)大學(xué)學(xué)報(bào),2009,31(1):112-115. GU Xianghong.Theoretical analysis of unsteady-state thermal conduction of infinite plate[J].Journal of Shenyang University of Technology,2009,31(1):112-115.

[16] 郭緒強(qiáng),榮淑霞,楊繼濤,等.基于PR狀態(tài)方程的黏度模型[J].石油學(xué)報(bào),1999,20(3):56-61. GUO Xuqiang,RONG Shuxia,YANG Jitao,et al.The viscosity model based on PR state equation[J].Acta Petrolei Sinica,1999,20(3):56-61.

[17] 雷宇.稠油黏度預(yù)測(cè)模型研究[D].成都:西南石油大學(xué),2012.

[18] 馮國(guó)強(qiáng),陳爽,郭緒強(qiáng),等.油藏流體黏度預(yù)測(cè)模型的對(duì)比分析[J].油氣地質(zhì)與采收率,2007,14(5):80-83. FENG Guoqiang,CHEN Shuang,GUO Xuqiang,et al.Comparative analysis of reservoir fluid viscosity prediction model[J].Petroleum Geology and Recovery Efficiency,2007,14(5):80-83.

[19] 王欣,趙法軍,劉江,等.稠油比熱容導(dǎo)熱系數(shù)及黏溫物性參數(shù)變化研究[J].化工科技,2012,20(3):10-13. WANG Xin,ZHAO Fajun,LIU Jiang,et al.Research on thermal transmissivity heat coefficient and temperature-viscosity parameters changes of heavy oil[J].Science & Technology in Chemical Industry,2012,20(3):10-13.

責(zé)任編輯：賀元旦

Study on Influence of High Frequency Vibration on Viscosity and Temperature Parameters of Heavy Oil in Porous Elastic Medium

PU Chunsheng1,2,XU Jiaxiang1,2,LIU Xi3,ZHENG Liming1,2,LIU Jing1,2

(1.Faculty of Petroleum Engineering,China University of Petroleum (East China),Qingdao 266580,Shandong,China;2.State Key Laboratory of Heavy Oil Processing,China University of Petroleum (East China),Qingdao 266580,Shandong,China;3.Western Oil Production Plant,Yanchang Oilfield Company,Yan'an 717500,Shaanxi,China)

In order to quantitatively describe the thermodynamic viscosity-reducing and plug-removing mechanism of high-frequency vibration excitation near wellbore,the steady and unsteady coupling dynamic models of heavy oil temperature and viscosity variation near wellbore of heavy oil reservoir caused by high-frequency vibration excitation are established respectively based on the energy transmission,loss and transformation kinetics of high-frequency vibration wave in poroelastic medium.The models consist of energy attenuation and transformation equation,the differential equation of inhomogeneous heat source conduction in cylindrical wall and the temperature-viscosity equation of heavy oil,and they are solved using coupling iterative method to reveal the main controlling factors and laws of high-frequency vibration influencing the temperature and viscosity variation of heavy oil in poroelastic medium.The results show that the vibration of 20 kHz has a better effect within 0.6m;the higher the frequency of vibration,the better the viscosity-reducing and plug-removing effect of vibration,but the smaller the effect range of vibration;compared with vibration frequency,vibration amplitude has greater effect;it takes about 400 s for the effect of vibration to transit from non-steady state to steady state.

recovery of heavy oil;high-frequency vibration;viscosity-reducing and plug-removing

2016-08-30

國(guó)家“863”計(jì)劃重大專項(xiàng)課題“超大功率超聲波油井增油技術(shù)及其裝置研究”(編號(hào):2007AA06Z227)；山東省自然科學(xué)基金“稠油熱波耦合輔助層內(nèi)化學(xué)催化裂解降黏技術(shù)研究”(編號(hào):ZR2010EM014)；中國(guó)石油化工股份有限公司重大科技專項(xiàng)項(xiàng)目“稠油注蒸汽熱波耦合催化裂解降黏技術(shù)研究”(編號(hào):[2011]075)

蒲春生(1959-)，男，教授、博士生導(dǎo)師，主要從事特種油氣藏改造與物理-化學(xué)增產(chǎn)增注理論與技術(shù)研究。 E-mail:chshpu@163.com

10.3969/j.issn.1673-064X.2016.06.008

TE345

1673-064X(2016)06-0054-06

A