資源描述：
《油藏?cái)?shù)值模擬大作業(yè)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、1差分方程組線性化二維空間油水兩相流體流動(dòng)方程組二維空間油水兩相流體流動(dòng)的差分方程組油相：水相：其中：令：上式兩邊同乘以，則油水兩相的差分方程為：油相：水相：令：上兩式可化為：差分方程線性化：方程右端項(xiàng)處理：油相：水相：利用IMPES方法，壓力取隱式形式，飽和度和毛管壓力取顯示形式。令進(jìn)一步可化簡(jiǎn)為：其中：油方程/Co2＋水方程/Cw2，得：（1）式中：如果網(wǎng)格塊上有井，需加上產(chǎn)量項(xiàng)，有兩種情況：1.定產(chǎn)量，直接在左端加上，注入為正，生產(chǎn)為負(fù);2.定井底流壓，處理如下：分別將產(chǎn)水量和產(chǎn)油量加到對(duì)應(yīng)的油、水方程中，即：2求解差分方程首先，對(duì)網(wǎng)格進(jìn)行編號(hào)
2、，見(jiàn)表2.1。表2.1網(wǎng)格編號(hào) 51015 16111621 4552 27121722263138465357593813182327323947545860491419242833404855 20252934414956 303542 364350 374451 其中，18號(hào)井是注水井，11號(hào)井和46號(hào)井為生產(chǎn)井。將網(wǎng)格編號(hào)后，就可計(jì)算出公式（1）的各個(gè)系數(shù)，并將其賦值到60*60的矩陣方程中，從而可計(jì)算出各網(wǎng)格的壓力變換量，最終求得不同生產(chǎn)時(shí)間的油水兩相壓力分布和各
3、相飽和度分布，以及生產(chǎn)井的生產(chǎn)數(shù)據(jù)。（其中矩陣方程的系數(shù)形狀見(jiàn)圖1）圖2.1系數(shù)矩陣形狀利用C語(yǔ)言編制程序?qū)崿F(xiàn)以上步驟，詳細(xì)程序見(jiàn)附錄。3求解結(jié)果（1）分別解出100d時(shí)的各網(wǎng)格的油相壓力及各網(wǎng)格的含水飽和度，見(jiàn)表3.1和表3.2。并根據(jù)表中數(shù)據(jù)分別繪制出地層壓力分布和含水飽和度分布圖，見(jiàn)圖3.1和圖3.2。表3.1100d時(shí)油相壓力分布 5737.545731.545739.20 5793.845738.575737.525842.105902.66 5873.965848.08 5746.585844.995942.44585
4、0.476056.926081.266012.145953.055826.945849.025954.825957.565849.455850.155951.325936.066059.816162.366153.015957.535950.125853.405955.855956.645851.815852.365954.855958.716161.446183.266173.966060.406056.155955.60 5961.686162.086189.776184.646062.236058.055956.47 620
5、4.126205.326104.27 6236.436106.216167.08 6257.256107.276167.54 表3.2100d時(shí)含水飽和度分布 0.46460.46240.4686 0.46720.47030.00000.46710.4739 0.45970.4654 0.46510.46630.46220.46490.46650.46830.46790.46450.00000.46710.46220.46170.46430.46220.46240.46560.46500.46
6、560.46520.46330.46150.46340.46310.46370.46320.46290.46340.46440.46540.46510.46330.46380.46290.4638 0.45940.46830.46390.46440.46410.46300.4632 0.46550.46420.4646 0.46310.46390.4640 0.46240.46340.4630 圖3.1300d時(shí)油相壓力分布圖3.2300d時(shí)含水飽和度分布（2）分別解出300d時(shí)的各網(wǎng)格的油相
7、壓力及各網(wǎng)格的含水飽和度，見(jiàn)表3.3和表3.4。并根據(jù)表中數(shù)據(jù)分別繪制出地層壓力分布和含水飽和度分布圖，見(jiàn)圖3.3和圖3.4。表3.3300d時(shí)油相壓力分布 5459.425450.945462.57 5466.245451.195436.055467.565484.22 5460.025480.20 5472.995470.785467.605480.375491.135498.935401.475488.285450.425481.755490.175494.195478.885482685481.085489.165495.
8、915400.815402.875494.875483.775488.195491.685494.31548
