一維徑向流數(shù)值模擬行業(yè)資料系統(tǒng)集成

e=10MPa，邊界定產(chǎn)Q=15m3/d。求各點(diǎn)網(wǎng)格點(diǎn)在不同換為等距離x左邊的微分方程(2-11)進(jìn)行差分求解。方程(2p(t0)(1-7)(2)邊界：1)定產(chǎn)邊界：pQrrrw2即為所求的一維徑向單相流的數(shù)學(xué)模型。方程中的未知量為p(r,ke=10MPa，邊界定產(chǎn)Q=15m3/d。求各點(diǎn)網(wǎng)格點(diǎn)在不同換為等距離x左邊的微分方程(2-11)進(jìn)行差分求解。方程(2p(t0)(1-7)(2)邊界：1)定產(chǎn)邊界：pQrrrw2即為所求的一維徑向單相流的數(shù)學(xué)模型。方程中的未知量為p(r,ktt--1一維徑向單相流數(shù)學(xué)模型方程為(1-4)即為所求的一維徑向單相流的數(shù)學(xué)模型。方程中的未知量為，則上述方程可簡(jiǎn)化為：(rp)Cp(1-4)方程為(1-4)邊界條件下，方程(2-15)對(duì)應(yīng)于i=0到n的各個(gè)網(wǎng)格所構(gòu)成cm3/s；滲流微分方程(1-4)，則上述方程可簡(jiǎn)化為：(rp)Cp(1-4)方程為(1-4)邊界條件下，方程(2-15)對(duì)應(yīng)于i=0到n的各個(gè)網(wǎng)格所構(gòu)成cm3/s；滲流微分方程(1-4)與初始條件、邊界條件一起，幾乎重合。這也說明了生產(chǎn)之初壓力下降速度較快，但由于是定壓邊wwwf--pQei為：-可修-.--rrr于是數(shù)學(xué)模型(2-4)可轉(zhuǎn)換為：2(n則：i1iii1i式(2-15)即為一維徑向流時(shí)的差分方程-3MPa-1，原始?jí)毫=10MPa，最大模擬時(shí)間t=36pn=pe已知，因此只需要求第0到n-1個(gè)網(wǎng)格點(diǎn)的壓力。如上1w為：-可修-.--rrr于是數(shù)學(xué)模型(2-4)可轉(zhuǎn)換為：2(n則：i1iii1i式(2-15)即為一維徑向流時(shí)的差分方程-3MPa-1，原始?jí)毫=10MPa，最大模擬時(shí)間t=36pn=pe已知，因此只需要求第0到n-1個(gè)網(wǎng)格點(diǎn)的壓力。如上1w21w32wenww1--滲流微分方程(1-4)與初始條件、邊界條件一起，構(gòu)成了一維徑向單相流問壓力分布，以及井底流壓p或產(chǎn)量。2差分方程的建立劃分采用不等距網(wǎng)格，即井底附近網(wǎng)格劃分密一些，遠(yuǎn)離井底要疏一些。2iirr為把一維徑向坐標(biāo)r轉(zhuǎn)換為直角坐標(biāo)x，需要找到r與x的對(duì)應(yīng)關(guān)系。由式回一系列的結(jié)果，繪制網(wǎng)格劃分示意圖、各網(wǎng)格點(diǎn)在不同時(shí)刻壓力分劃分密一些，遠(yuǎn)離井底要疏一些。在此選取等比級(jí)數(shù)網(wǎng)格，即：a,r(2-2)這樣實(shí)現(xiàn)了井底附近網(wǎng)格小，而遠(yuǎn)離井底處網(wǎng)格壓大的附近壓力變化慢的特點(diǎn)，網(wǎng)格劃分采用不等距網(wǎng)格，即井底附近網(wǎng)格x為方程w回一系列的結(jié)果，繪制網(wǎng)格劃分示意圖、各網(wǎng)格點(diǎn)在不同時(shí)刻壓力分劃分密一些，遠(yuǎn)離井底要疏一些。在此選取等比級(jí)數(shù)網(wǎng)格，即：a,r(2-2)這樣實(shí)現(xiàn)了井底附近網(wǎng)格小，而遠(yuǎn)離井底處網(wǎng)格壓大的附近壓力變化慢的特點(diǎn)，網(wǎng)格劃分采用不等距網(wǎng)格，即井底附近網(wǎng)格x為方程wwwwwnwewrnrnw--令wwxwxwwwrr.附近壓力變化慢的特點(diǎn)，網(wǎng)格劃分采用不等距網(wǎng)格，即井底附近網(wǎng)格-5)可看出，r與x之間的對(duì)應(yīng)關(guān)系為：lnx于是：wrr.(緩，整個(gè)油藏逐步達(dá)到穩(wěn)態(tài)。通過以上圖像分析可知，得到的變化趨布圖和不同時(shí)刻各網(wǎng)格點(diǎn)的壓力沿徑向的分布圖。Solve函數(shù)主1r2x21pCp22附近壓力變化慢的特點(diǎn)，網(wǎng)格劃分采用不等距網(wǎng)格，即井底附近網(wǎng)格-5)可看出，r與x之間的對(duì)應(yīng)關(guān)系為：lnx于是：wrr.(緩，整個(gè)油藏逐步達(dá)到穩(wěn)態(tài)。通過以上圖像分析可知，得到的變化趨布圖和不同時(shí)刻各網(wǎng)格點(diǎn)的壓力沿徑向的分布圖。Solve函數(shù)主1r2x21pCp22wkp111n2112nCk--222C2pw令Mr2令布圖和不同時(shí)刻各網(wǎng)格點(diǎn)的壓力沿徑向的分布圖。Solve函數(shù)主越遠(yuǎn)），壓力下降幅度越小，下降的速度也越慢，在外邊界處壓力保r(2-2)這樣實(shí)現(xiàn)了井底附近網(wǎng)格小，而遠(yuǎn)離井底處網(wǎng)格壓大的kh2)定流壓邊界：p(r,t)p(t0)(1-9)式中，rQ布圖和不同時(shí)刻各網(wǎng)格點(diǎn)的壓力沿徑向的分布圖。Solve函數(shù)主越遠(yuǎn)），壓力下降幅度越小，下降的速度也越慢，在外邊界處壓力保r(2-2)這樣實(shí)現(xiàn)了井底附近網(wǎng)格小，而遠(yuǎn)離井底處網(wǎng)格壓大的kh2)定流壓邊界：p(r,t)p(t0)(1-9)式中，rQeepQ--3一維徑向單相流模擬事例3.1模擬條件與要求wepQpQ-5)邊界條件包括外邊界和邊界。相應(yīng)的外邊界條件如下：(1)p(t0)(1-7)(2)邊界：1)定產(chǎn)邊界：pQrrrw2其中，主程序Main主要作用是輸入地層、流體參數(shù)以及初始和邊找到r與x的對(duì)應(yīng)關(guān)系。由式(2-2)可得：.-可修-.--rppndpdpep0pdd-5)邊界條件包括外邊界和邊界。相應(yīng)的外邊界條件如下：(1)p(t0)(1-7)(2)邊界：1)定產(chǎn)邊界：pQrrrw2其中，主程序Main主要作用是輸入地層、流體參數(shù)以及初始和邊找到r與x的對(duì)應(yīng)關(guān)系。由式(2-2)可得：.-可修-.--rppndpdpep0pddppd112221pnpn2d21pppd--1ppi01nn4計(jì)算程序框圖一維徑向流程序框圖如圖2所示。p(t0)(1-7)(2)邊界：1)定產(chǎn)邊界：pQrrrw2p(t0)(1-7)(2)邊界：1)定產(chǎn)邊界：pQrrrw2上述過程，將不等距的徑向坐標(biāo)r轉(zhuǎn)換成了等距離的x坐標(biāo)，而且將界條件，設(shè)置與模擬時(shí)間相關(guān)的參數(shù)，通過調(diào)用Solve函數(shù)，返外邊界：1)封閉外邊界：rrre2)定壓外邊界：p(r,t)--圖2一維徑向流程序框圖5模擬結(jié)果分析e=10MPa，邊界定產(chǎn)Q=15m3/d。求各點(diǎn)網(wǎng)格點(diǎn)在不同3e=10MPa，邊界定產(chǎn)Q=15m3/d。求各點(diǎn)網(wǎng)格點(diǎn)在不同3中給定的條件，繪制的網(wǎng)格分布劃分示意圖如圖3所示。圖3網(wǎng)格件，該類邊界條件一般形式為：p(r,t)p下面主要構(gòu)建在上述圖1不等距r坐標(biāo)與等距x坐標(biāo)之間的轉(zhuǎn)換已知r，r和網(wǎng)格數(shù)n時(shí)--0-50-100-150-200-250PPM90增加a圖4各網(wǎng)格點(diǎn)在不同時(shí)刻壓力分布），劃分示意圖各網(wǎng)格節(jié)點(diǎn)在不同時(shí)刻的壓力分布如圖4所示：圖4各網(wǎng)劃分示意圖各網(wǎng)格節(jié)點(diǎn)在不同時(shí)刻的壓力分布如圖4所示：圖4各網(wǎng)勢(shì)符合一維徑向定壓邊界線填充圖和各網(wǎng)格點(diǎn)壓力沿徑向分布圖，分別如圖5和圖6所示。圖越遠(yuǎn)），壓力下降幅度越小，下降的速度也越慢，在外邊界處壓力保--0-100-200-2000-100-200-200M9.59000-100-200圖54個(gè)時(shí)間節(jié)點(diǎn)處壓力分布等值線填充圖圖64時(shí)間節(jié)點(diǎn)各網(wǎng)格點(diǎn)壓力沿徑向分布圖