亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        用非結(jié)構(gòu)化網(wǎng)格模擬器模擬聚合物驅(qū)

        2010-10-13 08:00:18編譯周長江西南石油大學(xué)研究生部
        石油石化節(jié)能 2010年7期
        關(guān)鍵詞:水相摩爾結(jié)構(gòu)化

        編譯:周長江 (西南石油大學(xué)研究生部)

        審校:施雷庭 (油氣藏地質(zhì)與開發(fā)工程國家重點(diǎn)實(shí)驗(yàn)室·西南石油大學(xué))

        用非結(jié)構(gòu)化網(wǎng)格模擬器模擬聚合物驅(qū)

        編譯:周長江 (西南石油大學(xué)研究生部)

        審校:施雷庭 (油氣藏地質(zhì)與開發(fā)工程國家重點(diǎn)實(shí)驗(yàn)室·西南石油大學(xué))

        描述了用非結(jié)構(gòu)化網(wǎng)格模擬聚合物驅(qū)的研究進(jìn)展。聚合物驅(qū)模擬涉及到很多專業(yè)問題:聚合物溶液黏度 (具有剪切稀釋或剪切增稠的特性)是聚合物濃度、剪切速率、礦化度和油藏溫度的函數(shù);聚合物的流度控制能力受近井地帶高速剪切的影響;聚合物的熱力學(xué)降解是聚合物濃度、礦化度和油藏溫度的函數(shù);聚合物的滯留會影響其流動特征。除此之外,文中還提供了方程的數(shù)值解法及其應(yīng)用實(shí)例。

        聚合物驅(qū) 模擬 非結(jié)構(gòu)化網(wǎng)格

        1 引言

        聚合物驅(qū)已經(jīng)在油田大規(guī)模應(yīng)用,其概念很簡單,即在注入水中加入少量水溶性聚合物以增加注入水的黏度,實(shí)現(xiàn)有利的流度比,從而提高水驅(qū)效果;但是,其機(jī)理是很復(fù)雜的。聚合物溶液黏度受聚合物濃度、礦化度、溫度及剪切速率等因素的影響;吸附和滯留降低了聚合物濃度和水相滲透率;降解作用使聚合物喪失其增黏能力。因此,合理的聚合物驅(qū)設(shè)計和評價必須結(jié)合油藏特征、室內(nèi)實(shí)驗(yàn)、油藏模擬、工藝技術(shù)及礦場試驗(yàn)。其中,準(zhǔn)確的油藏模擬是聚合物驅(qū)設(shè)計和效果評價的關(guān)鍵。

        從20世紀(jì)60年代后期開始,很多學(xué)者進(jìn)行了聚合物驅(qū)數(shù)值模擬研究。1968年,Zeito提出了一個簡單的聚合物驅(qū)數(shù)學(xué)模型:水相黏度是聚合物濃度的函數(shù)。1969年,Jewett和Shurz根據(jù)近60個聚合物驅(qū)項(xiàng)目,開發(fā)了一個兩相、多層、能夠模擬一維介質(zhì)和五點(diǎn)井網(wǎng)的模擬器。1970年,Slater和Farouq Ali采用阻力系數(shù)經(jīng)驗(yàn)公式和差分方法模擬五點(diǎn)井網(wǎng)聚合物驅(qū)。Bondor和 Hirasaki(1972年)以及Hirasaki和Pope(1974年)開發(fā)的數(shù)學(xué)模型認(rèn)為聚合物只溶于水中,他們用修正的Blake -Kozeny模型并考慮了聚合物吸附、滲透率下降系數(shù)和聚合物溶液非牛頓流變特性的影響。1976年,Vela、Peaceman和 Sandvik開發(fā)的模擬器考慮了聚合物滯留、不可入孔隙體積、剪切降解和礦化度的影響。許多學(xué)者還考慮了化學(xué)降解和熱力學(xué)降解對聚合物溶液特性的影響。

        20世紀(jì)80年代以來,油藏模擬變得越來越復(fù)雜,對計算的準(zhǔn)確性要求也越來越高。通過對工藝技術(shù)和流線的研究,現(xiàn)代油藏模擬力求實(shí)現(xiàn)對細(xì)小裂縫、復(fù)雜井眼軌跡等多方面的影響因素進(jìn)行模擬。與常規(guī)網(wǎng)格系統(tǒng)相比,非結(jié)構(gòu)化網(wǎng)格系統(tǒng)能夠獲得更多、更準(zhǔn)確的地質(zhì)和工程信息。對常規(guī)網(wǎng)格系統(tǒng)的研究比較成熟,而關(guān)于非結(jié)構(gòu)化網(wǎng)格系統(tǒng)的研究卻較少。2001年,Beckner等人為ExxonMobil公司開發(fā)了非結(jié)構(gòu)化網(wǎng)格模擬器,并對一個擁有復(fù)雜地質(zhì)裂縫和多個油藏、采用常規(guī)工藝技術(shù)生產(chǎn)的油田進(jìn)行了模擬。研究結(jié)果表明:大約10萬個非結(jié)構(gòu)化網(wǎng)格就能實(shí)現(xiàn)相當(dāng)于160萬個正交網(wǎng)格才能實(shí)現(xiàn)的地質(zhì)裂縫的準(zhǔn)確模擬,所以該模擬器簡單易用,大大縮減了模擬時間。

        常規(guī)油藏模擬的基本方程見附錄A(附錄見原文,下同)。這些方程中包括體積守恒方程和組分摩爾守恒方程??傮w積守恒方程A-8用來求解壓力,各相體積守恒方程A-7用來求解相飽和度,摩爾守恒方程A-18用來求解各組分摩爾數(shù)。

        2 聚合物驅(qū)的物理模型

        正如引言中提到的那樣,聚合物驅(qū)機(jī)理是非常復(fù)雜的。EMpower是 ExxonMobil公司在MARS模擬器基礎(chǔ)上開發(fā)的新型模擬器。其特點(diǎn)包括:用非結(jié)構(gòu)化網(wǎng)格模擬聚合物驅(qū);聚合物溶液性質(zhì)與礦化度、剪切增稠和聚合物降解等因素有關(guān)。該模擬器結(jié)合油藏中流體的流動、布井方式、流線和工藝技術(shù),能夠準(zhǔn)確有效地模擬復(fù)雜油藏。

        2.1 聚合物溶液黏度

        聚合物溶液是非牛頓流體,在低剪切速率下,黏度隨剪切速率的增加而降低。但是,一些學(xué)者已經(jīng)證明了聚合物溶液在多孔介質(zhì)中具有明顯的剪切增稠特性。2008年,Delshad等人證明了剪切增稠是德布拉數(shù)的函數(shù)。EMpower既能模擬聚合物的剪切稀釋,又能模擬剪切增稠。除剪切速率外,聚合物溶液黏度還是聚合物濃度、礦化度和溫度的函數(shù)。

        2.2 傳導(dǎo)率修正

        通常,近井地帶的剪切速率高于油藏中的剪切速率。因此,從井底到油藏中剪切速率的變化需要特殊處理,尤其是當(dāng)井眼位于大網(wǎng)格塊中時。在本文的研究過程中,根據(jù)近井地帶剪切速率的變化對傳導(dǎo)率進(jìn)行修正,這可以通過模擬井眼周圍壓降對黏度的影響來實(shí)現(xiàn)。

        2.3 不可入孔隙體積

        不可入孔隙體積 (IPV)描述了聚合物不能進(jìn)入的孔隙空間。IPV可理解為:孔壁排斥效應(yīng)和孔喉排斥效應(yīng)。在孔壁排斥效應(yīng)中,聚合物溶液流經(jīng)所有孔隙,但聚合物集中在高速區(qū)而遠(yuǎn)離孔隙表面。在孔喉排斥效應(yīng)中,聚合物只流經(jīng)大孔隙 (含一部分中孔隙)而不流經(jīng)小孔隙。在 EMpower中,考慮到IPV的影響,對聚合物的流速進(jìn)行了調(diào)整。IPV與巖石類型和滲透率有關(guān)。

        2.4 聚合物吸附

        吸附和滯留是聚合物在多孔介質(zhì)中損失的主要原因,與滲透率 (代替孔隙表面積)有關(guān)。EMpower可將聚合物吸附分為可逆和不可逆兩種情況進(jìn)行模擬。

        2.5 殘余阻力系數(shù) (RRF)

        聚合物吸附可引起孔隙縮小,從而降低水相滲透率。水相滲透率下降是聚合物吸附量的函數(shù)。通常,殘余阻力系數(shù)只降低水相滲透率,幾乎不影響油相滲透率。

        2.6 聚合物降解

        聚合物不僅會發(fā)生熱力學(xué)降解,還會與地層水中的金屬陽離子反應(yīng)而沉降。不同的油田,對聚合物耐鹽性有不同的要求。本文在模擬聚合物降解時考慮了溫度和各種金屬陽離子的影響。

        3 數(shù)學(xué)方程

        模擬時間步的目的是為了獲得在該時間步結(jié)束時的壓力、組分?jǐn)?shù)量及相飽和度,可根據(jù)巖石和流體的性質(zhì)獲得這些參數(shù)的初始值。文中建立了多個數(shù)學(xué)方程,通過IMPES、耦合-隱含、順序及IMPSAT算法來求解這些未知數(shù)。在解出壓力、飽和度及其他組分的摩爾數(shù)之后,再求解聚合物的摩爾數(shù)。方程的選擇需要綜合考慮穩(wěn)定性和計算機(jī)的運(yùn)算能力。

        根據(jù)體積守恒方程A-8可推導(dǎo)出求解壓力變化的線性方程。線性方程A-15的算法決定了求解組分摩爾流量的復(fù)雜程度。

        3.1 IMPES算法

        對于IMPES和順序算法,壓力方程如式 (1)所示:

        由于其穩(wěn)定性差使用得較少。在IMPES算法中,壓力隨時間變化,而飽和度是固定不變的。當(dāng)流速較高、前緣飽和度變化較快的時候,不能很好地模擬油藏實(shí)際情況。

        3.2 耦合-隱含算法

        1980年,Coats提出了耦合-隱含 (CI)算法,可同時準(zhǔn)確求解壓力和組分摩爾數(shù)。與IMPES算法相比,它更穩(wěn)定,但運(yùn)算速度較慢。CI算法的組分流動方程如下所示:

        組分的摩爾守恒方程見A-15。對式 (2)相飽和度求導(dǎo) (因?yàn)榉匠讨酗柡投仁且蜃兞?便可得相體積守恒方程A-7。體積守恒方程A-8用來求解壓力。

        3.3 順序算法

        1973年,Spillette、Hillestad和Stone提出了一種介于IMPES算法和CI算法之間的順序算法。1986年,Watts用同樣的方法建立了組分模型。與CI算法同時求解壓力和摩爾數(shù)不同,順序算法是先求解壓力,然后依次求解飽和度和摩爾數(shù)。首先,用順序算法將壓力方程線性化以求解壓力變化;然后將飽和度方程線性化,求出飽和度變化;再將組分的摩爾守恒方程A-15線性化,便可計算組分?jǐn)?shù)變化;最后,根據(jù)最新的組成計算下一時間步的飽和度。

        3.4 IMPSAT算法

        1989年,Quandalle和Savary首次提出了IMPSAT算法,稱之為隱式壓力、隱式飽和度算法。1999年,Cao和Aziz在傳導(dǎo)方程中顯式計算相密度和組分的摩爾分?jǐn)?shù),研究表明,壓力與組分摩爾分?jǐn)?shù)之間的相關(guān)度要比壓力與相飽和度之間的相關(guān)度弱。根據(jù)定義,IMPSAT算法同時隱式求解壓力和飽和度,而顯式求解組分摩爾分?jǐn)?shù),因此,它相似于但又不同于分兩步求解壓力和飽和度的順序算法。在這幾種算法中,隱式程度由低到高的順序?yàn)镮MPES、順序、IMPSA T、CI。與之對應(yīng),IMPSAT算法在運(yùn)算上的穩(wěn)定性比IMPES算法和順序算法要好,但比CI算法差。因?yàn)镮MPSA T算法每個節(jié)點(diǎn)可以求解3個未知數(shù),所以比IMPES算法 (每個節(jié)點(diǎn)求解1個未知數(shù))和順序算法 (每個節(jié)點(diǎn)求解1+2個未知數(shù))的運(yùn)算成本高,但與全隱式算法相比,對計算機(jī)內(nèi)存和CPU的要求相對較低。

        4 聚合物摩爾守恒

        為模擬聚合物驅(qū)油過程,還需增加聚合物和鹽的摩爾守恒方程。式 (3)所列方程包含水相中溶解的聚合物、巖石表面吸附的聚合物、降解損失以及水相中作為示蹤劑的聚合物的量。該方程包括黏度和吸附對水相流度的影響以及IPV對聚合物流速的影響。

        式中 m——聚合物組分;

        a——水相;

        r——巖石;

        d——降解后的聚合物;

        κ——熱降解常數(shù);

        Nm,a——溶解聚合物的摩爾數(shù);

        Nm,r——巖石表面吸附的聚合物摩爾數(shù);

        Nm,d——降解的聚合物摩爾數(shù);

        F——模擬低滲透巖石對聚合物的“封堵”系數(shù);

        f——不可入孔隙體積的函數(shù)。上式可簡化為:

        式中,t表示聚合物總摩爾數(shù)。在每個網(wǎng)格中,將聚合物吸附和降解看作點(diǎn)源和點(diǎn)匯。顯式算法中,在時間步開始時求解這些參數(shù)并假設(shè)其在時間步長中保持為常數(shù)。降解常數(shù)κ由用戶自定義。先計算降解再計算吸附。吸附是水相中聚合物濃度的函數(shù),并符合Langmuir等溫吸附。殘余阻力系數(shù)是聚合物在巖石表面吸附量的函數(shù),反映了該網(wǎng)格聚合物的吸附量和降低水相滲透率的能力。獲得吸附降解后的聚合物濃度,便可計算水相黏度,從而通過方程 A-4可計算相流度。通常,f=1.0-IPV (IPV由用戶自定義),并存在區(qū)域性變化。F是模擬低滲透巖石對聚合物的“封堵”系數(shù),是網(wǎng)格滲透率的函數(shù)。雖然聚合物摩爾數(shù)可通過顯式或隱式求解,但還是應(yīng)根據(jù)穩(wěn)定性來決定采用哪種方式求解。假設(shè)水相密度不受聚合物濃度的影響。

        5 模擬結(jié)果

        5.1 模擬巖心驅(qū)替實(shí)驗(yàn)

        將EMpower模擬聚合物驅(qū)的采收率情況與Masuda等人1992年的巖心驅(qū)替實(shí)驗(yàn)結(jié)果進(jìn)行比較。Masuda等人測量了 HPAM溶液在玻璃微珠人造多孔介質(zhì)模型中的流變性 (圖1),并通過巖心驅(qū)替實(shí)驗(yàn)測量了聚合物驅(qū)采收率提高幅度。在不調(diào)整模型參數(shù)的情況下,EMpower模擬結(jié)果與實(shí)驗(yàn)結(jié)果的擬合度較高。實(shí)驗(yàn)和模擬結(jié)果如圖2和圖3所示。

        圖1 HPAM溶液 (200 ppm)的流變曲線

        圖2 采收率隨注入PV數(shù)的變化曲線

        5.2 非結(jié)構(gòu)化網(wǎng)格模型

        2001年,Beckner等人為 ExxonMobil公司開發(fā)了非結(jié)構(gòu)化網(wǎng)格模擬器 EMpower。對復(fù)雜結(jié)構(gòu)井(如大斜度井、水平井和多分支井)、斷層和尖滅進(jìn)行網(wǎng)格劃分是比較復(fù)雜的,但使用非結(jié)構(gòu)化網(wǎng)格就比較簡單。

        圖3 含油率隨注入PV數(shù)的變化曲線

        圖4 五點(diǎn)井網(wǎng)非結(jié)構(gòu)化網(wǎng)格示意圖

        圖5 巖石和流體性質(zhì)

        非結(jié)構(gòu)化網(wǎng)格允許使用多種網(wǎng)格技術(shù)對五點(diǎn)井網(wǎng)模型進(jìn)行模擬 (圖4)。網(wǎng)格結(jié)構(gòu)對聚合物驅(qū)效果影響很大,而非結(jié)構(gòu)化網(wǎng)格就降低了這樣的影響。當(dāng)聚合物溶液黏度低于油相黏度 (如聚合物驅(qū)替稠油)或者高黏聚合物驅(qū)轉(zhuǎn)后續(xù)水驅(qū)時,容易發(fā)生黏性指進(jìn),應(yīng)特別注意網(wǎng)格的方位敏感性。五點(diǎn)井網(wǎng)模型聚合物驅(qū)實(shí)例的網(wǎng)格結(jié)構(gòu)如圖4(a)所示。模型尺寸為500 ft×50 ft×50 ft(1 ft=30.48 cm)。聚合物溶液的流變性如圖5所示。聚合物遵循Langmuir等溫吸附,最大吸附量為20μg/g砂。IPV為0.2。水相滲透率下降了20%,但油相滲透率不受影響。模型孔隙度為26%,滲透率為2 D。油水相對滲透率曲線如圖5所示。油藏初始含水飽和度為束縛水飽和度,初始壓力為2 140 psia(1 psi=6.895 kPa)。為了簡化,不考慮礦化度、聚合物降解和油水壓縮性的影響 (體積系數(shù)為1),并且只模擬二維的情況。聚合物注入濃度為1 500 ppm(1 ppm=10-6),注聚時間為10年。產(chǎn)油井井底壓力為1 500 psia。五點(diǎn)井網(wǎng)模型復(fù)雜非結(jié)構(gòu)化網(wǎng)格中含油飽和度和聚合物濃度如圖6所示。目前已開展針對實(shí)際地質(zhì)模型的模擬,運(yùn)算速度較快。

        圖6 注聚5 730天后的含油飽和度和聚合物濃度分布

        6 結(jié)論

        非結(jié)構(gòu)化網(wǎng)格可以詳細(xì)地模擬聚合物驅(qū)過程,并有效地模擬復(fù)雜油藏。

        聚合物驅(qū)模擬包括聚合物的流動特性和驅(qū)替效果,其驅(qū)替特征如下:

        ◇聚合物溶液黏度是聚合物濃度、剪切速率、礦化度和油藏溫度的函數(shù);

        ◇聚合物溶液黏度具有剪切稀釋或剪切增稠的特性;

        ◇傳導(dǎo)系數(shù)受近井地帶高速剪切的影響;

        ◇聚合物降解是聚合物濃度、離子濃度和油藏溫度的函數(shù);

        ◇聚合物滯留對聚合物流動特性有重要影響;

        ◇可使用多種數(shù)值算法 (如順序算法,IMPSAT等)和多種模型 (如黑油模型、組分模型和熱采模型)進(jìn)行聚合物驅(qū)模擬。

        資料來源于美國《SPE 118985》

        10.3969/j.issn.1002-641X.2010.7.001

        2009-04-08)

        猜你喜歡
        水相摩爾結(jié)構(gòu)化
        促進(jìn)知識結(jié)構(gòu)化的主題式復(fù)習(xí)初探
        戰(zhàn)場上的雕塑家——亨利摩爾
        河北畫報(2020年10期)2020-11-26 07:20:56
        結(jié)構(gòu)化面試方法在研究生復(fù)試中的應(yīng)用
        海上中高滲透率砂巖油藏油水相滲曲線合理性綜合分析技術(shù)
        更 正
        西方摩爾研究概觀
        地下水流速與介質(zhì)非均質(zhì)性對于重非水相流體運(yùn)移的影響
        基于圖模型的通用半結(jié)構(gòu)化數(shù)據(jù)檢索
        用三辛胺和磷酸三丁酯萃取、銨溶液反萃取鉬的研究
        濕法冶金(2014年3期)2014-04-08 01:04:51
        基于軟信息的結(jié)構(gòu)化轉(zhuǎn)換
        中文字幕有码人妻在线| 亚洲AV无码成人精品区天堂| 国产主播一区二区在线观看| 新久久国产色av免费看| 日日摸天天摸97狠狠婷婷| 国产人澡人澡澡澡人碰视频| 中国女人a毛片免费全部播放| 国产精品伦理久久一区| 久久亚洲av成人无码电影a片| 激情偷乱人伦小说视频在线| 色优网久久国产精品| 国产主播性色av福利精品一区| 久久午夜无码鲁丝片午夜精品| 日本动态120秒免费| 亚洲色图在线视频免费观看 | 亚洲视频免费一区二区| 亚洲av日韩av天堂一区二区三区| 手机看片福利盒子久久青| 狼人综合干伊人网在线观看| av一区二区三区在线| 免费国产裸体美女视频全黄| 日韩久久久黄色一级av| 亚洲av调教捆绑一区二区三区| 黄色三级视频中文字幕| 日本在线一区二区三区视频观看 | 精品国产一区二区三区av麻| 巨茎中出肉欲人妻在线视频| 亚洲免费不卡| 无码AV午夜福利一区| 亚洲不卡av一区二区三区四区| 中文字幕日本人妻久久久免费 | 97久久综合区小说区图片专区| 午夜少妇高潮在线观看视频| 国产免费艾彩sm调教视频| 国产欧美乱夫不卡无乱码| 伊人久久大香线蕉综合av| 日本孕妇潮喷高潮视频| 熟妇人妻无码中文字幕| 北岛玲精品一区二区三区| 五月天中文字幕日韩在线| 国产精品无码专区av在线播放|